Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,906,608 G→C 8.3% intergenic (+25/+35) ygcG → / ← eno TPM domain protein, putative phosphatase/enolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,906,6080GC 8.3% 286.4 / 15.8 96intergenic (+25/+35)ygcG/enoTPM domain protein, putative phosphatase/enolase
Reads supporting (aligned to +/- strand):  ref base G (38/50);  new base C (5/3);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-01

AGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGA  >  NC_000913/2906459‑2906756
                                                                                                                                                     |                                                                                                                                                     
aGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAACGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCg                                                                                                                                                       <  1:1275893/150‑1 (MQ=255)
       agCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTg                                                                                                                                                <  1:1498236/150‑1 (MQ=255)
            cTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑tt                                                                                                                                          <  1:630536/150‑1 (MQ=255)
              ggTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCA‑tttt                                                                                                                                        <  2:1256882/150‑1 (MQ=255)
              ggTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑tttt                                                                                                                                        <  1:785749/150‑1 (MQ=255)
              ggTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑tttt                                                                                                                                        >  2:472353/1‑150 (MQ=255)
               gTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑ttttt                                                                                                                                       <  2:1585960/150‑1 (MQ=255)
               gTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGCTGGCGGCGCATCCGGCCGCTGTTAACTCAACCATAAAAAATGACAGCCCGGAGGCGGGCA‑ttttt                                                                                                                                       >  2:1100594/1‑150 (MQ=255)
                tAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCA‑tttttt                                                                                                                                      <  2:1204387/150‑2 (MQ=255)
                tAGAAGCTATCATTCAGACAACAGCGGTTCCGCTGGAGGCGCAGATCGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCA‑tttttt                                                                                                                                      <  1:94548/150‑2 (MQ=255)
                     gTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCa                                                                                                                                 >  1:638099/1‑150 (MQ=255)
                         tCATTCAGACAACAGCGGTTCCGGTTGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCTGCCCGGAGGCTGGCA‑TTTTTAAATCAGATa                                                                                                                             <  2:341106/150‑1 (MQ=255)
                          cATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATaa                                                                                                                            >  1:1461535/1‑150 (MQ=255)
                             tCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAagt                                                                                                                         <  1:122743/150‑1 (MQ=255)
                              cAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAagtc                                                                                                                        >  1:17257/1‑150 (MQ=255)
                               aGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAagtca                                                                                                                       >  2:1103618/1‑150 (MQ=255)
                                  caacaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAagtcagtc                                                                                                                    >  1:1469172/1‑150 (MQ=255)
                                   aacaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCt                                                                                                                   >  1:1176337/1‑150 (MQ=255)
                                   aacaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCt                                                                                                                   >  1:482442/1‑150 (MQ=255)
                                    acaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCAGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCACTCtt                                                                                                                  >  1:1448848/1‑150 (MQ=255)
                                       gCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATg                                                                                                               <  2:503460/150‑1 (MQ=255)
                                           gtCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTg                                                                                                           <  2:701627/149‑1 (MQ=255)
                                            tCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTgg                                                                                                          <  1:1570927/150‑1 (MQ=255)
                                              cGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGcc                                                                                                        >  1:225567/1‑150 (MQ=255)
                                               gCTGGAGGCTGAGATAGCTGCGGATTTAGCGTCGGGGGCGGTGATTCTGGAGGTGGCGGCGCATCCGGCGGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAATTCAGATAAAGTCAGTCTTATGCCTGGCCt                                                                                                       <  1:387153/150‑1 (MQ=255)
                                                cTGGTGGCTCAGATAGCGGCGGTTTTAGCGGCGGGGGCTGTTCTTCTTGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATATAGTCAGTCTTATGCCTGGCCtt                                                                                                      <  1:1090520/150‑1 (MQ=255)
                                                cTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATTCCTGGCCtt                                                                                                      >  1:1456449/1‑150 (MQ=255)
                                                  ggAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTg                                                                                                    <  2:534652/150‑1 (MQ=255)
                                                  ggAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTg                                                                                                    <  2:541362/150‑1 (MQ=255)
                                                  ggAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTg                                                                                                    <  2:609415/150‑1 (MQ=255)
                                                   gAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGa                                                                                                   <  1:898565/150‑1 (MQ=255)
                                                     ggCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGTGGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAt                                                                                                  <  2:1529684/150‑1 (MQ=255)
                                                     ggCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtc                                                                                                 <  2:17257/150‑1 (MQ=255)
                                                     ggCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtc                                                                                                 <  1:1709955/150‑1 (MQ=255)
                                                     ggCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtc                                                                                                 <  1:1365007/150‑1 (MQ=255)
                                                     ggCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCAT‑AAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtct                                                                                                >  2:1204339/1‑150 (MQ=255)
                                                      gCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCACTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtct                                                                                                <  2:139404/150‑1 (MQ=255)
                                                       cTCAGATAGCGGCGGATTTAGGGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTAAGTCTTATGCCTGGCCTTTGAtctc                                                                                               <  2:182426/150‑1 (MQ=255)
                                                       cTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtctc                                                                                               <  2:236035/150‑1 (MQ=255)
                                                        tCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATctct                                                                                              <  2:301904/150‑1 (MQ=255)
                                                        tCAGAACGCGGCGGAGTACGCGGCGGGGGCGGCTCTGCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATctct                                                                                              <  2:832669/150‑1 (MQ=255)
                                                              aGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGa                                                                                        <  1:476063/150‑1 (MQ=255)
                                                               gcggcgGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGAc                                                                                       <  1:1123357/150‑1 (MQ=255)
                                                                 ggcggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGGAACTCAACCATAAAAAATGCCAGACCGGAGCCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCGGGCCTTTGATCTCTTTACGACCg                                                                                     >  2:1103261/1‑150 (MQ=255)
                                                                     gATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGt                                                                                 >  1:1160034/1‑150 (MQ=255)
                                                                     gATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGt                                                                                 >  2:1281613/1‑150 (MQ=255)
                                                                              gcggGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCtt                                                                        >  1:541582/1‑150 (MQ=255)
                                                                                 ggggCGGTTCTTCTTGAGGTGGAGGCTCATCCGGACGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCttttt                                                                     <  2:169882/150‑1 (MQ=255)
                                                                                  gggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTc                                                                    >  1:1398285/1‑150 (MQ=255)
                                                                                   ggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTACCTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCGCTTTACGACCGTTGTACGGTGCTTTTTCg                                                                   >  1:1136279/1‑150 (MQ=255)
                                                                                   ggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGGCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCg                                                                   <  1:203010/150‑1 (MQ=255)
                                                                                    gCGGTTCTTCTGGAGGTGGCGGCACATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGc                                                                  <  1:96487/150‑1 (MQ=255)
                                                                                     cGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGACA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGcc                                                                 >  2:1051140/1‑150 (MQ=255)
                                                                                      ggTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGccc                                                                <  2:1054534/150‑1 (MQ=255)
                                                                                       gTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCa                                                               >  1:319264/1‑150 (MQ=255)
                                                                                           ttctGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGc                                                           >  2:986537/1‑150 (MQ=255)
                                                                                               ggAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttct                                                       >  1:1284776/1‑150 (MQ=255)
                                                                                                  ggTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCg                                                    <  1:1025620/150‑1 (MQ=255)
                                                                                                    tGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGAt                                                  <  1:1062257/150‑1 (MQ=255)
                                                                                                      gcggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATAc                                                <  2:1461535/150‑1 (MQ=255)
                                                                                                      gcggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATAc                                                <  2:245577/150‑1 (MQ=255)
                                                                                                       gggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACg                                               <  2:429277/149‑1 (MQ=255)
                                                                                                       cggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACg                                               >  1:1002520/1‑150 (MQ=255)
                                                                                                         gcgcATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGaa                                             <  1:1620722/150‑1 (MQ=255)
                                                                                                            cATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGATTTTTCGCCCAGAGCTTCTTCGAAACGAATCa                                          >  2:1043140/1‑150 (MQ=255)
                                                                                                             aTCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAg                                         <  1:929373/150‑1 (MQ=255)
                                                                                                                  gCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGt                                    >  1:1183723/1‑150 (MQ=255)
                                                                                                                  gCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGt                                    <  1:188324/150‑1 (MQ=255)
                                                                                                                  gCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGt                                    <  2:611000/150‑1 (MQ=255)
                                                                                                                    cGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTg                                  <  2:482442/150‑1 (MQ=255)
                                                                                                                    cGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTg                                  <  2:961340/150‑1 (MQ=255)
                                                                                                                    cGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTg                                  <  2:1176337/150‑1 (MQ=255)
                                                                                                                     gCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGt                                 >  2:686953/1‑150 (MQ=255)
                                                                                                                       tGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTAt                               >  1:436307/1‑150 (MQ=255)
                                                                                                                       tGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTAt                               <  1:244083/150‑1 (MQ=255)
                                                                                                                        ggTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTAtt                              <  1:898555/150‑1 (MQ=255)
                                                                                                                         gTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTAttt                             <  2:319264/150‑1 (MQ=255)
                                                                                                                           aaCTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAg                           <  2:156491/150‑1 (MQ=255)
                                                                                                                             cTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGGCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCa                         <  1:1103618/150‑1 (MQ=255)
                                                                                                                              tCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCaa                        >  1:1173512/1‑150 (MQ=255)
                                                                                                                                 aCCATAAAAAATGCCAGCCCGGAGGCTGTCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAacac                     <  1:1423109/150‑1 (MQ=255)
                                                                                                                                  gatttaaaaaTGCCAGCCTCCGGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTTGTTGTATTTAGCAACACg                    >  2:272826/6‑150 (MQ=255)
                                                                                                                                    aTAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGt                  >  2:396221/1‑150 (MQ=255)
                                                                                                                                     ttaaaaaTGCCAGCCTCCGGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTc                 >  2:108347/3‑150 (MQ=255)
                                                                                                                                      aaaaaaTGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCa                >  1:462774/1‑150 (MQ=255)
                                                                                                                                       aaaaaTGCCAGCCTCCGGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTAAGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGCTCAg               >  2:848749/1‑150 (MQ=255)
                                                                                                                                        aaaaTGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAAGCAGCTGGTTGTTTTTAGCAACACGGTCAGa              >  1:742577/1‑150 (MQ=255)
                                                                                                                                          aaTGCCAGCCTCCGGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCACCACGGTCAGAAc            >  2:540971/1‑150 (MQ=255)
                                                                                                                                           aTGCCAGCCTCCGGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACg           >  1:1667541/1‑150 (MQ=255)
                                                                                                                                                aGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCGTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTGTTTCGCCCAGAGCTTAGTCGATACGAAGCAGGTGGTGGTTTTTAGCAACCGGGTCAGAACGGCTCa      >  2:285065/1‑150 (MQ=255)
                                                                                                                                                aGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGc         >  1:1560777/1‑147 (MQ=255)
                                                                                                                                                aGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGc         <  2:1560777/147‑1 (MQ=255)
                                                                                                                                                aGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCa      >  1:1490932/1‑150 (MQ=255)
                                                                                                                                                  cccGGAGGCTGGCA‑TTTTTATATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATa    >  1:431753/1‑150 (MQ=255)
                                                                                                                                                   ccGTAGGCTGGAA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAg   <  2:1448848/150‑1 (MQ=255)
                                                                                                                                                   ccGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAg   >  1:198067/1‑150 (MQ=255)
                                                                                                                                                    cGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGAGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACTGTCAGAACGGCTCATAGa  >  2:1293035/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCA‑TTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGA  >  NC_000913/2906459‑2906756

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