Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,101,541 G→C 5.5% P28R (CCG→CGG)  yggN ← DUF2884 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,101,5410GC 5.5% 328.1 / 6.7 105P28R (CCG→CGG) yggNDUF2884 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (46/53);  new base C (4/2);  total (50/55)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 111 total aligned reads displayed.

TGATAATCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGT  >  NC_000913/3101394‑3101688
                                                                                                                                                   |                                                                                                                                                   
tGATAATCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCGCGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACgcggc                                                                                                                                                   >  1:949323/1‑150 (MQ=255)
tGATAATCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCAGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACgcggc                                                                                                                                                   >  1:1351350/1‑150 (MQ=255)
 gATAATCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGTATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGTTTGCCGTTTTCGCCCTTCACCTGCTCGGTTTGCGGGCTGACAATCGCATCGTCACgcggcg                                                                                                                                                  >  1:215564/1‑150 (MQ=255)
   tAATCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTg                                                                                                                                                <  2:214318/150‑1 (MQ=255)
    aaTCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGAGTGa                                                                                                                                               <  2:1537189/150‑1 (MQ=255)
            cccTGCTCGCGCGGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCCGCGTGACGATgtga                                                                                                                                       <  1:407129/149‑4 (MQ=255)
             ccTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACt                                                                                                                                      <  2:601213/150‑1 (MQ=255)
             ccTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCCGCGTGACGATgtgat                                                                                                                                      <  2:920414/150‑5 (MQ=255)
              cTGCTAGCTCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTg                                                                                                                                     <  1:219544/150‑1 (MQ=255)
                gcgcgcgCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGt                                                                                                                                   <  1:232781/147‑1 (MQ=255)
                gCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGt                                                                                                                                   <  1:1035335/150‑1 (MQ=255)
                gCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGt                                                                                                                                   >  1:478331/1‑150 (MQ=255)
                 cTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGGTATACATCACGTTGCCGTCTGGCGTGAGCACCAGATTGCCGTTTTCGCCCTTCACCTGCCCGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTa                                                                                                                                  <  1:1096998/150‑1 (MQ=255)
                   cgcgCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGt                                                                                                                                >  2:1362260/1‑150 (MQ=255)
                       cTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGc                                                                                                                            >  2:1441177/1‑150 (MQ=255)
                         gggCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGGCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGt                                                                                                                          <  2:1711499/150‑1 (MQ=255)
                              gcATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCg                                                                                                                     >  1:1357681/1‑150 (MQ=255)
                               cATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCgg                                                                                                                    <  1:199048/150‑1 (MQ=255)
                               cATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCgg                                                                                                                    >  2:1616137/1‑150 (MQ=255)
                                aTTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGt                                                                                                                   >  2:982863/1‑150 (MQ=255)
                                  tCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCa                                                                                                                 <  2:411335/150‑1 (MQ=255)
                                    aGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACAAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCAtt                                                                                                               >  1:862389/1‑150 (MQ=255)
                                     gggAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTg                                                                                                              >  1:691527/1‑150 (MQ=255)
                                      ggAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGc                                                                                                             <  1:735519/150‑1 (MQ=255)
                                      ggAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGc                                                                                                             <  2:1350128/150‑1 (MQ=255)
                                        aaTATTGCTTACCGTTATACATCACGTTGCCGGCTGGCGTGAGCACCAGAGTGCCGTTTTCGCCCTTCACCTGCCCGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCg                                                                                                           <  1:88974/150‑1 (MQ=255)
                                        aaTATTGCTAACCGTTATACATCACGTTGCCGTCTGGCGTGATCACGAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCg                                                                                                           <  2:349752/150‑1 (MQ=255)
                                           atTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCa                                                                                                        >  1:452739/1‑150 (MQ=255)
                                            ttGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCAc                                                                                                       <  1:1011985/150‑1 (MQ=255)
                                                 tACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGaaa                                                                                                  >  1:1351830/1‑150 (MQ=255)
                                                  aCCGTTATACATCACGTTGCCGGCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAg                                                                                                 <  2:1609864/150‑1 (MQ=255)
                                                   ccGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGt                                                                                                >  2:1685034/1‑150 (MQ=255)
                                                    cGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCCCGCGGCGTGACGCTGCACTGGTGGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGtg                                                                                               >  1:893666/1‑150 (MQ=255)
                                                    cGTTATACATCACGGTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACCTCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGtg                                                                                               <  2:55454/150‑1 (MQ=255)
                                                      ttATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGtgct                                                                                             <  1:341893/150‑1 (MQ=255)
                                                       tataCATCCCGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACCATCACATCGTCACGCGGCGGGACGCTGCACTGGTAGTCGGCGTGAGCGGGCATTGCCGTCACTGAAAGtgctg                                                                                            <  2:646445/150‑1 (MQ=255)
                                                         taCATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCTCGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTgcc                                                                                          >  2:1341763/1‑150 (MQ=255)
                                                         taCATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGGGACGCGGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTgcc                                                                                          <  1:992422/150‑1 (MQ=255)
                                                           cATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAAACACATCGTCACGCGGCGTGACGCTGCAATGGTAGTCGGCGTGAGCGGTCATTGCGGTCACTGAAATTGATgccgc                                                                                        >  2:1212579/1‑150 (MQ=255)
                                                            aTCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTgccgcc                                                                                       <  1:1700898/150‑1 (MQ=255)
                                                               aCGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCcagc                                                                                    <  2:1150985/150‑1 (MQ=255)
                                                                     ccGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCAtt                                                                              >  1:990673/1‑150 (MQ=255)
                                                                      cGTCTGGCGTGATCACCAGATTGGCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTTACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCAttt                                                                             <  2:1409082/150‑1 (MQ=255)
                                                                        tCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCTGGCGGCAGCActttc                                                                           <  2:238753/150‑2 (MQ=255)
                                                                           ggCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTgcgc                                                                        >  1:626029/1‑150 (MQ=255)
                                                                              gTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGcatc                                                                     >  1:1004888/1‑150 (MQ=255)
                                                                               tGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGcatca                                                                    >  2:1604348/1‑150 (MQ=255)
                                                                               tGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGcatca                                                                    >  2:602989/1‑150 (MQ=255)
                                                                                   cACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGt                                                                >  2:314390/1‑150 (MQ=255)
                                                                                   cACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGt                                                                >  1:1078945/1‑150 (MQ=255)
                                                                                   cACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGt                                                                >  1:1576424/1‑150 (MQ=255)
                                                                                   cACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGGAGTCGGCGTGCGCGGTCCTTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGt                                                                >  1:1073385/1‑150 (MQ=255)
                                                                                   cACCAGATTGCCGTTGTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGt                                                                <  1:443292/150‑1 (MQ=255)
                                                                                     ccAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCt                                                              >  1:650808/1‑150 (MQ=255)
                                                                                     ccAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGAAGTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCt                                                              >  2:1298264/1‑150 (MQ=255)
                                                                                     ccAGATTGCCGGATTCGCCCTTCACCTGCACTGTTTGCGCGCTGACAATCACATCGTCACGCGGCGTGACGCGGCACTGGTCGGCGGCGTGAGCGGTCATTGCCGTCACTGACAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCt                                                              <  2:296664/150‑1 (MQ=255)
                                                                                       agagTGCCGTTTGCGCCCTTCACCTGCTAGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCTTTGCCGTCACTTTAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTc                                                            <  2:1117351/150‑1 (MQ=255)
                                                                                        gATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTcc                                                           <  2:1728629/150‑1 (MQ=255)
                                                                                        gATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTcc                                                           >  1:1722147/1‑150 (MQ=255)
                                                                                           tGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTc                                                        <  2:352821/150‑1 (MQ=255)
                                                                                           tGCCGTTTTCGCCCGTGACCTGCTCGGGTTGCGGGCTGACAATCACATCGTCACGCGGCGGGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTc                                                        <  1:305357/150‑1 (MQ=255)
                                                                                            gCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACCTCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCCTCCTAGTCTTCCCGCa                                                       >  2:330534/1‑150 (MQ=255)
                                                                                            gCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCa                                                       >  1:269196/1‑150 (MQ=255)
                                                                                             ccGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCaa                                                      >  1:688911/1‑150 (MQ=255)
                                                                                                     gCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAAt                                              <  1:184285/150‑1 (MQ=255)
                                                                                                      ctcTTGACCTGCACGGTTTGCGGGCTGACATTCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGAGGTCATTGCCGTCACTGAAATTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATc                                             <  1:267491/148‑1 (MQ=255)
                                                                                                      cccTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATc                                             >  1:651809/1‑150 (MQ=255)
                                                                                                       ccTTCACCTGCACGGTTTGCGGGCTGACACTCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCa                                            <  1:1096900/150‑1 (MQ=255)
                                                                                                       ccTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCa                                            >  2:25962/1‑150 (MQ=255)
                                                                                                       ccTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCAATGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATcc                                            >  2:275973/1‑149 (MQ=255)
                                                                                                           cACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGcc                                        >  1:1092279/1‑150 (MQ=255)
                                                                                                             ccTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAg                                      >  1:1310483/1‑150 (MQ=255)
                                                                                                              cTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTATGCGCATCAGAATCTTCCCTCACGAAAAAATCCGGCCAGc                                     >  1:826679/1‑150 (MQ=255)
                                                                                                               tGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTTCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCa                                    >  1:836609/1‑150 (MQ=255)
                                                                                                                   cGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTg                                >  2:397877/1‑150 (MQ=255)
                                                                                                                    ggTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGc                               >  1:950242/1‑150 (MQ=255)
                                                                                                                     gTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCt                              <  2:1357681/150‑1 (MQ=255)
                                                                                                                       ttGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTgg                            <  2:1554405/150‑1 (MQ=255)
                                                                                                                           gggCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcctg                        <  2:1246061/150‑3 (MQ=255)
                                                                                                                           gggCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcctg                        <  1:898543/150‑3 (MQ=255)
                                                                                                                              cTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAAGCAGGCCAGCAAATGCTGGcctgatt                     <  1:83465/150‑6 (MQ=255)
                                                                                                                               tGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcctgattt                    <  2:1004888/150‑7 (MQ=255)
                                                                                                                                gACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAtt                   >  1:1484185/1‑150 (MQ=255)
                                                                                                                                  cAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                 >  1:713600/1‑150 (MQ=255)
                                                                                                                                    aTCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATa               >  2:1435929/1‑150 (MQ=255)
                                                                                                                                       gcagcctCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTAATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGGATCATAGTCTTCCCTCAAGAAAAAATCAGGCCACCATTTGCTGGCACCAGATTGATAACa            >  2:170130/7‑150 (MQ=255)
                                                                                                                                       gcagcGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATATTGCCCATCAAATAATTCCCTCAAGAAAAAAACAGGCCCGCATTTGATGGCCCCAGATGGATAACa            >  2:866958/5‑150 (MQ=255)
                                                                                                                                        cagcGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAATAATCAGGCCAGCATTTGCTGGCCCCAGATTGATATCaa           >  1:1492540/4‑150 (MQ=255)
                                                                                                                                        cagcGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaa           >  1:1485709/4‑150 (MQ=255)
                                                                                                                                         ctCGGCACGCGGCCTGACGCTGCCCTGGTAGCCGGCGTGAGCGGGCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGCCCAGCATTTGCTGGCCCCAGATTGATAACaac          <  2:862389/149‑2 (MQ=255)
                                                                                                                                          tCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg         <  1:422834/150‑1 (MQ=255)
                                                                                                                                          tCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg         <  1:1162341/150‑1 (MQ=255)
                                                                                                                                          tCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg         <  1:1362260/150‑1 (MQ=255)
                                                                                                                                           cGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGt        <  2:1576424/150‑1 (MQ=255)
                                                                                                                                            gTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTa       <  1:330534/150‑1 (MQ=255)
                                                                                                                                              cACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcg     <  2:691527/150‑1 (MQ=255)
                                                                                                                                              cACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTACCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcg     <  1:80494/150‑1 (MQ=255)
                                                                                                                                               aCGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATAGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgc    <  1:1095158/150‑1 (MQ=255)
                                                                                                                                                cgcgGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgcg   >  1:611649/1‑150 (MQ=255)
                                                                                                                                                 gcggcgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGACCCAGATTGATAACAAAGTACGCGt  >  1:457079/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGATAATCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGT  >  NC_000913/3101394‑3101688

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