Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,289,857 A→T 5.6% N375I (AAT→ATT)  PP_1126 → hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,289,8570AT 5.6% 651.7 / 19.5 198N375I (AAT→ATT) PP_1126hydrolase
Reads supporting (aligned to +/- strand):  ref base A (101/86);  new base T (6/5);  total (107/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

Only 100 of 202 total aligned reads displayed.

GCACCCGTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCG  >  NC_002947/1289708‑1290004
                                                                                                                                                     |                                                                                                                                                    
gCACCCGTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCaa                                                                                                                                                      >  1:4828616/1‑150 (MQ=255)
gCACCCGTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCaa                                                                                                                                                      >  1:5014408/1‑150 (MQ=255)
 cACCCGTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAAt                                                                                                                                                     >  1:4773959/1‑150 (MQ=255)
  aCCCGTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATg                                                                                                                                                    >  2:3308095/1‑150 (MQ=255)
     cGTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCt                                                                                                                                                >  1:1816802/1‑150 (MQ=255)
      gTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCtc                                                                                                                                               >  1:4714677/1‑150 (MQ=255)
         cTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCtcgtc                                                                                                                                            <  1:5402181/150‑1 (MQ=255)
          tATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGAGGCATTGCCAtcgtcg                                                                                                                                           >  2:409400/1‑150 (MQ=255)
           aTTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGc                                                                                                                                          >  1:429963/1‑150 (MQ=255)
            ttCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCa                                                                                                                                         >  1:4585922/1‑150 (MQ=255)
            ttCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCa                                                                                                                                         >  1:4969973/1‑150 (MQ=255)
            ttCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCa                                                                                                                                         >  1:2232661/1‑150 (MQ=255)
             tCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCaa                                                                                                                                        <  1:2090229/150‑1 (MQ=255)
              cccTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAg                                                                                                                                       <  1:1673855/150‑1 (MQ=255)
               ccTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAgg                                                                                                                                      <  2:3139824/150‑1 (MQ=255)
               ccTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAgg                                                                                                                                      <  2:1096134/150‑1 (MQ=255)
               ccTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGAGGCATTGCCAtcgtcgttgtt                                                                                                                                      <  1:445558/150‑6 (MQ=255)
                cTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAggg                                                                                                                                     >  1:2377507/1‑150 (MQ=255)
                cTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAggg                                                                                                                                     >  1:2377886/1‑150 (MQ=255)
                cTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAggg                                                                                                                                     >  1:1934250/1‑150 (MQ=255)
                 tGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCGAGGgc                                                                                                                                    >  1:2820985/1‑150 (MQ=255)
                 tGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGgc                                                                                                                                    >  2:4366163/1‑150 (MQ=255)
                 tGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGgc                                                                                                                                    <  1:5295062/150‑1 (MQ=255)
                          aCGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCa                                                                                                                           >  2:1495160/1‑150 (MQ=255)
                          aCGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCa                                                                                                                           <  2:789191/150‑1 (MQ=255)
                          aCGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCa                                                                                                                           >  1:2761912/1‑150 (MQ=255)
                           cGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCaa                                                                                                                          <  1:2509364/150‑1 (MQ=255)
                           cGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCaa                                                                                                                          <  1:2328360/150‑1 (MQ=255)
                            gCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAg                                                                                                                         <  2:1517928/150‑1 (MQ=255)
                             ccGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGc                                                                                                                        >  1:544198/1‑150 (MQ=255)
                                  gcggcAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAgggg                                                                                                                   >  1:89335/1‑150 (MQ=255)
                                  gcggcAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAgggg                                                                                                                   >  1:2031785/1‑150 (MQ=255)
                                    ggcAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGgcg                                                                                                                 >  2:1652421/1‑150 (MQ=255)
                                     gcAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGgcgc                                                                                                                                  >  1:86666/1‑132 (MQ=255)
                                     gcAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGgcgc                                                                                                                <  1:3629868/150‑1 (MQ=255)
                                     gcAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGgcgc                                                                                                                <  1:3257033/150‑1 (MQ=255)
                                       aGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTg                                                                                                              >  2:3051392/1‑150 (MQ=255)
                                       aGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTg                                                                                                              >  1:2906042/1‑150 (MQ=255)
                                       aGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTg                                                                                                              <  1:467183/150‑1 (MQ=255)
                                       aGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGAGGCATTGCCAtcgtcg                                                                                                                                           <  2:419845/121‑1 (MQ=255)
                                         aCGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGt                                                                                                            >  2:4699135/1‑150 (MQ=255)
                                          cggcggGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGtt                                                                                                           <  2:2203669/150‑1 (MQ=255)
                                           ggcgggcgTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTc                                                                                                          >  2:183383/1‑150 (MQ=255)
                                             cgggcgTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGa                                                                                                        <  2:2341471/150‑1 (MQ=255)
                                               ggcgTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACt                                                                                                      <  2:259195/150‑1 (MQ=255)
                                                 cgTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTc                                                                                                    <  1:1633004/150‑1 (MQ=255)
                                                   tACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTa                                                                                                  <  2:2143599/150‑1 (MQ=255)
                                                             gCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTgg                                                                                                             >  1:4838826/1‑129 (MQ=255)
                                                              ccGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAAcc                                                                                       <  1:3473348/150‑1 (MQ=255)
                                                              ccGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAAcc                                                                                       <  1:2167758/150‑1 (MQ=255)
                                                                gggCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTg                                                                                     <  2:4787101/150‑1 (MQ=255)
                                                                  gcgcGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCa                                                                                   <  2:2232798/150‑1 (MQ=255)
                                                                   cgcgAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAg                                                                                  <  2:3645305/150‑1 (MQ=255)
                                                                    gcgAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCaa                                                                                                                                <  2:1775226/103‑1 (MQ=255)
                                                                    gcgAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGTCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGc                                                                                 <  2:4822625/150‑1 (MQ=255)
                                                                            aTGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTAt                                                                         <  1:2368267/150‑1 (MQ=255)
                                                                            aTGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTAt                                                                         >  1:261599/1‑150 (MQ=255)
                                                                            aTGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTAt                                                                         >  1:301167/1‑150 (MQ=255)
                                                                            aTGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTAt                                                                         >  1:1425295/1‑150 (MQ=255)
                                                                               ttCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGc                                                                      >  2:1886766/1‑150 (MQ=255)
                                                                                tCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCa                                                                     >  1:4204954/1‑150 (MQ=255)
                                                                                 cAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCaa                                                                    >  2:4274065/1‑150 (MQ=255)
                                                                                 cAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGACGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCaa                                                                    <  2:2954327/150‑1 (MQ=255)
                                                                                     gCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAAccc                                                                <  1:156226/150‑1 (MQ=255)
                                                                                       cTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTg                                                              >  1:2719143/1‑150 (MQ=255)
                                                                                        ttGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGc                                                             >  2:5157468/1‑150 (MQ=255)
                                                                                             cGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGcc                                                        >  1:1734527/1‑150 (MQ=255)
                                                                                              gCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGccc                                                       <  2:379890/150‑1 (MQ=255)
                                                                                                      cGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGa                                               >  2:3000035/1‑150 (MQ=255)
                                                                                                      cGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGa                                               <  1:2493209/150‑1 (MQ=255)
                                                                                                        aTGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAgg                                             <  1:5262213/150‑1 (MQ=255)
                                                                                                          gCACTGGCCGCTTTCTATGGCGAAATCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGgc                                           >  1:1859578/1‑150 (MQ=255)
                                                                                                           cACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGgcc                                          <  1:4869283/150‑1 (MQ=255)
                                                                                                            aCTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAg                                              >  2:3903450/1‑145 (MQ=255)
                                                                                                            aCTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAg                                              <  1:3903450/145‑1 (MQ=255)
                                                                                                            aCTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGgcca                                         <  2:3881041/150‑1 (MQ=255)
                                                                                                             cTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGgccag                                        <  1:330484/150‑1 (MQ=255)
                                                                                                              tGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGgccagc                                       <  2:2173881/150‑1 (MQ=255)
                                                                                                                   gCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCagccagcc                                  >  1:2408722/1‑150 (MQ=255)
                                                                                                                    cTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCg                                 <  2:1425295/150‑1 (MQ=255)
                                                                                                                     tttCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGa                                <  1:3710393/150‑1 (MQ=255)
                                                                                                                     tttCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGa                                <  1:3659078/150‑1 (MQ=255)
                                                                                                                      ttCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGAt                               <  2:4828630/150‑1 (MQ=255)
                                                                                                                       tCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGAtt                              <  2:4544450/150‑1 (MQ=255)
                                                                                                                         tATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTgc                            <  1:4438871/150‑1 (MQ=255)
                                                                                                                           tGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTgcgc                          <  2:429793/150‑1 (MQ=255)
                                                                                                                           tGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTgcgc                          <  2:893588/150‑1 (MQ=255)
                                                                                                                           tGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTgcgc                          <  2:4668293/150‑1 (MQ=255)
                                                                                                                             gCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTgcgcgc                        >  1:4230189/1‑150 (MQ=255)
                                                                                                                               gAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGcgccg                      >  2:4935278/1‑150 (MQ=255)
                                                                                                                                aaaGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGcgccgc                     <  2:5299069/150‑1 (MQ=255)
                                                                                                                                     ccttgcgacgaGGCATTGCCATCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTg                >  2:2891504/6‑150 (MQ=255)
                                                                                                                                       ttgcgacgaGGCATTGCCATCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGccgccg                   >  1:3833512/4‑145 (MQ=255)
                                                                                                                                       ttgcgacgaGGCATTGCCATCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCt              >  1:2673128/4‑150 (MQ=255)
                                                                                                                                         acgacgATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTgg            >  2:4599199/1‑150 (MQ=255)
                                                                                                                                             cgaTGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGc        >  1:4880759/1‑150 (MQ=255)
                                                                                                                                              gaTGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGcc       >  2:741771/1‑150 (MQ=255)
                                                                                                                                                tGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGa     >  1:1136839/1‑150 (MQ=255)
                                                                                                                                                 ggCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGAc    >  2:3734920/1‑150 (MQ=255)
                                                                                                                                                   cAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCg  >  2:5205106/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCACCCGTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCG  >  NC_002947/1289708‑1290004

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: