Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,736,916 A→G 6.2% intergenic (+33/‑238) bamD → / → raiA BamABCDE complex OM biogenesis lipoprotein/cold shock protein associated with 30S ribosomal subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,736,9160AG 6.2% 459.1 / 11.3 141intergenic (+33/‑238)bamD/raiABamABCDE complex OM biogenesis lipoprotein/cold shock protein associated with 30S ribosomal subunit
Reads supporting (aligned to +/- strand):  ref base A (74/58);  new base G (3/6);  total (77/64)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 141 total aligned reads displayed.

CCAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACAAAAA  >  NC_000913/2736769‑2737063
                                                                                                                                                    |                                                                                                                                                   
ccAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAg                                                                                                                                                    >  1:1074592/1‑149 (MQ=255)
 cAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGc                                                                                                                                                   >  1:1170430/1‑149 (MQ=255)
  aGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCt                                                                                                                                                  >  3:1065320/1‑149 (MQ=255)
  aGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCt                                                                                                                                                  <  2:1470834/149‑1 (MQ=255)
  aGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCt                                                                                                                                                  >  2:1414294/1‑149 (MQ=255)
    gCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGc                                                                                                                                                >  1:1290569/1‑149 (MQ=255)
        cgcgTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCAAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCTCAAGGGCTGCCGtt                                                                                                                                            <  1:76829/149‑1 (MQ=255)
         gcgTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGttt                                                                                                                                           >  3:625732/1‑149 (MQ=255)
         gcgTGATGAGCTGCGGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACGTAACAGAAACCTGAAACAC‑AAAACGGCAGCTCAAGGGCTGCCGttt                                                                                                                                           <  4:966302/149‑1 (MQ=255)
           gTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGttttt                                                                                                                                         >  4:535136/1‑149 (MQ=255)
            tGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGtttttt                                                                                                                                        <  2:820445/149‑1 (MQ=255)
            tGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGtttttt                                                                                                                                        <  3:1107205/149‑1 (MQ=255)
                gcgcTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTAtt                                                                                                                                    <  1:1053038/149‑1 (MQ=255)
                  gcTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAAACGGCAGCTCAAGGGCTGCCGttttgtgtttc                                                                                                                                   <  2:1487947/149‑8 (MQ=38)
                    tGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGt                                                                                                                                <  4:1008357/149‑1 (MQ=255)
                     gCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTc                                                                                                                               <  1:1414294/149‑1 (MQ=255)
                     gCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTc                                                                                                                               <  1:815405/149‑1 (MQ=255)
                        gCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGt                                                                                                                            >  2:626524/1‑149 (MQ=255)
                          tGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTtg                                                                                                                          >  2:1519929/1‑149 (MQ=255)
                           taTGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTtgt                                                                                                                         <  3:136952/148‑1 (MQ=255)
                              ggAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGaa                                                                                                                      >  1:425920/1‑149 (MQ=255)
                               gAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGt                                                                                                                                >  4:250430/1‑138 (MQ=255)
                               gAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGt                                                                                                                                <  3:250430/138‑1 (MQ=255)
                                    tGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGaa                                                                                                                <  3:1201031/149‑1 (MQ=255)
                                    tGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGaa                                                                                                                <  4:1422138/149‑1 (MQ=255)
                                       aTACCGTCAGATGCAGATTAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATCCATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGaa                                                                                                                >  2:1176411/1‑146 (MQ=255)
                                       aTACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGaa                                                                                                                <  1:1176411/146‑1 (MQ=255)
                                       aTACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCg                                                                                                             >  2:501873/1‑149 (MQ=255)
                                       aTACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCg                                                                                                             >  2:303305/1‑149 (MQ=255)
                                        tACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGa                                                                                                            >  3:520515/1‑149 (MQ=255)
                                         aCCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGAt                                                                                                           <  2:750122/149‑1 (MQ=255)
                                          ccGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGAtt                                                                                                          <  4:247323/149‑1 (MQ=255)
                                               agatgcagatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTc                                                                                                     <  4:1291049/149‑1 (MQ=255)
                                               agatgcagatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTc                                                                                                     <  3:265857/149‑1 (MQ=255)
                                                  tgcagatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGt                                                                                                                                <  1:736445/119‑1 (MQ=255)
                                                   gcagatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTa                                                                                                 >  1:64026/1‑149 (MQ=255)
                                                   gcagatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTa                                                                                                 >  1:940888/1‑149 (MQ=255)
                                                     agatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGCCATCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGCTGTGATACTGAAGCGATTTAGTCGCTATc                                                                                               >  4:548513/1‑149 (MQ=255)
                                                      gatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCg                                                                                              >  2:743651/1‑149 (MQ=255)
                                                        tgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGAt                                                                                            <  2:940888/149‑1 (MQ=255)
                                                         gAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGAtc                                                                                           >  4:899995/1‑149 (MQ=255)
                                                               gcAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCATTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCtcatca                                                                                     >  4:1247297/1‑149 (MQ=255)
                                                                   gCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAata                                                                                 <  3:557984/149‑1 (MQ=255)
                                                                   gCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAata                                                                                 >  2:940373/1‑149 (MQ=255)
                                                                   gCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAata                                                                                 <  2:64026/149‑1 (MQ=255)
                                                                    cTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAatat                                                                                <  2:425920/149‑1 (MQ=255)
                                                                           aGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTc                                                                                                                                   <  3:309821/91‑1 (MQ=255)
                                                                           aGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTc                                                                                                                                   >  4:309821/1‑91 (MQ=255)
                                                                              aGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCttt                                                                      >  1:137356/1‑149 (MQ=255)
                                                                               gCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑GAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTg                                                                                                                  >  2:841785/1‑104 (MQ=255)
                                                                               gCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑GAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTg                                                                                                                  <  1:841785/104‑1 (MQ=255)
                                                                               gCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTg                                                                     <  4:1021147/149‑1 (MQ=255)
                                                                               gCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTg                                                                     <  3:743643/149‑1 (MQ=255)
                                                                                 gAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTgag                                                                   >  1:506361/1‑149 (MQ=255)
                                                                                  aaaaTCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTgaga                                                                  >  2:511175/1‑149 (MQ=255)
                                                                                   aaaTCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATgg                                                                              >  4:554353/1‑136 (MQ=255)
                                                                                   aaaTCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATgg                                                                              <  3:554353/136‑1 (MQ=255)
                                                                                         tcGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCt                                                           >  2:1268932/1‑149 (MQ=255)
                                                                                         tcGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCt                                                           >  1:54581/1‑149 (MQ=255)
                                                                                            ccgcAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCaa                                                        >  3:1025063/1‑149 (MQ=255)
                                                                                            ccgcAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCaa                                                        <  1:809568/149‑1 (MQ=255)
                                                                                            ccgcAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCaa                                                        <  1:911444/149‑1 (MQ=255)
                                                                                              gcAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGt                                                      >  3:353870/1‑149 (MQ=255)
                                                                                                        gcaATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACAt                                            <  2:1170430/149‑1 (MQ=255)
                                                                                                               aTAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGAtt                                                                                                          >  1:284472/1‑80 (MQ=255)
                                                                                                               aTAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACa                                             <  1:1209455/141‑1 (MQ=255)
                                                                                                               aTAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACa                                             >  2:1209455/1‑141 (MQ=255)
                                                                                                                 aaCAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCaaa                                                                                   >  3:915002/1‑101 (MQ=255)
                                                                                                                 aaCAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGTTATTCCTCAAGTAAAAAAACATCTATTCCta                                   >  2:952097/1‑148 (MQ=255)
                                                                                                                 aaCAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTg                                   >  3:248393/1‑149 (MQ=255)
                                                                                                                 aaCAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTg                                   >  4:14501/1‑149 (MQ=255)
                                                                                                                  aCAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAg                                                                                                                                                    <  4:44347/35‑1 (MQ=255)
                                                                                                                  aCAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAg                                                                                                                                                    >  3:44347/1‑35 (MQ=255)
                                                                                                                  aCAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGc                                  >  2:284223/1‑149 (MQ=255)
                                                                                                                    aGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTc                                      >  2:167438/1‑143 (MQ=255)
                                                                                                                    aGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGa                                >  2:649274/1‑149 (MQ=255)
                                                                                                                     gAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGAt                               <  1:1555606/149‑1 (MQ=255)
                                                                                                                       aaCCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGAttt                             >  1:213061/1‑149 (MQ=255)
                                                                                                                        aCCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGCTTtc                            >  4:1415215/1‑149 (MQ=255)
                                                                                                                        aCCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTtc                            <  2:213061/149‑1 (MQ=255)
                                                                                                                         ccTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTtct                           <  4:919980/149‑1 (MQ=255)
                                                                                                                           tGAAACAC‑AAAACGGCAGCCCTTGGGCTGCCGATTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTca                         >  4:313098/1‑149 (MQ=255)
                                                                                                                           tGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTca                         <  2:832323/149‑1 (MQ=255)
                                                                                                                               acac‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaa                     <  1:649274/149‑1 (MQ=255)
                                                                                                                                 ac‑aAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaaaa                   >  3:202115/1‑149 (MQ=255)
                                                                                                                                  c‑aAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaaaaa                  <  2:1074592/149‑1 (MQ=255)
                                                                                                                                  c‑aAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaaaaa                  <  4:1164457/149‑1 (MQ=255)
                                                                                                                                    aaaaCGGCAGCTCAAGGGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACATAAAAg                 >  3:897753/1‑149 (MQ=255)
                                                                                                                                    aaaaCGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAg                 >  2:1409536/1‑149 (MQ=255)
                                                                                                                                    aaaaCGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAg                 >  1:872705/1‑149 (MQ=255)
                                                                                                                                        cGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTc             <  1:1569318/149‑1 (MQ=255)
                                                                                                                                          gCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGt           >  4:389243/1‑149 (MQ=255)
                                                                                                                                          gCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGt           >  4:1173232/1‑149 (MQ=255)
                                                                                                                                            aGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTg         >  2:956449/1‑149 (MQ=255)
                                                                                                                                               ccTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGAt               <  4:60656/140‑1 (MQ=255)
                                                                                                                                               ccTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACa      >  3:62410/1‑149 (MQ=255)
                                                                                                                                               ccTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACa      >  1:191737/1‑149 (MQ=255)
                                                                                                                                                cTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACaa     >  2:1007417/1‑149 (MQ=255)
                                                                                                                                                 ttGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACaaa    >  4:707561/1‑149 (MQ=255)
                                                                                                                                                   gAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACaaaaa  <  4:1025063/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACAAAAA  >  NC_000913/2736769‑2737063

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑