Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,378,351 C→G 5.2% S395T (AGC→ACC)  dnaK ← chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,378,3510CG 5.2% 522.7 / 13.5 155S395T (AGC→ACC) dnaKchaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (67/80);  new base G (4/4);  total (71/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 160 total aligned reads displayed.

GCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCG  >  NC_002947/5378202‑5378495
                                                                                                                                                     |                                                                                                                                                
gCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGc                                                                                                                                                  >  2:59882/1‑150 (MQ=255)
gCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGc                                                                                                                                                  >  1:1525341/1‑150 (MQ=255)
gCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGc                                                                                                                                                  >  2:1864030/1‑150 (MQ=255)
gCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGc                                                                                                                                                  >  2:4348437/1‑150 (MQ=255)
 cTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCt                                                                                                                                                 <  2:5069861/150‑1 (MQ=255)
   cACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGa                                                                                                                                               <  1:5005821/150‑1 (MQ=255)
    aCCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGAc                                                                                                                                              <  2:799997/150‑1 (MQ=255)
      catGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGACAACACCTGCGACTTCTTGGTCGGGATGCAGCTCTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCTATAACCAGTGTCAGCGGGCTGACGt                                                                                                                                            >  2:3026450/3‑150 (MQ=255)
       ttGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTc                                                                                                                                           <  2:2556631/150‑1 (MQ=255)
       gtGCAGCACGTGACTGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTc                                                                                                                                           <  2:2844680/149‑1 (MQ=255)
                gTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcagcagcag                                                                                                                                  >  1:4121886/1‑150 (MQ=255)
                   aaTGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGAc                                                                                                                               >  1:4514407/1‑150 (MQ=255)
                   aaTGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGAc                                                                                                                               <  2:4298642/150‑1 (MQ=255)
                     tGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGTGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGt                                                                                                                             <  1:3821940/150‑1 (MQ=255)
                       gTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCt                                                                                                                           >  1:4482157/1‑150 (MQ=255)
                       gTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCt                                                                                                                           <  2:410824/150‑1 (MQ=255)
                       gTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCt                                                                                                                           >  2:3326148/1‑150 (MQ=255)
                       gTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAGGCGGGTGAACACGCCACCGAGCGTTTCAATGCCGAGCGAAAGCGGGGTCACGTCGAGCAGCAGAACGTcc                                                                                                                           <  1:3937228/150‑2 (MQ=255)
                         cACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCttt                                                                                                                         <  2:2311095/150‑1 (MQ=255)
                         cACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCttt                                                                                                                         <  1:290077/150‑1 (MQ=255)
                         cACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCttt                                                                                                                         <  2:1226198/150‑1 (MQ=255)
                          aCGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTc                                                                                                                        <  1:2643237/150‑1 (MQ=255)
                          aCGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTc                                                                                                                        >  1:966681/1‑150 (MQ=255)
                          aCGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTc                                                                                                                        >  1:597905/1‑150 (MQ=255)
                          aCGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTc                                                                                                                        >  1:138938/1‑150 (MQ=255)
                            ggCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCAc                                                                                                                      <  2:4581226/150‑1 (MQ=255)
                             gcgcTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACg                                                                                                                     <  2:48050/150‑1 (MQ=255)
                             gcgcTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACg                                                                                                                     >  2:2027438/1‑150 (MQ=255)
                             gcgcTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACg                                                                                                                     <  2:3980200/150‑1 (MQ=255)
                               gcTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTc                                                                                                                   >  2:3186892/1‑150 (MQ=255)
                               gcTCGGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCGTTCACGTc                                                                                                                   >  1:106887/1‑150 (MQ=255)
                                 tctGGTTGTCGTCGGCAGTCGATAACACCTGCGACGTCTTGGTCGGGATGGTGGTGTACTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGATTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGacac                                                                                                                 >  2:4188759/1‑150 (MQ=255)
                                 tctGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCAc                                                                                                                 >  1:737254/1‑150 (MQ=255)
                                     attGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAGGCGGGTGAACACGCCACCGAGCGTTTCAATGCCGAGCGAAAGCGGGGTCACGTCGAGCAGCAGAACGTCCTTGACGTCGccctg                                                                                                             >  1:1168354/2‑147 (MQ=255)
                                      ttGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAGGCGGGTGAACACGCCACCGAGCGTTTCAATGCCGAGCGAAAGCGGGGTCACGTCGAGCAGCAGAACGTCCTTGACGTCGccctgc                                                                                                            >  2:2678208/1‑146 (MQ=255)
                                       tGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAGGCGGGTGAACACGCCACCGAGCGTTTCAATGCCGAGCGAAAGCGGGGTCACGTCGAGCAGCAGAACGTCCTTGACGTCGccctgca                                                                                                           <  1:4066852/150‑6 (MQ=255)
                                       tGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCa                                                                                                           <  2:347614/150‑1 (MQ=255)
                                         tcgtcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTTTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCATGTCACCGGCCAGt                                                                                                         >  2:2381026/1‑150 (MQ=255)
                                          cgtcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCCGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTa                                                                                                        <  2:3025280/150‑1 (MQ=255)
                                            tcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTAcc                                                                                                      <  1:2199144/150‑1 (MQ=255)
                                            tcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTAcc                                                                                                      <  1:2869276/150‑1 (MQ=255)
                                                     gagaACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTgg                                                                                             >  1:3770570/1‑150 (MQ=255)
                                                     gagaACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTgg                                                                                             >  2:1118788/1‑150 (MQ=255)
                                                       gaACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGAt                                                                                           <  2:3236576/150‑1 (MQ=255)
                                                            ccTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGgcag                                                                                      <  2:1566880/150‑1 (MQ=255)
                                                             cTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGgcagc                                                                                     >  1:1844019/1‑150 (MQ=255)
                                                                 gACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCAccc                                                                                 <  1:2930302/150‑1 (MQ=255)
                                                                   cttcttGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCAt                                                                               <  1:4032677/150‑1 (MQ=255)
                                                                    ttcttGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATg                                                                              <  2:4424048/150‑1 (MQ=255)
                                                                      gttGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGc                                                                            <  1:4296584/149‑1 (MQ=255)
                                                                       ttGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCa                                                                           >  2:1888084/1‑150 (MQ=255)
                                                                         ggTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAAc                                                                         >  2:2720840/1‑150 (MQ=255)
                                                                           tCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACgg                                                                       >  2:1409580/1‑150 (MQ=255)
                                                                           tCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACgg                                                                       <  1:4676415/150‑1 (MQ=255)
                                                                              ggATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCtt                                                                    >  1:214415/1‑150 (MQ=255)
                                                                              ggATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCtt                                                                    <  1:232066/150‑1 (MQ=255)
                                                                                  ggtggtGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGGTTCGAGACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTc                                                                <  2:1491682/150‑1 (MQ=255)
                                                                                    tggtgTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCg                                                              <  1:4610888/150‑1 (MQ=255)
                                                                                        gtTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCGACGGCTTCGTCCGGGtt                                                          >  1:626489/1‑150 (MQ=255)
                                                                                           cttctCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGAc                                                       >  1:530248/1‑150 (MQ=255)
                                                                                            ttctCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGgca                                                                                       >  1:1945466/1‑117 (MQ=255)
                                                                                            ttctCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACg                                                      <  2:4113857/150‑1 (MQ=255)
                                                                                              ctcGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTc                                                    <  2:1247760/150‑1 (MQ=255)
                                                                                                cGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCt                                                  <  2:266141/150‑1 (MQ=255)
                                                                                                 gATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCtt                                                 >  2:5098158/1‑150 (MQ=255)
                                                                                                   tCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTgc                                               >  2:1044428/1‑150 (MQ=255)
                                                                                                    cAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTgcg                                              <  1:607731/150‑1 (MQ=255)
                                                                                                     aGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTgcgc                                             <  1:237251/150‑1 (MQ=255)
                                                                                                         cgGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGctt                                         <  1:323349/150‑1 (MQ=255)
                                                                                                           gTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGTCTTCGTCCGGGTTGACGTCCTTGCGCGcttct                                       >  1:3931811/1‑150 (MQ=255)
                                                                                                              atcaCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGgca                                                                                       >  1:1602590/1‑99 (MQ=255)
                                                                                                              atcaCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTa                                    >  2:3527158/1‑150 (MQ=255)
                                                                                                              atcaCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTa                                    >  1:1615552/1‑150 (MQ=255)
                                                                                                               tcaCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTAc                                   <  2:4380889/150‑1 (MQ=255)
                                                                                                               tcaCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTAc                                   >  1:919305/1‑150 (MQ=255)
                                                                                                                caCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTAcc                                  >  1:2033246/1‑150 (MQ=255)
                                                                                                                   gccgccCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaa                               <  2:4597939/150‑1 (MQ=255)
                                                                                                                   gccgccCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaa                               <  1:215661/150‑1 (MQ=255)
                                                                                                                   gccgccCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaa                               <  2:3456544/150‑1 (MQ=255)
                                                                                                                       cccATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCCGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGaagaag                           >  2:2581274/1‑150 (MQ=255)
                                                                                                                       cccATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGaagaag                           >  1:4047339/1‑150 (MQ=255)
                                                                                                                          aTGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCg                        >  1:3510738/1‑150 (MQ=255)
                                                                                                                          aTGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCg                        >  1:5012284/1‑150 (MQ=255)
                                                                                                                          aTGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCg                        >  1:1598672/1‑150 (MQ=255)
                                                                                                                               ttCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCAc                   >  1:4343552/1‑150 (MQ=255)
                                                                                                                                 cGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCActt                 >  2:3816138/1‑150 (MQ=255)
                                                                                                                                  gATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCActtc                <  2:2138916/150‑1 (MQ=255)
                                                                                                                                  gATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCActtc                >  1:118201/1‑150 (MQ=255)
                                                                                                                                  gATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCActtc                >  2:2399187/1‑150 (MQ=255)
                                                                                                                                  gATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCActtc                >  1:4905235/1‑150 (MQ=255)
                                                                                                                                  gATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCACACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCActtc                >  1:1006742/1‑150 (MQ=255)
                                                                                                                                     aCCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACttcttt             <  1:2137254/150‑1 (MQ=255)
                                                                                                                                     aCCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACttcttt             <  1:59882/150‑1 (MQ=255)
                                                                                                                                       ccAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACttctttct           >  2:1705464/1‑150 (MQ=255)
                                                                                                                                       ccAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACttctttct           <  1:3026450/150‑1 (MQ=255)
                                                                                                                                       ccAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACAGAAGAAGTCGGCCACttctttct           >  1:3503339/1‑150 (MQ=255)
                                                                                                                                             tCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCa     <  2:1615552/150‑1 (MQ=255)
                                                                                                                                                gCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCg  >  1:1674610/1‑150 (MQ=255)
                                                                                                                                                gCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCg  >  2:1126129/1‑150 (MQ=255)
                                                                                                                                                gCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCg  >  2:622213/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCG  >  NC_002947/5378202‑5378495

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