Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,706,980 T→C 5.1% Y140H (TAT→CAT) ‡ PP_3267 → Clp protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,706,9800TC 5.1% 1277.0 / 33.1 393Y140H (TAT→CAT) ‡PP_3267Clp protease
Reads supporting (aligned to +/- strand):  ref base T (172/201);  new base C (12/8);  total (184/209)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 397 total aligned reads displayed.

GCAATATCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATCGCGTCA  >  NC_002947/3706831‑3707129
                                                                                                                                                     |                                                                                                                                                     
 cAATATCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTa                                                                                                                                                      <  1:6044623/150‑1 (MQ=255)
  aaTATCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTAt                                                                                                                                                     <  1:7614786/150‑1 (MQ=255)
  aaTATCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTAt                                                                                                                                                     >  1:8271296/1‑150 (MQ=255)
      tCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcgc                                                                                                                                                 >  1:8808972/1‑150 (MQ=255)
       ccTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCa                                                                                                                                                <  1:3235985/150‑1 (MQ=255)
       ccTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCa                                                                                                                                                >  2:5518853/1‑150 (MQ=255)
       ccTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCCGGATTACGACCTGCGCCTGTATGCGCa                                                                                                                                                >  2:2024452/1‑150 (MQ=255)
         tctTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGa                                                                                                                                              <  1:6305012/150‑1 (MQ=255)
           ttCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCATACAGGCGCAGGTc                                                                                                                                            <  2:8700629/150‑1 (MQ=255)
              cgggCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTGTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTg                                                                                                                                         <  2:70282/147‑1 (MQ=255)
              cTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTg                                                                                                                                         <  2:1317506/150‑1 (MQ=255)
              cTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTg                                                                                                                                         >  2:8031359/1‑150 (MQ=255)
               ttGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGATCTGCGCATACAGGCGCAGGTCGtaa                                                                                                                                        <  1:1994653/150‑3 (MQ=255)
               ttGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTgg                                                                                                                                        >  2:4536213/1‑150 (MQ=255)
                 gCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGATCTGCGCATACAGGCGCAGGTCGtaatc                                                                                                                                      <  1:8534872/150‑5 (MQ=255)
                     gagaGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAgc                                                                                                                                  <  2:4560551/150‑1 (MQ=255)
                            aCGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAgcgcctttc                                                                                                                           <  1:2483859/150‑6 (MQ=255)
                            aCGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGaa                                                                                                                           >  1:7816426/1‑150 (MQ=255)
                              gCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAg                                                                                                                         <  1:3091411/150‑1 (MQ=255)
                               cACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAgg                                                                                                                        >  2:8782212/1‑150 (MQ=255)
                               cACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAgg                                                                                                                        >  2:77087/1‑150 (MQ=255)
                                aCGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGgc                                                                                                                       >  2:6070611/1‑150 (MQ=255)
                                aCGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGgc                                                                                                                       >  1:489836/1‑150 (MQ=255)
                                aCGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGgc                                                                                                                       >  1:6618596/1‑150 (MQ=255)
                                aCGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGgc                                                                                                                       <  2:6691941/150‑1 (MQ=255)
                                    cTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCgg                                                                                                                   >  2:7473031/1‑150 (MQ=255)
                                        cAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGa                                                                                                               <  2:4342017/150‑1 (MQ=255)
                                               aGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTg                                                                                                        <  1:4492151/150‑1 (MQ=255)
                                               aGCAAGTTTCTTGTCCCCCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTg                                                                                                        <  2:6163771/150‑1 (MQ=255)
                                                gCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTgg                                                                                                       <  1:4224036/150‑1 (MQ=255)
                                                gCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTgg                                                                                                       >  2:2405729/1‑150 (MQ=255)
                                                 cAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGa                                                                                                      <  2:4975562/150‑1 (MQ=255)
                                                 cAAGTTTCTTTGCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGa                                                                                                      <  2:7302460/150‑1 (MQ=255)
                                                        cTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACcg                                                                                               >  1:7994063/1‑150 (MQ=255)
                                                          tttCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACcgcg                                                                                             >  2:8533135/1‑150 (MQ=255)
                                                          tttCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTCGGATGTACcgcg                                                                                             >  1:2109860/1‑150 (MQ=255)
                                                               aCCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATc                                                                                        <  1:3348520/150‑1 (MQ=255)
                                                               aCCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATc                                                                                        <  1:8147441/150‑1 (MQ=255)
                                                               aCCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATc                                                                                        >  2:6466322/1‑150 (MQ=255)
                                                               aCCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATc                                                                                        <  2:4162893/150‑1 (MQ=255)
                                                                  cGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTaa                                                                                     <  1:7961463/150‑1 (MQ=255)
                                                                    tttCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATg                                                                                   <  1:2284374/150‑1 (MQ=255)
                                                                      tCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAg                                                                                                              <  1:8047556/121‑1 (MQ=255)
                                                                      tCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGc                                                                                 <  2:7557688/150‑1 (MQ=255)
                                                                      tCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGc                                                                                 <  1:5832623/150‑1 (MQ=255)
                                                                        aTTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCt                                                                               >  2:2890906/1‑150 (MQ=255)
                                                                        aTTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCt                                                                               >  2:1547382/1‑150 (MQ=255)
                                                                         ttGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTa                                                                              >  1:5570662/1‑150 (MQ=255)
                                                                         ttGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGGTATCTAATGGCCTa                                                                              >  1:1001759/1‑150 (MQ=255)
                                                                              ccTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCaa                                                                         <  1:3461247/150‑1 (MQ=255)
                                                                              ccTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCaa                                                                         <  1:6881502/150‑1 (MQ=255)
                                                                              ccTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCaa                                                                         <  1:3778685/150‑1 (MQ=255)
                                                                              ccTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCaa                                                                         <  2:3863653/150‑1 (MQ=255)
                                                                              ccTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCaa                                                                         <  1:1334812/150‑1 (MQ=255)
                                                                               cTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAg                                                                        <  2:4504766/150‑1 (MQ=255)
                                                                                       gTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTgg                                                                <  2:3471268/150‑1 (MQ=255)
                                                                                             tGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGcgccgc                                                          <  1:1069179/150‑1 (MQ=255)
                                                                                             tGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGcgccgc                                                          >  2:6288203/1‑150 (MQ=255)
                                                                                               gACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCgg                                                        <  1:1785673/150‑1 (MQ=255)
                                                                                               gACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCgg                                                        >  2:8327304/1‑150 (MQ=255)
                                                                                               gACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCgg                                                        <  2:480579/150‑1 (MQ=255)
                                                                                                aCCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCggg                                                       >  2:2722418/1‑150 (MQ=255)
                                                                                                 ccATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGa                                                      <  1:2886719/150‑1 (MQ=255)
                                                                                                 ccATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGa                                                      >  2:3928389/1‑150 (MQ=255)
                                                                                                       gcgGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCt                                                <  1:1428999/150‑1 (MQ=255)
                                                                                                       gcgGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCt                                                <  2:2129762/150‑1 (MQ=255)
                                                                                                        cgGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTg                                               <  2:7400034/150‑1 (MQ=255)
                                                                                                         ggATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTgg                                              <  2:2550684/150‑1 (MQ=255)
                                                                                                           aTGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGaa                                            <  1:5850163/150‑1 (MQ=255)
                                                                                                              tCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGgc                                                                                                                       <  1:5510137/72‑1 (MQ=255)
                                                                                                              tCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCCGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAAc                                         <  1:7357543/150‑1 (MQ=255)
                                                                                                                cGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACgg                                       <  1:3647715/150‑1 (MQ=255)
                                                                                                                 gAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACggg                                      <  1:7130187/150‑1 (MQ=255)
                                                                                                                  aGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGga                                     <  1:4115000/150‑1 (MQ=255)
                                                                                                                  aGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGga                                     <  1:5450135/150‑1 (MQ=255)
                                                                                                                  aGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGga                                     <  2:7913229/150‑1 (MQ=255)
                                                                                                                     aTGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACg                                        <  2:8438886/144‑1 (MQ=255)
                                                                                                                       gCGATGAGCCTGGATTACGACCTGCGCCTATATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATc                                <  1:5160835/150‑1 (MQ=255)
                                                                                                                        cGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATcc                               >  2:8185721/1‑150 (MQ=255)
                                                                                                                        cGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATcc                               >  1:6726460/1‑150 (MQ=255)
                                                                                                                         gATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCa                              <  2:5550456/150‑1 (MQ=255)
                                                                                                                              gCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGaaa                         <  2:4628607/150‑1 (MQ=255)
                                                                                                                                cTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATc                       <  2:5269962/150‑1 (MQ=255)
                                                                                                                                cTGGATTACGACCTGCGCCTGTAGGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATc                       <  2:3059547/150‑1 (MQ=255)
                                                                                                                                  ggATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcga                     <  1:7250950/150‑1 (MQ=255)
                                                                                                                                   gATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGTTCCATGAAATcgac                    <  2:6879991/150‑1 (MQ=255)
                                                                                                                                   gATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgac                    <  1:4537041/150‑1 (MQ=255)
                                                                                                                                   gATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgac                    <  2:1813739/150‑1 (MQ=255)
                                                                                                                                    aTTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgacg                   >  2:1038609/1‑150 (MQ=255)
                                                                                                                                     ttACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgacga                  <  2:7429447/150‑1 (MQ=255)
                                                                                                                                     ttACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgacga                  <  1:7487995/150‑1 (MQ=255)
                                                                                                                                     ccaCGATCTGCGCATACAGGCGCAGGTCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGTTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgacga                  >  1:2166461/3‑150 (MQ=255)
                                                                                                                                       aCGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACgaga                <  1:8682637/150‑1 (MQ=255)
                                                                                                                                         gACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATg              >  1:336852/1‑150 (MQ=255)
                                                                                                                                          aCCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGc             <  2:3252951/150‑1 (MQ=255)
                                                                                                                                            cTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCg           >  1:2699069/1‑150 (MQ=255)
                                                                                                                                                gcCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATcg       >  2:193296/1‑150 (MQ=255)
                                                                                                                                                  cTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATcgcg     >  1:3980148/1‑150 (MQ=255)
                                                                                                                                                    gTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATCGCGTc   >  1:6999342/1‑150 (MQ=255)
                                                                                                                                                    gTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATCGCGTc   >  2:3349724/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAATATCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATCGCGTCA  >  NC_002947/3706831‑3707129

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: