Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,998,723 A→T 10.7% L267H (CTT→CAT)  rffG ← dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,998,7230AT10.7% 214.2 / 16.0 75L267H (CTT→CAT) rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base T (4/4);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

GCTCGATCTTGCGTGCATCGATGGCATAACGGCGGTCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATCTTCGACATACAGCCAGTCGCGGATTTGCTGG  >  NC_002947/1998574‑1998869
                                                                                                                                                     |                                                                                                                                                   
gCTCGATCTTGCGTGCATCGATGGCATAACGGCGGTCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGaa                                                                                                                                                     <  2:1255166/150‑1 (MQ=255)
gCTCGATCTTGCGTGCATCGATGGCATAACGGCGGTCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGaa                                                                                                                                                     >  2:306306/1‑150 (MQ=255)
  tCGATCTTGCGTGCATCGATGGCATAACGGCGGTCATGGCCTGGGCGGTCGGTTACGTACGCGATGAGTTCTGCATACGGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCCTATGGTAGGGAAgg                                                                                                                                                   <  2:819870/150‑1 (MQ=255)
     aTCTTGCGTGCATCGATGGCATAACGGCGGTCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCTTCTCGACGAGATGGcac                                                                                                                                                <  1:2102628/150‑3 (MQ=255)
      tCTTGCGTGCATCGATGGCATAACGGCGGTCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGc                                                                                                                                               >  1:1570479/1‑150 (MQ=255)
            gTGCATCGATGGCATAGCGGCGGTCATGGCCTGGGCGGTCGGTTACGTAGGCTATGAGTTCTGCATACTGTTCGACAGGCTCGCCCGTCTTCTGTTTTATTACCTGGCGCTATACGGCGGGTGCCATCTCGTGGAGCAGGCCGCagagcg                                                                                                                                         >  1:1991120/1‑148 (MQ=255)
            gTGCATCGATGGCATAACTGCGGTTTTGGCGTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGNgtg                                                                                                                                         <  2:1496264/150‑1 (MQ=255)
                      ggCATAACGGCGGTCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCAcaat                                                                                                                              <  2:1488019/150‑1 (MQ=255)
                           aaCGGCGGTCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGc                                                                                                                                               >  1:1113805/1‑129 (MQ=255)
                           aaCGGCGGTCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGc                                                                                                                                               <  2:1113805/129‑1 (MQ=255)
                            aCGGCGGTCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCAcaatgtcaat                                                                                                                        <  2:1781016/150‑1 (MQ=255)
                                cggTCATGGCCTGGGCGGTCGGTTACGTAGGTAATGATACGGGGATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTg                                                                                                                    <  2:2018408/150‑1 (MQ=255)
                                cggTCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTg                                                                                                                    <  1:1290405/150‑1 (MQ=255)
                                cggTCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTg                                                                                                                    <  2:1570479/150‑1 (MQ=255)
                                cggGCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTTGCGCGATCACTCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTg                                                                                                                    <  1:428055/150‑1 (MQ=255)
                                 ggTCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTgg                                                                                                                   >  2:1849361/1‑150 (MQ=255)
                                 ggTCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTgg                                                                                                                   >  2:1850409/1‑150 (MQ=255)
                                  gTCCTGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGAAAACTATACGTAAAACATTTGATGCTTCTTTTTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGc                                                                                                                  <  1:877501/150‑1 (MQ=255)
                                  gTCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGc                                                                                                                  <  2:1193387/150‑1 (MQ=255)
                                   tCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGCTTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCt                                                                                                                 >  2:192192/1‑150 (MQ=255)
                                         ccTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCAcaatgt                                                                                                                            >  2:440704/1‑133 (MQ=255)
                                         ccTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCAcaatgt                                                                                                                            <  1:440704/133‑1 (MQ=255)
                                         ccTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCa                                                                                                           <  1:1850409/150‑1 (MQ=255)
                                         ccTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCa                                                                                                           <  1:1849361/150‑1 (MQ=255)
                                           tGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCAtt                                                                                                         <  1:1782792/150‑1 (MQ=255)
                                            gggCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTGGGCTTTTTCATTc                                                                                                        >  1:949973/1‑150 (MQ=255)
                                              gCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCa                                                                                                      <  1:306306/150‑1 (MQ=255)
                                                                 aTGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTg                                                                                                                                                <  2:1644294/90‑1 (MQ=255)
                                                                 aTGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTg                                                                                                                                                >  1:1644294/1‑90 (MQ=255)
                                                                 aTGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGtc                                                                                   <  2:266418/150‑1 (MQ=255)
                                                                 aTGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGtc                                                                                   >  2:311751/1‑150 (MQ=255)
                                                                      ttCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGcc                                                                              >  1:2391831/1‑150 (MQ=255)
                                                                      ttCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGcc                                                                              >  1:199592/1‑150 (MQ=255)
                                                                      ttCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGcc                                                                              >  2:1471467/1‑150 (MQ=255)
                                                                          gCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAAc                                                                          >  1:812992/1‑150 (MQ=255)
                                                                                    cGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCtt                                                                >  2:2208743/1‑150 (MQ=255)
                                                                                     gACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTc                                                               <  1:1471467/150‑1 (MQ=255)
                                                                                     gACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCGGCAGAGTGTACGCACAATGTCCGTGGTGGCTTTTTCATTCCAGCCGCCAAGGTTGTACGTCTCGCCGAACGCACCGGCTTc                                                               >  1:25837/1‑150 (MQ=255)
                                                                                        aGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCTTCTCGACGAGATGGca                                                                                                                                                 <  2:1461715/66‑2 (MQ=38)
                                                                                        aGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCTTCTCGACGAGATGGca                                                                                                                                                 <  2:1461746/66‑2 (MQ=38)
                                                                                        aGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCTTCTCGACGAGATGGca                                                                                                                                                 >  1:1461746/1‑65 (MQ=38)
                                                                                        aGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCTTCTCGACGAGATGGca                                                                                                                                                 >  1:1461715/1‑65 (MQ=38)
                                                                                          gCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCAAAGCTGTTATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTa                                                          >  1:1332984/1‑150 (MQ=255)
                                                                                            tCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTAcg                                                        >  1:1869308/1‑150 (MQ=255)
                                                                                               ccGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGa                                                     <  1:465105/150‑1 (MQ=255)
                                                                                                     ttctGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATg                                               <  2:1291875/150‑1 (MQ=255)
                                                                                                      tctGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGc                                              >  1:179725/1‑150 (MQ=255)
                                                                                                        tgattgatTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGccc                                            >  1:977763/1‑150 (MQ=255)
                                                                                                        tgattgatTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGccc                                            <  1:1568611/150‑1 (MQ=255)
                                                                                                        tgattgatTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGccc                                            <  1:1091735/150‑1 (MQ=255)
                                                                                                          attgattACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTc                                                               <  1:1855178/129‑1 (MQ=255)
                                                                                                          attgattACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTc                                                               >  2:1855178/1‑129 (MQ=255)
                                                                                                          attgattACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGa                                          <  2:199592/150‑1 (MQ=255)
                                                                                                           ttgattACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAAc                                                                          <  1:943268/117‑1 (MQ=255)
                                                                                                           ttgattACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAAc                                                                          >  2:943268/1‑117 (MQ=255)
                                                                                                               ttACCTGGCGCGATACCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGt                                     <  2:1991120/150‑1 (MQ=255)
                                                                                                                   cTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATc                                 >  1:1399179/1‑150 (MQ=255)
                                                                                                                      gcgcGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATCTTc                              <  2:179725/150‑1 (MQ=255)
                                                                                                                          gATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATCTTCGACa                          <  2:812992/150‑1 (MQ=255)
                                                                                                                                     agccttctcGACGAGATGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATCTTCGACATACAGCCAGTc               >  1:2214324/6‑150 (MQ=255)
                                                                                                                                      gCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATCTTCGACATACAGCCAGTcg              >  1:461689/1‑150 (MQ=255)
                                                                                                                                       ccttctcGACGAGATGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCa                                       >  1:1422736/4‑124 (MQ=255)
                                                                                                                                       ccttctcGACGAGATGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCa                                       <  2:1422736/121‑1 (MQ=255)
                                                                                                                                       ccATCTCGTCTAAAATGCTGCAGATTGTACTCACAATTTCAATGCTAACTTTTTCATTCCAGCCGCCAATATTATACGTCTCCACGAACGCAACTATTTCCATTCCGGGACGATTGCCGGTAATGTCATCTTCTACATTCAGCCATTcgc             <  2:1110458/150‑1 (MQ=255)
                                                                                                                                        cATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGCACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATCTTCGACATACAGCCAGTcgcg            >  2:70465/1‑150 (MQ=255)
                                                                                                                                         aTCTCGTCGATAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATCTTCGACATACAGCCAGTCGCgg           >  2:1912503/1‑150 (MQ=255)
                                                                                                                                          tctcGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATCTTCGACATACAGCCAGTCGCGGa          >  1:545333/1‑150 (MQ=255)
                                                                                                                                          tctcGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCCACGGATGCCCAAACCGTGATCGTCGACATACAGCCAGCcgcgtt          >  2:514719/1‑148 (MQ=255)
                                                                                                                                          aatcgtcgAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCa                                                             >  2:1827524/3‑99 (MQ=255)
                                                                                                                                          aatcgtcgAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCa                                                             <  1:1827524/97‑1 (MQ=255)
                                                                                                                                              gtcgAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATCTTCGACATACAGCCAGTCGCGGATTTg      >  2:1563509/1‑150 (MQ=255)
                                                                                                                                               tcgAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATCTTCGACATACAGCCAGTCGCGGATTTGc     <  2:949973/150‑1 (MQ=255)
                                                                                                                                               tcgAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATCTTCGACATACAGCCAGTCGCGGATTTGc     >  1:18655/1‑150 (MQ=255)
                                                                                                                                                cgAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATCTTCGACATACAGCCAGTCGCGGATTTGCt    <  2:18655/150‑1 (MQ=255)
                                                                                                                                                 gagaAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATCTTCGACATACAGc                    >  1:1446561/1‑133 (MQ=255)
                                                                                                                                                 gagaAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATCTTCGACATACAGc                    <  2:1446561/133‑1 (MQ=255)
                                                                                                                                                  agaAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATCTTCGACATACAGCCAGTCGCGGATTTGCTgg  <  1:1563509/150‑1 (MQ=255)
                                                                                                                                                    aaGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGAt                                                <  2:2136972/102‑1 (MQ=255)
                                                                                                                                                    aaGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGAt                                                >  1:2136972/1‑102 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCTCGATCTTGCGTGCATCGATGGCATAACGGCGGTCATGGCCTGGGCGGTCGGTTACGTAGGCGATGAGTTCTGCATACTGTTCGACAGGCTCGCCGGTCTTCTGATTGATTACCTGGCGCGATGCCGCAGGTGCCATCTCGTCGAGAAGGCTGCAGAGTGTACGCACAATGTCAATGTTGGCTTTTTCATTCCAGCCGCCAATATTGTACGTCTCGCCGAACGCACCGGCTTCCAGTACGCGACGGATGCCCGAGCAGTGATCTTCGACATACAGCCAGTCGCGGATTTGCTGG  >  NC_002947/1998574‑1998869

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: