Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,736,914 T→A 6.5% intergenic (+31/‑240) 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,9140TA 6.5% 583.1 / 23.9 170intergenic (+31/‑240)bamD/raiABamABCDE complex OM biogenesis lipoprotein/cold shock protein associated with 30S ribosomal subunit
Reads supporting (aligned to +/- strand):  ref base T (74/83);  new base A (2/9);  total (76/94)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

Only 100 of 171 total aligned reads displayed.

CCCAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACAAA  >  NC_000913/2736768‑2737061
                                                                                                                                                  |                                                                                                                                                   
cccAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGa                                                                                                                                                   <  2:582791/149‑1 (MQ=255)
cccAGGCTACGCGTGATGCGCCGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGACACCTGAATCACAAAACGGCAGCCCTTGa                                                                                                                                                   <  2:1570170/149‑1 (MQ=255)
 ccAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAg                                                                                                                                                  >  3:1139451/1‑149 (MQ=255)
 ccAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAg                                                                                                                                                  >  4:391348/1‑149 (MQ=255)
 ccAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAg                                                                                                                                                  >  2:1617934/1‑149 (MQ=255)
 ccAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAg                                                                                                                                                  >  1:1415830/1‑149 (MQ=255)
 ccAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAg                                                                                                                                                  >  3:502587/1‑149 (MQ=255)
     gCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGc                                                                                                                                              >  2:2239/1‑149 (MQ=255)
         cgcgTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCTCAAGGGCTGCCGtt                                                                                                                                          <  1:940144/149‑1 (MQ=255)
         cgcgTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCTCAAGGGCTGCCGtt                                                                                                                                          <  1:750194/149‑1 (MQ=255)
          gcgTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCTCAAGGGCTGCCGttt                                                                                                                                         <  1:1052059/149‑1 (MQ=255)
          gcgTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCTCAAGGGCTGCCGttt                                                                                                                                         <  1:229006/149‑1 (MQ=255)
          gcgTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGttt                                                                                                                                         >  3:1551176/1‑149 (MQ=255)
           cgTGATGCGCTTCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGtttt                                                                                                                                        <  3:913428/149‑1 (MQ=255)
           cgTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGtttt                                                                                                                                        <  3:386708/149‑1 (MQ=255)
            gtaaTGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGttttt                                                                                                                                       <  1:480461/146‑1 (MQ=255)
            gTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGttttt                                                                                                                                       <  3:797349/149‑1 (MQ=255)
            gTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGttttt                                                                                                                                       <  3:268029/149‑1 (MQ=255)
            gTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGttttt                                                                                                                                       >  3:99927/1‑149 (MQ=255)
            gTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGttttt                                                                                                                                       <  3:295977/149‑1 (MQ=255)
             tGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCTCAAGGGCTGCCGttttgt                                                                                                                                      <  3:697327/149‑3 (MQ=255)
             tGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGtttttt                                                                                                                                      <  3:391348/149‑1 (MQ=255)
             tGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGtttttt                                                                                                                                      <  3:371290/149‑1 (MQ=255)
                tGCGCTGCCGCTGATTGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTAt                                                                                                                                   <  1:931586/149‑1 (MQ=255)
                      gCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTc                                                                                                                             <  4:386335/149‑1 (MQ=255)
                          cTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTCGAGCTGCCGTTTTTTTATTCTGTCAGtt                                                                                                                         <  3:1447711/149‑1 (MQ=255)
                          cTGATGGAAAATGCATAACGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGtt                                                                                                                         <  4:1663780/149‑1 (MQ=255)
                             aTGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCTCAAGGGCTGCCGtttt                                                                                                                                        >  3:1402465/1‑131 (MQ=255)
                             aTGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCTCAAGGGCTGCCGtttt                                                                                                                                        <  4:1402465/131‑1 (MQ=255)
                              tGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGa                                                                                                                     >  2:85597/1‑149 (MQ=255)
                                gAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGaaa                                                                                                                   >  1:378958/1‑149 (MQ=255)
                                 aaaaTGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAAc                                                                                                                  >  3:1296671/1‑149 (MQ=255)
                                   aaTGCATCCCGTCAGATGCAGATGCATGCGCACGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACATCAGCAATAGATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTg                                                                                                                >  1:1116854/1‑149 (MQ=255)
                                   aaTGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTg                                                                                                                <  3:1252910/149‑1 (MQ=255)
                                     tGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGaa                                                                                                              <  2:1116854/149‑1 (MQ=255)
                                      gCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTc                                                                                                                                 <  1:1066314/129‑1 (MQ=255)
                                        aTACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCg                                                                                                           >  4:835523/1‑149 (MQ=255)
                                         tACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGATACTGAAGCGAt                                                                                                         <  3:1459031/149‑1 (MQ=255)
                                          aCCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGTCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGAt                                                                                                         <  4:1237192/149‑1 (MQ=255)
                                          aCCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGAAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGAt                                                                                                         <  3:958115/149‑1 (MQ=255)
                                          aCCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCAAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGAt                                                                                                         <  3:238745/149‑1 (MQ=255)
                                           ccGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGttt                                                                                                        <  1:635597/149‑1 (MQ=255)
                                           ccGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGAtt                                                                                                        <  2:420400/149‑1 (MQ=255)
                                            cGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGAttt                                                                                                       >  3:1314935/1‑149 (MQ=255)
                                              tCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAg                                                                                                     <  2:660078/149‑1 (MQ=255)
                                                agatgcagatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTc                                                                                                   >  2:1652085/1‑149 (MQ=255)
                                                 gatgcagatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCg                                                                                                  >  2:1244283/1‑149 (MQ=255)
                                                  atgcagatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGc                                                                                                 <  3:835523/149‑1 (MQ=255)
                                                      agatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATc                                                                                             >  1:1027647/1‑149 (MQ=255)
                                                       gatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCg                                                                                            >  2:1301293/1‑149 (MQ=255)
                                                        atgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGa                                                                                           >  4:1079013/1‑149 (MQ=255)
                                                        atgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGa                                                                                           <  4:363314/149‑1 (MQ=255)
                                                         tgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGACTAAATCGCttcagtt                                                                                          <  3:587954/149‑7 (MQ=255)
                                                           aaTGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGAtct                                                                                        <  4:1529687/149‑1 (MQ=255)
                                                             tGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCtca                                                                                      <  1:1193503/149‑1 (MQ=255)
                                                                gcAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCtcatca                                                                                   >  2:750227/1‑149 (MQ=255)
                                                                      tGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATg                                                                             >  1:544799/1‑149 (MQ=255)
                                                                      tGAAAAAGTAGCGAAAATCATCGCCACAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCt                                                                                                                                <  2:117930/98‑1 (MQ=255)
                                                                       gAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATgg                                                                            <  3:1055251/149‑1 (MQ=255)
                                                                             gTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCt                                                                      >  4:1519060/1‑149 (MQ=255)
                                                                              tAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCtt                                                                     <  1:1177916/149‑1 (MQ=255)
                                                                               aGCGAAAATCATCGCCGCAAACAGCAGCAATACATAGCAGAAACCTGAAACGCAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCttt                                                                    <  4:1139451/149‑1 (MQ=255)
                                                                               aGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCttt                                                                    <  1:1586688/149‑1 (MQ=255)
                                                                               aGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCttt                                                                    >  3:1019092/1‑149 (MQ=255)
                                                                                gCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGAt                                                                                                         >  4:1304350/1‑111 (MQ=255)
                                                                                 cGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTga                                                                  <  3:1432669/149‑1 (MQ=255)
                                                                                 cGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTga                                                                  <  1:1326653/149‑1 (MQ=255)
                                                                                 cGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTga                                                                  <  4:1359687/149‑1 (MQ=255)
                                                                                          tcGCCGCAAACAGCAGCAATACATAACAGAAACCTTAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCt                                                         <  1:1457729/149‑1 (MQ=255)
                                                                                           cgccgcAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTc                                                        <  3:1079013/149‑1 (MQ=255)
                                                                                                     agcagcAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCaa                                                      <  2:1137740/141‑1 (MQ=255)
                                                                                                       cagcaATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTc                                                           <  2:313803/134‑1 (MQ=255)
                                                                                                       cagcaATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTc                                                           >  1:313803/1‑134 (MQ=255)
                                                                                                          caATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGa                                                                                                               <  1:1015615/79‑1 (MQ=255)
                                                                                                          caATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGa                                                                                                               >  2:1015615/1‑79 (MQ=255)
                                                                                                          caATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACAtc                                         <  4:543117/149‑1 (MQ=255)
                                                                                                                aTAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTcc                                   >  2:413357/1‑149 (MQ=255)
                                                                                                                    cAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCg                               >  3:313326/1‑149 (MQ=255)
                                                                                                                     aGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGa                              >  3:1662404/1‑149 (MQ=255)
                                                                                                                     aGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGa                              >  4:1432732/1‑149 (MQ=255)
                                                                                                                     aGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGa                              >  1:544705/1‑149 (MQ=255)
                                                                                                                      gAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGAt                             >  3:1358293/1‑149 (MQ=255)
                                                                                                                          ccTGAAACACAAAACGGCAGCCCTTGAGATGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTtct                         <  1:351951/149‑1 (MQ=255)
                                                                                                                           cTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTtctc                        <  1:1199756/149‑1 (MQ=255)
                                                                                                                            tGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTca                       >  3:565896/1‑149 (MQ=255)
                                                                                                                            tGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTca                       >  4:699336/1‑149 (MQ=255)
                                                                                                                             gAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTcac                      >  2:1379076/1‑149 (MQ=255)
                                                                                                                                   caAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaaaaa                >  3:32794/1‑149 (MQ=255)
                                                                                                                                    aaaaCGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAatat                                                                              <  4:601163/86‑1 (MQ=255)
                                                                                                                                    aaaaCGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAatat                                                                              >  3:601163/1‑86 (MQ=255)
                                                                                                                                    aaaaCGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaaaaa                >  1:872718/1‑148 (MQ=255)
                                                                                                                                    aaaaCGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAg               >  1:1453195/1‑149 (MQ=255)
                                                                                                                                    aaaaCGGCAGACCTGGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaaaaa                <  2:872718/148‑1 (MQ=255)
                                                                                                                                     aaaCGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGa              <  1:1617934/149‑1 (MQ=255)
                                                                                                                                      aaCGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGAt             <  4:1019092/149‑1 (MQ=255)
                                                                                                                                           cAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaaaaa                >  2:1029728/1‑141 (MQ=255)
                                                                                                                                           cAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGtt        <  4:1136459/149‑1 (MQ=255)
                                                                                                                                             gCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGa      <  1:1652085/149‑1 (MQ=255)
                                                                                                                                                cTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACaa   >  4:499964/1‑149 (MQ=255)
                                                                                                                                                 ttGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACaaa  >  2:1382167/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCCAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACAAA  >  NC_000913/2736768‑2737061

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