Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,444,193 A→C 7.6% T28P (ACC→CCC)  tamB → translocation and assembly module subunit TamB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,444,1930AC 7.6% 356.0 / 16.5 117T28P (ACC→CCC) tamBtranslocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base A (55/53);  new base C (6/3);  total (61/56)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 121 total aligned reads displayed.

TTTGCCGTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCT  >  NC_000913/4444044‑4444341
                                                                                                                                                     |                                                                                                                                                    
tttGCCGTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGca                                                                                                                                                      <  1:272053/150‑1 (MQ=255)
   gCCGTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGcacca                                                                                                                                                   <  2:674764/150‑1 (MQ=255)
    ccGTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGcaccac                                                                                                                                                  >  1:1475178/1‑150 (MQ=255)
       ttccGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAg                                                                                                                                               <  2:1499488/148‑1 (MQ=255)
         ccGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCg                                                                                                                                             >  1:438950/1‑150 (MQ=255)
           ggTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGTGCCCACCAgaaac                                                                                                                                           <  1:524873/150‑5 (MQ=255)
           ggTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGTGCCCACCAgaaac                                                                                                                                           <  1:527580/150‑5 (MQ=255)
            ggcgcgGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGcc                                                                                                                                          <  2:1361225/148‑1 (MQ=255)
            gTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGcc                                                                                                                                          <  2:835651/150‑1 (MQ=255)
              cgcgGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTg                                                                                                                                        <  1:496293/150‑1 (MQ=255)
              cgcgGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTg                                                                                                                                        >  1:1391762/1‑150 (MQ=255)
                cgGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCa                                                                                                                                      <  1:1505176/150‑1 (MQ=255)
                   aTAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGc                                                                                                                                           <  1:400909/142‑1 (MQ=255)
                   aTAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGc                                                                                                                                           >  2:400909/1‑142 (MQ=255)
                    tAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTg                                                                                                                                  >  2:674936/1‑150 (MQ=255)
                    tAAAGACGAACACGGGTTACAGGTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTg                                                                                                                                  <  1:949705/150‑1 (MQ=255)
                     aaaGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTgg                                                                                                                                 <  1:1422218/150‑1 (MQ=255)
                       aGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTa                                                                                                                               <  2:484610/150‑1 (MQ=255)
                       aGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTa                                                                                                                               <  2:694368/150‑1 (MQ=255)
                       aGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTa                                                                                                                               >  1:1094866/1‑150 (MQ=255)
                       aGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTa                                                                                                                               >  1:484610/1‑150 (MQ=255)
                        gACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCCCCAGCGGCCTGCATCTGGTAt                                                                                                                              >  2:1277268/1‑150 (MQ=255)
                        gACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTAt                                                                                                                              >  2:211688/1‑150 (MQ=255)
                             acacGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTaaa                                                                                                                         >  1:151405/1‑150 (MQ=255)
                              cacGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAg                                                                                                                        >  1:354308/1‑150 (MQ=255)
                                cGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAgcg                                                                                                                      <  1:211688/150‑1 (MQ=255)
                                  ggTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAgcggc                                                                                                                    <  1:503508/150‑1 (MQ=255)
                                   gTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAgcggcg                                                                                                                   <  2:463910/150‑1 (MQ=255)
                                      aCAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGAt                                                                                                                >  1:1238513/1‑150 (MQ=255)
                                       cAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATc                                                                                                               <  1:966778/150‑1 (MQ=255)
                                         gTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGc                                                                                                             <  1:1065048/150‑1 (MQ=255)
                                         gTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGc                                                                                                             <  1:16489/150‑1 (MQ=255)
                                            ttACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTGAAGCGGCGGATCGCTgg                                                                                                          >  1:385945/1‑150 (MQ=255)
                                              aCATCGGTCGGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGt                                                                                                        >  2:316887/1‑150 (MQ=255)
                                               cATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTg                                                                                                       >  2:485181/1‑150 (MQ=255)
                                                 tCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTg                                                                                                           <  1:297594/144‑1 (MQ=255)
                                                 tCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTg                                                                                                           >  2:297594/1‑144 (MQ=255)
                                                 tCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGcc                                                                                                     >  1:361187/1‑150 (MQ=255)
                                                 tCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGcc                                                                                                     >  1:1265395/1‑150 (MQ=255)
                                                         gggCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTgg                                                                                             <  2:12806/150‑1 (MQ=255)
                                                          ggCCAGAATTATGAGTTTATGGAAAAACATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCCCCCGCGGCCTGCATCTGGGATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGa                                                                                            >  2:1502303/1‑150 (MQ=255)
                                                          ggCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGa                                                                                            <  2:159764/150‑1 (MQ=255)
                                                          ggCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGa                                                                                            >  2:841945/1‑150 (MQ=255)
                                                           gCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGat                                                                                           >  2:1102795/1‑150 (MQ=255)
                                                            ccAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGata                                                                                          <  2:930727/150‑1 (MQ=255)
                                                              aGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATAtt                                                                                        >  1:1090359/1‑150 (MQ=255)
                                                                    aTGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAg                                                                                  >  2:64031/1‑150 (MQ=255)
                                                                          ttATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCAcc                                                                            <  1:1480664/150‑1 (MQ=255)
                                                                           tATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACcg                                                                           >  1:427352/1‑150 (MQ=255)
                                                                            aTGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACcgg                                                                          >  2:651337/1‑150 (MQ=255)
                                                                              ggAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGat                                                                                           <  1:1276607/131‑1 (MQ=255)
                                                                              ggAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGat                                                                                           >  2:1276607/1‑131 (MQ=255)
                                                                               gAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACcggcgg                                                                       >  1:952752/1‑150 (MQ=255)
                                                                                aaaaaaaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCggcggc                                                                      >  1:791517/1‑150 (MQ=255)
                                                                                    aaaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGgc                                                                  >  2:1431409/1‑150 (MQ=255)
                                                                                     actCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTGCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGgcg                                                                 <  1:588037/148‑1 (MQ=255)
                                                                                     aaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGgcg                                                                 <  1:651337/150‑1 (MQ=255)
                                                                                          gCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAtc                                                            >  1:101062/1‑150 (MQ=255)
                                                                                          gCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAtc                                                            <  2:1037492/150‑1 (MQ=255)
                                                                                           ccTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTGCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGGTATTGGCAAGGTCACCGGCGGCTGGCGCGAtct                                                           <  1:778228/150‑1 (MQ=255)
                                                                                            cTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAtctc                                                          >  1:221281/1‑150 (MQ=255)
                                                                                                   tGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGt                                                   >  1:340861/1‑150 (MQ=255)
                                                                                                    ggTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTc                                                  >  2:436583/1‑150 (MQ=255)
                                                                                                       tatcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGa                                               >  1:110099/1‑150 (MQ=255)
                                                                                                        atcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAc                                              <  2:151405/150‑1 (MQ=255)
                                                                                                         tcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcg                                             <  1:769317/150‑1 (MQ=255)
                                                                                                          cgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgt                                            <  1:496462/150‑1 (MQ=255)
                                                                                                           gttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCCCCTTGTCTGAcgtt                                           >  2:1005124/1‑150 (MQ=255)
                                                                                                             tatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAtc                                                            <  2:1156373/131‑1 (MQ=255)
                                                                                                             tatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAtc                                                            >  1:1156373/1‑131 (MQ=255)
                                                                                                               tcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgttcgtt                                       <  2:1067044/150‑1 (MQ=255)
                                                                                                                cTTACTGGTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTACAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTa                                      <  2:1135505/150‑1 (MQ=255)
                                                                                                                  tACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATg                                    <  2:1485104/150‑1 (MQ=255)
                                                                                                                     tGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGagc                                 >  2:1452027/1‑150 (MQ=255)
                                                                                                                        tGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGagcagc                              >  2:517359/1‑150 (MQ=255)
                                                                                                                          cTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCa                            >  2:1091419/1‑150 (MQ=255)
                                                                                                                            gggATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAgg                          >  1:88004/1‑150 (MQ=255)
                                                                                                                              gATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCg                        >  2:123848/1‑150 (MQ=255)
                                                                                                                                  ggTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGc                    >  1:1317199/1‑150 (MQ=255)
                                                                                                                                  ggTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGc                    <  1:1497867/150‑1 (MQ=255)
                                                                                                                                   gTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCg                   >  2:28528/1‑150 (MQ=255)
                                                                                                                                   gTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCg                   >  2:1459501/1‑150 (MQ=255)
                                                                                                                                    tGGCGTTTCTGGTGGGCCCCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATAGTGGCAAGGTCACCGGCGGCTGGCGCGCTCTCACCTTGTCGGACGTTCGGTATGAGCAGCCAGGCGTGGCgg                  >  1:601393/1‑150 (MQ=255)
                                                                                                                                    tGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCgg                  <  1:316887/150‑1 (MQ=255)
                                                                                                                                    tGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCgg                  <  2:1090359/150‑1 (MQ=255)
                                                                                                                                      gCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTa                <  1:64031/150‑1 (MQ=255)
                                                                                                                                      gCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTa                <  2:427352/150‑1 (MQ=255)
                                                                                                                                      gCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTa                <  1:841945/150‑1 (MQ=255)
                                                                                                                                      aggccgcTGGTGGTGCCCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTa                >  2:131281/7‑150 (MQ=255)
                                                                                                                                       ggccgcTGGTGGTGCCCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTaa               >  1:812652/6‑150 (MQ=255)
                                                                                                                                       cGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTaa               >  2:160021/1‑150 (MQ=255)
                                                                                                                                       cGTATCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTaa               >  1:1258873/1‑150 (MQ=255)
                                                                                                                                        gccgcTGGTGGTGCCCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTaaa              <  1:131281/146‑1 (MQ=255)
                                                                                                                                         ccgcTGGTGGTGCCCACCAGCGGCCTGCATCTGCTATTTAAAGCGGCGCATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTaaaa             >  2:1153065/4‑150 (MQ=255)
                                                                                                                                          cgcTGGTGGTGCCCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAg            >  2:1407640/3‑150 (MQ=255)
                                                                                                                                             tggtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCgg         <  2:54375/150‑1 (MQ=255)
                                                                                                                                               gtggGCACCACCAGCGGCCTGCATCTGGTATTTGAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCAGCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGc       <  2:385945/150‑1 (MQ=255)
                                                                                                                                               gtggGCACCACCAGCGGCCTGCATCTGGTATTTGAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGAGGTTTTATATATATAAGCTTTTCTGACTTTCGTTATGAGAAGCAATGCGTTTCGGTAAAAGAggga       <  2:380373/150‑2 (MQ=255)
                                                                                                                                                  ggCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAAt    >  2:899813/1‑150 (MQ=255)
                                                                                                                                                    caccacCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCt  >  1:1416870/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTTGCCGTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCT  >  NC_000913/4444044‑4444341

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: