Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 311,480 (T)8→7 intergenic (‑144/‑632) ecpR ← / → ykgL putative transcriptional regulator for the ecp operon/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913311,4730T.97.5% 368.4 / 3.3 81intergenic (‑137/‑639)ecpR/ykgLputative transcriptional regulator for the ecp operon/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (28/51);  total (29/52)
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.94e-01

TTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGTTTCATAAGGTTATTGGCAATATAACGTTCAAT  >  NC_000913/311325‑311618
                                                                                                                                                    |                                                                                                                                                 
ttGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGtt                                                                                                                                                  <  1:1078983/150‑1 (MQ=255)
    cATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGttttttt                                                                                                                                             >  1:385443/1‑151 (MQ=255)
         cACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTcccg                                                                                                                                         <  1:87533/150‑2 (MQ=255)
         cACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTcccg                                                                                                                                         >  1:923173/1‑149 (MQ=255)
             actTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAg                                                                                                                                   <  2:665416/151‑1 (MQ=255)
                ttCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAgtg                                                                                                                                 <  1:888628/150‑1 (MQ=255)
                      aCCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAatat                                                                                                                          <  1:915015/151‑1 (MQ=255)
                          gTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGa                                                                                                                      <  1:1005595/151‑1 (MQ=255)
                            aaTTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGaa                                                                                                                     <  1:144994/150‑1 (MQ=255)
                            aaTTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAAc                                                                                                                    >  2:96333/1‑151 (MQ=255)
                                    cAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAgg                                                                                                            <  1:727792/151‑1 (MQ=255)
                                    cAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAgg                                                                                                            >  1:1188838/1‑151 (MQ=255)
                                        tGCATTCAGTTTCCTTAGTAGTCCAGAAGTAAATAACAGTCGAAGTACATATCAATATAGTTTTAACAACTTTAATTTTAATTATTTGCAATTCGGCTTAAATTTATGTTTTTTTCCCGATCAGTGTCAATATGTGAACGGGAAAAGAAt                                                                                                         <  2:435883/150‑1 (MQ=255)
                                         gCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTc                                                                                                       >  2:1087223/1‑151 (MQ=255)
                                          cATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTc                                                                                                       >  1:679092/1‑150 (MQ=255)
                                               aGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGt                                                                                                 <  1:180602/151‑1 (MQ=255)
                                                        aGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATTGTTAAGTCAat                                                                                        <  2:863041/151‑1 (MQ=255)
                                                         gTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAata                                                                                       <  1:1163818/151‑1 (MQ=255)
                                                         gTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAata                                                                                       <  2:632368/151‑1 (MQ=255)
                                                           aaaCCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAatat                                                                                      <  2:1009608/150‑1 (MQ=255)
                                                             aCCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATAttt                                                                                    >  1:104867/1‑150 (MQ=255)
                                                              ccGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAg                                                                                  <  2:381759/151‑1 (MQ=255)
                                                                 gAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGaa                                                                                >  2:926457/1‑150 (MQ=255)
                                                                  aaGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGaaaa                                                                              <  2:923190/151‑1 (MQ=255)
                                                                   aGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGaaaaa                                                                             <  2:198972/151‑1 (MQ=255)
                                                                     tAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACt                                                                           >  2:154435/1‑151 (MQ=255)
                                                                      aaaTAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTg                                                                          >  1:583420/1‑151 (MQ=255)
                                                                       aaTAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGa                                                                         <  1:1085002/151‑1 (MQ=255)
                                                                       aaTAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGa                                                                         <  2:861691/151‑1 (MQ=255)
                                                                       aaTAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGa                                                                         <  1:1039275/151‑1 (MQ=255)
                                                                        aTAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGaa                                                                        <  1:564090/151‑1 (MQ=255)
                                                                        aTAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGa                                                                         >  1:1075444/1‑150 (MQ=255)
                                                                            gATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAata                                                                    <  2:1093371/151‑1 (MQ=255)
                                                                               aCGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATAtt                                                                  <  2:162603/150‑1 (MQ=255)
                                                                                cGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGc                                                                >  1:538684/1‑151 (MQ=255)
                                                                                cGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGc                                                                <  2:287452/151‑1 (MQ=255)
                                                                                 gtagtaCATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGGAAATATTGCt                                                               >  2:124575/1‑151 (MQ=255)
                                                                                   agtaCATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCttt                                                             <  2:858737/151‑1 (MQ=255)
                                                                                   agtaCATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCttt                                                             <  2:385446/151‑1 (MQ=255)
                                                                                    gtaCATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCtttt                                                            <  1:986351/151‑1 (MQ=255)
                                                                                        atatCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTg                                                         >  2:669829/1‑150 (MQ=255)
                                                                                          atCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGAt                                                       <  2:174281/150‑1 (MQ=255)
                                                                                           tCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATa                                                      <  1:1055885/150‑1 (MQ=255)
                                                                                            cAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATaa                                                     <  1:365243/150‑1 (MQ=255)
                                                                                              aTAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATaagaa                                                  <  1:943468/151‑1 (MQ=255)
                                                                                              aTAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATaaga                                                   >  1:1048442/1‑150 (MQ=255)
                                                                                               tAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATaaga                                                   >  1:743847/1‑149 (MQ=255)
                                                                                                aaaGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATaagaag                                                 <  2:1078716/150‑1 (MQ=255)
                                                                                                 aaGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAgaagaa                                               <  1:213576/151‑1 (MQ=255)
                                                                                                   ggTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGaaaa                                             <  2:1163012/151‑1 (MQ=255)
                                                                                                   ggTTAAACGAGATTAAGGTTAATTATTTTCGATTGGGGTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGAGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGaaa                                              >  1:939934/1‑150 (MQ=255)
                                                                                                      tAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATTAGTGTCATTATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGaaaa                                             >  1:787660/1‑148 (MQ=255)
                                                                                                      tAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGaaaa                                             >  2:628385/1‑148 (MQ=255)
                                                                                                      tAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAg                                           <  2:97060/150‑1 (MQ=255)
                                                                                                       aaaCGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGAc                                         <  1:208167/151‑1 (MQ=255)
                                                                                                        aaCGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAATAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGa                                          >  2:317654/1‑149 (MQ=255)
                                                                                                        aaCGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGAc                                         >  2:1170952/1‑150 (MQ=255)
                                                                                                        aaCGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGAc                                         >  1:162583/1‑150 (MQ=255)
                                                                                                        aaCGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGAc                                         >  2:235079/1‑150 (MQ=255)
                                                                                                         aCGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACAt                                       <  2:161721/151‑1 (MQ=255)
                                                                                                         aCGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACAt                                       <  1:408145/151‑1 (MQ=255)
                                                                                                          cGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACAtt                                      >  1:203543/1‑151 (MQ=255)
                                                                                                           gagaTTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACAttt                                     <  1:928178/151‑1 (MQ=255)
                                                                                                           gagaTTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACAttt                                     <  2:24450/151‑1 (MQ=255)
                                                                                                           gagaTTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACAttt                                     <  2:603485/151‑1 (MQ=255)
                                                                                                            agaTTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTg                                    <  1:687321/151‑1 (MQ=255)
                                                                                                             gaTTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTgg                                   <  1:454590/151‑1 (MQ=255)
                                                                                                             gaTTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTgg                                   <  2:454903/151‑1 (MQ=255)
                                                                                                              aTTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGt                                  <  1:774143/151‑1 (MQ=255)
                                                                                                                tAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGtt                                 <  1:1010800/150‑1 (MQ=255)
                                                                                                                 aaGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGTTTc                               <  1:1014273/151‑1 (MQ=255)
                                                                                                                          tattTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGTTTCATAAGGtt                       >  1:750290/1‑150 (MQ=255)
                                                                                                                           attTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGTTTCATAAGGTTa                      <  1:682731/150‑1 (MQ=255)
                                                                                                                           attTTCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGTTTCATAAGGTTa                      <  1:321855/150‑1 (MQ=255)
                                                                                                                            ttttCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGTTTCATAAGGTTAtt                    >  2:689886/1‑151 (MQ=255)
                                                                                                                             tttCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGTTTCATAAGGTTAtt                    <  1:339414/150‑1 (MQ=255)
                                                                                                                              ttCGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGTTTCATAAGGTTATTg                   <  2:274312/150‑1 (MQ=255)
                                                                                                                                cGATTGGGCTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGTTTCATAAGGTTATTGGc                 >  1:853372/1‑150 (MQ=255)
                                                                                                                                        cTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGTTTCATAAGGTTATTGGCAATATAACg        >  1:423194/1‑151 (MQ=255)
                                                                                                                                        cTTAAATTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGTTTCATAAGGTTATTGGCAATATAACg        >  2:1008909/1‑151 (MQ=255)
                                                                                                                                           aaaTTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGTTTCATAAGGTTATTGGCAATATAACGtt      <  1:679299/150‑1 (MQ=255)
                                                                                                                                           aaaTTTATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGTTTCATAAGGTTATTGGCAATATAACGTTc     <  1:378894/151‑1 (MQ=255)
                                                                                                                                              tttATGTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGTTTCATAAGGTTATTGGCAATATAACGTTCAAt  >  2:668873/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAGTGTCAATATGTGAACGGAAAAGGAATTCAATGGTTAAGTCAATATTTAGAAAAACTGAAAATATTGCTTTTTTGATAAGAAGAAAAAGACATTTGGTTTCATAAGGTTATTGGCAATATAACGTTCAAT  >  NC_000913/311325‑311618

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: