Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,542,364 C→T 7.0% P134P (CCG→CCA)  galE ← UDP‑glucose 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,542,3640CT 7.0% 289.3 / 8.4 100P134P (CCG→CCA) galEUDP‑glucose 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (48/45);  new base T (4/3);  total (52/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

GGTTTCCGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGCG  >  NC_002947/3542220‑3542502
                                                                                                                                                |                                                                                                                                          
ggTTTCCGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGCCAGCGCCCGCTGCATTGAAATATCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGtt                                                                                                                                       <  1:623449/150‑1 (MQ=255)
    tctggCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCGTATGGGTTGagc                                                                                                                                   <  1:1952752/147‑3 (MQ=255)
    gccGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCTGGTGGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATc                                                                                                                                   <  2:211332/149‑1 (MQ=255)
       ggCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGATTACCGAGTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGc                                                                                                                                >  2:1327446/1‑150 (MQ=255)
          tCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGcccc                                                                                                                             >  2:196123/1‑150 (MQ=255)
            gTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCtt                                                                                                                           >  1:1548082/1‑150 (MQ=255)
            gTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCtt                                                                                                                           >  1:453568/1‑150 (MQ=255)
              gACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGATCAACCCGTATGGGTTGagcaaat                                                                                                                               <  2:1266106/144‑7 (MQ=255)
              gACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGATCAACCCGTATGGGTTGagcaaat                                                                                                                               >  1:1266106/1‑138 (MQ=255)
               aCGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGgt                                                                                                                        <  1:2352689/150‑1 (MQ=255)
                   tCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGgtgtgt                                                                                                                    <  2:416724/150‑1 (MQ=255)
                    cGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGtgtgtg                                                                                                                   <  1:1327446/150‑1 (MQ=255)
                     gCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGgtgtgt                                                                                                                    >  1:521372/1‑148 (MQ=255)
                     gCCTAGCTGGCCTTCCGGGTCAGCGACCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGgtgtgt                                                                                                                    <  2:521372/148‑1 (MQ=255)
                      ccTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCt                                                                                                                 >  2:909981/1‑150 (MQ=255)
                          gCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCat                                                                                                             <  2:1854692/150‑1 (MQ=255)
                          gCTGGCCTTCCGGGTCAGCGCCCGATGCATTGCACAAACTCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGCAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCat                                                                                                             >  1:1380166/1‑150 (MQ=255)
                           cTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATc                                                                                                                                   >  2:1869605/1‑127 (MQ=255)
                           cTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATc                                                                                                                                   <  1:1869605/127‑1 (MQ=255)
                             ggCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGG‑CCCTTGGTGTGTGCTTCatcgat                                                                                                         <  1:1039946/150‑1 (MQ=255)
                               ccTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatcgatc                                                                                                        <  1:1704980/150‑1 (MQ=255)
                                cTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatcgatcg                                                                                                       <  2:1578238/150‑1 (MQ=255)
                                cTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatcgatcg                                                                                                       >  2:2121094/1‑150 (MQ=255)
                                cTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatcgatcg                                                                                                       >  1:1251982/1‑150 (MQ=255)
                                  tCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGc                                                                                                     <  2:1101478/150‑1 (MQ=255)
                                    cGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCCACGGCTTGGTGTGTGCTTCATCGATCGGCAc                                                                                                   <  1:1881503/150‑1 (MQ=255)
                                      ggTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACgg                                                                                                                                          <  2:1305657/109‑1 (MQ=255)
                                      ggTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACgg                                                                                                                                          >  1:1305657/1‑109 (MQ=255)
                                          aGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTg                                                                                             <  1:2400108/150‑1 (MQ=255)
                                           gcgcCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatcgatc                                                                                                        >  1:171806/1‑138 (MQ=255)
                                           gcgcCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatcgatc                                                                                                        <  2:171806/138‑1 (MQ=255)
                                                  cTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTAAAGAATCTGCTCGACCATCCATTTGCTCACCCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGc                                                                                     <  1:26738/150‑1 (MQ=255)
                                                   tGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGcc                                                                                    <  1:306743/150‑1 (MQ=255)
                                                      aTTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGG‑CCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGtat                                                                                >  2:2176142/1‑150 (MQ=255)
                                                        ttgaaTAGGGCAAGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGtata                                                                               <  2:1317610/147‑1 (MQ=255)
                                                           aaTAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATAcgg                                                                            >  1:342493/1‑150 (MQ=255)
                                                            aTAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCTAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCTATCGGCACGTACTGCGGGTTGCCGTATAcggc                                                                           <  1:564601/150‑1 (MQ=255)
                                                            aTAGCGCAGGCTTACCGACGTCAAACCATAGGGCCGGGCGAAGTCCTCAAGAATCTGCTCGAACATCCATTTGCTCAACCGTTGCTGGTTTGCCTGTCCCTTGTTGTGTGTTTCAGCGTTCGGCATGTACGGCGGGTTGCCGTATAcggc                                                                           <  2:272061/150‑1 (MQ=255)
                                                             gaGCGCAGTCGTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATAcggcg                                                                          >  2:2241735/2‑150 (MQ=255)
                                                              aGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACgg                                                                                                                                          >  2:1608228/1‑85 (MQ=255)
                                                              aGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACgg                                                                                                                                          <  1:1608228/85‑1 (MQ=255)
                                                                     gCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCCTGCAACGGCGCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCga                                                                  >  2:919699/1‑150 (MQ=255)
                                                                     gCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCga                                                                  >  1:259857/1‑150 (MQ=255)
                                                                        tACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCgaaga                                                               <  1:1168436/150‑1 (MQ=255)
                                                                           cGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCGTATGGGTTGa                                                                                                                                     <  2:2146017/77‑1 (MQ=255)
                                                                           cGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCGTATGGGTTGa                                                                                                                                     >  1:2146017/1‑77 (MQ=255)
                                                                           cGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCgaaga                                                               >  2:1579728/1‑147 (MQ=255)
                                                                           cGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCgaaga                                                               <  1:1579728/147‑1 (MQ=255)
                                                                           cGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGaaaa                                                            >  1:2340127/1‑150 (MQ=255)
                                                                           cGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGaaaa                                                            >  2:787098/1‑150 (MQ=255)
                                                                           cGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGaaaa                                                            >  1:460504/1‑150 (MQ=255)
                                                                             aCTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGACGAAAaca                                                          >  2:1144062/1‑150 (MQ=255)
                                                                             aCTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAaca                                                          >  2:540195/1‑150 (MQ=255)
                                                                                tCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCa                                                       <  1:196123/150‑1 (MQ=255)
                                                                                      cATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGcc                                                 >  1:1539596/1‑150 (MQ=255)
                                                                                       aTAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCCGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCt                                                <  2:1015427/150‑1 (MQ=255)
                                                                                           gCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAac                                                           <  1:97488/135‑1 (MQ=255)
                                                                                           gCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAac                                                           >  2:97488/1‑135 (MQ=255)
                                                                                              cGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCg                                         <  2:342493/150‑1 (MQ=255)
                                                                                                 tCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCa                                      <  2:1548082/150‑1 (MQ=255)
                                                                                                   gAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCTGGTTGCCGTATACGGCGGCACTCGAAGAAAACACTAGGTGCCTGATGCCGGCAtt                                    >  1:1208274/1‑150 (MQ=255)
                                                                                                    aaGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTg                                   >  1:1919704/1‑150 (MQ=255)
                                                                                                         cTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCAt                              >  2:2276685/1‑150 (MQ=255)
                                                                                                         cTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCAt                              >  2:1776289/1‑150 (MQ=255)
                                                                                                          tCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGCACGGACCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATg                             >  1:1879215/1‑150 (MQ=255)
                                                                                                          tCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGa                                  <  1:920075/145‑1 (MQ=255)
                                                                                                          tCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGa                                  >  2:920075/1‑145 (MQ=255)
                                                                                                          tCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATg                             >  1:2204752/1‑150 (MQ=255)
                                                                                                          tCAAGAATCTGCTCGACCATCCATTTACTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATg                             >  2:2395796/1‑150 (MQ=255)
                                                                                                          tCAAGAAGCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATg                             <  2:502839/150‑1 (MQ=255)
                                                                                                           cAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATgg                            <  1:2241735/150‑1 (MQ=255)
                                                                                                            aaGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGc                           >  1:2006572/1‑150 (MQ=255)
                                                                                                            aaGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGc                           >  1:395024/1‑150 (MQ=255)
                                                                                                                 tCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTgca                      >  2:526675/1‑150 (MQ=255)
                                                                                                                   tGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTgcagc                    >  1:1130543/1‑150 (MQ=255)
                                                                                                                     cTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGTGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGcagcag                  >  2:804313/1‑150 (MQ=255)
                                                                                                                     cTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGcagcag                  <  2:1130543/150‑1 (MQ=255)
                                                                                                                     cTCGACCATCCATTTGCTCAACCCATACGGCTTGCTCCGCGCTTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTTCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGcagcag                  <  2:1380166/150‑1 (MQ=255)
                                                                                                                      tCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGa                                  <  1:1604018/133‑1 (MQ=255)
                                                                                                                      tCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGa                                  >  2:1604018/1‑133 (MQ=255)
                                                                                                                      tCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAgg                 <  2:2204752/150‑1 (MQ=255)
                                                                                                                      tCGACCATCCATTTGCTCAACCCATACGCGCTGACAGGCGAACTGGCAGTTCATTGGCCCAGGGACACGGACTGTGTGTTGCCGTATAAGGCGACACTCATAGAAATCACCAGGTGCCTGATGCCTTCATTGACCATGGCCTGCAGCAgg                 >  1:1225170/1‑150 (MQ=255)
                                                                                                                        gACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAg                                                      <  2:1699111/111‑1 (MQ=255)
                                                                                                                        gACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAg                                                      >  1:1699111/1‑111 (MQ=255)
                                                                                                                            atccatTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAgg           >  1:1797146/1‑150 (MQ=255)
                                                                                                                             tccatTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAggg          <  1:1776289/150‑1 (MQ=255)
                                                                                                                              ccatTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGTGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGt         <  1:804313/150‑1 (MQ=255)
                                                                                                                              ccatTTACTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGt         <  1:2395796/150‑1 (MQ=255)
                                                                                                                               catTTGCTCAATCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTg        <  1:1077938/150‑1 (MQ=255)
                                                                                                                                atTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGTCCAATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTgg       <  1:2121094/150‑1 (MQ=255)
                                                                                                                                atTTGCTCAACCCATACGGGCTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTgg       >  1:1210794/1‑150 (MQ=255)
                                                                                                                                    gatCAACCCGTATGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGc   >  1:1569896/3‑150 (MQ=255)
                                                                                                                                    gatCAACCAGTATGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGc   >  2:1971646/3‑150 (MQ=255)
                                                                                                                                     cTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTAAACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGCg  >  2:1046182/1‑150 (MQ=255)
                                                                                                                                           ccATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTgca                      >  1:820241/1‑124 (MQ=255)
                                                                                                                                           ccATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTgca                      <  2:820241/124‑1 (MQ=255)
                                                                                                                                           ccATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTACCTGATGCCGGCATTGACCATGGCCTgcagc                    <  2:1675461/126‑1 (MQ=255)
                                                                                                                                           ccATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTACCTGATGCCGGCATTGACCATGGCCTgcagc                    >  1:1675461/1‑126 (MQ=255)
                                                                                                                                             aTACGGGTTGATCGGCCCCTTGGTGTGTGCTTCat                                                                                                             <  2:2201862/35‑1 (MQ=255)
                                                                                                                                             aTACGGGTTGATCGGCCCCTTGGTGTGTGCTTCat                                                                                                             >  1:2201862/1‑35 (MQ=255)
                                                                                                                                                |                                                                                                                                          
GGTTTCCGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGCG  >  NC_002947/3542220‑3542502

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