Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,718,784 G→T 5.1% W674L (TGG→TTG) ‡ PP_4171 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,718,7840GT 5.1% 790.9 / 21.6 233W674L (TGG→TTG) ‡PP_4171hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (115/106);  new base T (8/4);  total (123/110)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 243 total aligned reads displayed.

GTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCTTTGTC  >  NC_002947/4718635‑4718933
                                                                                                                                                     |                                                                                                                                                     
gTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGAGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCtgtg                                                                                                                                                       >  2:5103016/1‑150 (MQ=255)
  ggATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGgc                                                                                                                                                     <  1:553771/150‑1 (MQ=255)
      gcgcGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGgcagca                                                                                                                                                 <  2:867001/150‑1 (MQ=255)
      gcgcGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGgcagca                                                                                                                                                 <  2:2679057/150‑1 (MQ=255)
      gcgcGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGgcagca                                                                                                                                                 <  2:2505741/150‑1 (MQ=255)
        gcgAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCagcagc                                                                                                                                               >  2:5874862/1‑150 (MQ=255)
          gAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAgcgc                                                                                                                                             >  1:3212711/1‑150 (MQ=255)
          gAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAgcgc                                                                                                                                             >  2:4453186/1‑150 (MQ=255)
           aaTGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCGACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGCTGCCACAGCGcc                                                                                                                                            <  2:5608614/150‑1 (MQ=255)
            aTGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCt                                                                                                                                           <  1:5013032/150‑1 (MQ=255)
             tGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTg                                                                                                                                          <  1:1541274/150‑1 (MQ=255)
              ggggTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTgg                                                                                                                                         <  1:4862094/150‑1 (MQ=255)
               gggTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGa                                                                                                                                        <  2:5372096/150‑1 (MQ=255)
                      ggCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGaa                                                                                                                                 <  1:3290480/150‑1 (MQ=255)
                      ggCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGaa                                                                                                                                 <  2:6722585/150‑1 (MQ=255)
                       gCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAAc                                                                                                                                <  2:2832019/150‑1 (MQ=255)
                        cTGGTCGAGCGCGCTGTTCTATTCCGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAAcc                                                                                                                               <  2:6277387/150‑1 (MQ=255)
                         tGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAAccc                                                                                                                              >  2:3196212/1‑150 (MQ=255)
                          ggTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCACGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCa                                                                                                                             <  2:4390484/150‑1 (MQ=255)
                             cGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCt                                                                                                                          <  1:2172380/150‑1 (MQ=255)
                              gAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTc                                                                                                                         >  1:1397117/1‑150 (MQ=255)
                                      tGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGc                                                                                                                 >  2:2651002/1‑150 (MQ=255)
                                      tGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGc                                                                                                                 >  2:3872519/1‑150 (MQ=255)
                                      tGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGc                                                                                                                 <  1:4643463/150‑1 (MQ=255)
                                      tGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGc                                                                                                                 >  2:3232475/1‑150 (MQ=255)
                                       gTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCa                                                                                                                <  1:982265/150‑1 (MQ=255)
                                       gTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCa                                                                                                                <  2:4237162/150‑1 (MQ=255)
                                           tATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAAc                                                                                                            >  1:5709362/1‑150 (MQ=255)
                                           tATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAAc                                                                                                            >  2:6837741/1‑150 (MQ=255)
                                            aTTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAAcc                                                                                                           <  2:63927/150‑1 (MQ=255)
                                             ttACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCg                                                                                                          <  1:6401864/150‑1 (MQ=255)
                                             ttACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCg                                                                                                          >  1:4195302/1‑150 (MQ=255)
                                                cGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTa                                                                                                       >  2:1960989/1‑150 (MQ=255)
                                                cGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTa                                                                                                       <  2:432572/150‑1 (MQ=255)
                                                 ggTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTAc                                                                                                      >  2:4873322/1‑150 (MQ=255)
                                                 ggTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTGTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTAc                                                                                                      >  2:1200551/1‑150 (MQ=255)
                                                           cccTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTc                                                                                            <  1:5919605/150‑1 (MQ=255)
                                                            ccTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAgcgc                                                                                                                                             >  2:2964592/1‑100 (MQ=255)
                                                            ccTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGAAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCa                                                                                           >  2:3058112/1‑150 (MQ=255)
                                                             cTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCAc                                                                                          <  1:1232830/150‑1 (MQ=255)
                                                             cTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCAc                                                                                          <  2:1074244/150‑1 (MQ=255)
                                                              tCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCAcc                                                                                         <  2:2892250/150‑1 (MQ=255)
                                                              tCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCAcc                                                                                         >  2:1622051/1‑150 (MQ=255)
                                                                  ccAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGgc                                                                                                                                                     >  2:2097630/1‑86 (MQ=255)
                                                                   cAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtg                                                                                    >  2:3767777/1‑150 (MQ=255)
                                                                   cAGTACCGCTACGAGCGTCACGACATGCTACTGGCCAAGCGTTTTCGGGTGTATGGCAGCAAACAGGAGCAGGCGGAGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCCGCACTGCAGAACCGCTACGACAGCCGACACCGCCtc                                                                                    >  2:4791258/1‑149 (MQ=255)
                                                                    aGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtga                                                                                   <  1:6606974/150‑1 (MQ=255)
                                                                        ccGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgt                                                                               >  2:1163343/1‑150 (MQ=255)
                                                                          gCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgtga                                                                             <  1:6080811/150‑1 (MQ=255)
                                                                          gCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgtga                                                                             <  1:2878806/150‑1 (MQ=255)
                                                                           cTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgtgag                                                                            >  2:6515473/1‑150 (MQ=255)
                                                                                    tCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCt                                                                   >  2:6053742/1‑150 (MQ=255)
                                                                                        gACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAg                                                               <  1:5332913/150‑1 (MQ=255)
                                                                                          ccTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGtt                                                             <  1:267609/150‑1 (MQ=255)
                                                                                           cTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTc                                                            >  1:4445715/1‑150 (MQ=255)
                                                                                            tGCGCCTGGCCAAGCGTTTTCGGGTGGATGGCAGCAAACTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCt                                                           >  2:563936/1‑150 (MQ=255)
                                                                                               gcCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGAt                                                        >  1:6206028/1‑150 (MQ=255)
                                                                                               gcCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGAt                                                        >  2:352906/1‑150 (MQ=255)
                                                                                                ccTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATg                                                       >  1:1273852/1‑150 (MQ=255)
                                                                                                 cTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATgg                                                      >  2:1576856/1‑150 (MQ=255)
                                                                                                   ggCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGcc                                                    <  1:1470694/150‑1 (MQ=255)
                                                                                                    gCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGccc                                                   >  1:1046442/1‑150 (MQ=255)
                                                                                                       aaGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGCTGCCACAgcgc                                                                                                                                             <  1:4537624/57‑1 (MQ=37)
                                                                                                        aGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctgct                                               >  1:483553/1‑150 (MQ=255)
                                                                                                        aGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctgct                                               >  2:6234501/1‑150 (MQ=255)
                                                                                                        aGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctgct                                               >  2:3074415/1‑150 (MQ=255)
                                                                                                             tttCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCtgctgctgc                                          <  1:2417512/150‑1 (MQ=255)
                                                                                                                cGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGcc                                       <  1:3720373/150‑1 (MQ=255)
                                                                                                                 gggTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGc                                                                                       >  2:3314554/1‑101 (MQ=255)
                                                                                                                  ggTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgtgagc                                                                           <  2:1250193/112‑1 (MQ=255)
                                                                                                                  ggTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTg                                     <  1:3842774/150‑1 (MQ=255)
                                                                                                                   gtgtATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGt                                    <  1:2598394/150‑1 (MQ=255)
                                                                                                                   gtgtATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGt                                    <  2:40536/150‑1 (MQ=255)
                                                                                                                    tgtATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGtt                                   <  1:6640133/150‑1 (MQ=255)
                                                                                                                     gtATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTtg                                  >  2:4642921/1‑150 (MQ=255)
                                                                                                                          gcagcaACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTgctgct                             >  1:3145049/1‑150 (MQ=255)
                                                                                                                          gcagcaACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTgctgct                             >  2:2255489/1‑150 (MQ=255)
                                                                                                                              caACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGc                         <  2:3424434/150‑1 (MQ=255)
                                                                                                                               aaCCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGcc                        >  1:5881642/1‑150 (MQ=255)
                                                                                                                                aCCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCt                       <  1:2724412/150‑1 (MQ=255)
                                                                                                                                   tGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCtgtg                    >  1:5223563/1‑150 (MQ=255)
                                                                                                                                   tGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCtgtg                    <  2:4913806/150‑1 (MQ=255)
                                                                                                                                   tGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCtgtg                    >  2:1706497/1‑150 (MQ=255)
                                                                                                                                      aGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGcc                 <  2:6764172/150‑1 (MQ=255)
                                                                                                                                        cTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGt               <  1:1831306/150‑1 (MQ=255)
                                                                                                                                         tccaggcgCTGCTGCCACAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGtt              >  1:3618995/5‑150 (MQ=255)
                                                                                                                                           gcggcgCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTgg            >  1:4777206/1‑150 (MQ=255)
                                                                                                                                           gcggcgCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTgg            <  2:1518558/150‑1 (MQ=255)
                                                                                                                                            aggcgCTGCTGCCACAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTggg           >  1:1343060/2‑150 (MQ=255)
                                                                                                                                             ggcgCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGcc                                                    >  1:2090651/1‑108 (MQ=255)
                                                                                                                                             ggcgCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTtgctgc                              >  1:5784700/1‑130 (MQ=255)
                                                                                                                                             ggcgCTGCTGCCACAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTgggg          >  2:2331359/1‑150 (MQ=255)
                                                                                                                                               cgcTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGcc        <  1:3856826/150‑1 (MQ=255)
                                                                                                                                                gcTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCt       <  1:5088195/150‑1 (MQ=255)
                                                                                                                                                gcTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCt       >  2:6466874/1‑150 (MQ=255)
                                                                                                                                                gcTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCt       <  1:1321178/150‑1 (MQ=255)
                                                                                                                                                gcTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCt       <  2:2412701/150‑1 (MQ=255)
                                                                                                                                                  tgtgGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCttt     <  1:1664220/150‑1 (MQ=255)
                                                                                                                                                  tgtgGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCttt     <  2:5161702/150‑1 (MQ=255)
                                                                                                                                                     ggCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCTTTGTc  <  2:186551/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCTTTGTC  >  NC_002947/4718635‑4718933

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: