Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,292,4580TA46.8% 96.2 / 52.9 77intergenic (+93/+51)galU/hnsglucose‑1‑phosphate uridylyltransferase/global DNA‑binding transcriptional dual regulator H‑NS
Reads supporting (aligned to +/- strand):  ref base T (38/3);  new base A (0/36);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 8.44e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTT  >  NC_000913/1292309‑1292604
                                                                                                                                                      |                                                                                                                                                  
acCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAAt                                                                                                                                                    >  1:1010292/1‑150 (MQ=255)
  ccTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATaa                                                                                                                                                  >  2:781182/1‑150 (MQ=255)
     tGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATaaaaa                                                                                                                                               >  2:314465/1‑150 (MQ=255)
      ggCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAAt                                                                                                                                             <  1:720945/151‑1 (MQ=255)
        cACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATcc                                                                                                                                           >  2:772648/1‑151 (MQ=255)
        cACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATcc                                                                                                                                           >  2:1057401/1‑151 (MQ=255)
         aCGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATccc                                                                                                                                          <  2:1197266/151‑1 (MQ=255)
          cGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCccg                                                                                                                                         >  1:995441/1‑151 (MQ=255)
           ggAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCccgc                                                                                                                                        >  1:239018/1‑151 (MQ=255)
                        tGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGAtt                                                                                                                           >  1:1094535/1‑151 (MQ=255)
                         ggCTTGAAGAAGAGATGGGCATTAAGAAGAAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTAAAACCGA‑AAAAAAAAAAACCCCGCCGCTGGCGGGAttt                                                                                                                          <  2:127658/151‑1 (MQ=255)
                               aagaagAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGGTTTTaa                                                                                                                       >  2:975374/1‑148 (MQ=255)
                               aagaagAGAGGGGCTTAAAGAATTAACTTCCGTATCGTTTTTATCCACAAAAGGGCGTTGAGCAATCGACGCGTTTTTTTTATAGTTTATTTTTATAAAATTGT‑ATAAACCGGAAAAAAAAAAAACCCCGCCGCTGGCGGGATTTTaa                                                                                                                       <  1:370624/148‑1 (MQ=255)
                                agaagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCa                                                                                                                    >  2:618357/1‑150 (MQ=255)
                                agaagaGAGGGGCTAAAAAAAGAAACATCCGTATGGGTTTTACCCACAAACCGGCGTTGAGAACTCGACGCTTTTTTTTTATAGTTTATTTTTATTAAATTTTCTAAACCCGG‑AAAAAAAAAAATCCCGCCGCTGGCGGGAtttt                                                                                                                         <  1:1172425/145‑1 (MQ=255)
                                 gaagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCa                                                                                                                    >  2:1151351/1‑149 (MQ=255)
                                  aagaGATGGGAATAAAGAAGTAACACCCGTATGGGTTTTACCCACGAACCGGCTTGGAGAAACCGACGCTTTTTTTTTATAGTTTTTTTTTATTAAATTGTCTAAACCCGG‑AAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGt                                                                                                                 <  2:647228/151‑1 (MQ=255)
                                         gggAATAAAGAATAAACATCCTTATCGGTGTTATCCACAAAACGGCTTTGAGCAATCGACGCCTTTTTTTTATAGCTTTTTTTTATTAAATTTTCTTAACCCGA‑AAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCa                                                                                                              <  1:301950/147‑1 (MQ=255)
                                         ggacATTAAGAAGTAACACCCGTATCGTTGTTCCCCACAAACCGTGGTTGAGCAATCGCCGCCTTTTTTTTATCGCTTATTTTTATAAAATTTTCTAAAACCGG‑AAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGc                                                                                                               <  1:500939/143‑1 (MQ=255)
                                               tAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaaa                                                                                                     >  1:733443/1‑150 (MQ=255)
                                                aagaatTATCACCCTTCTCAGTATTACCCCCAAAACGCGTTTGACAAATCGCCGCTTTTTTTTTATAGCTTTTTCTTATAAATTTTTCTAACCCCGA‑ACAAAAAAAATACCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATc                                                                                                           <  1:1244950/143‑1 (MQ=255)
                                                aagaagTAACATCCGTATGGGTTTTATCCACGAAACGGCTTTGAGAAACCGACGCCTTTTTTTTATAGCTTATTCTTATTAAATTTTCTAAAACCGG‑AAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaaaa                                                                                                    <  2:1008996/150‑1 (MQ=255)
                                                       aaCATCCTTATGGTTGTTATCCAAAAAACGGCGTGGAGAAACCGCCGCTTTTTTTTTATAGTTTATTCTTATTAAATTGTCTAAACCCGA‑AAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTg                                                                                            <  1:825142/151‑1 (MQ=255)
                                                       aaCATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGattatt                                                                                             >  2:731229/1‑150 (MQ=255)
                                                        aCATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGc                                                                                           >  1:969230/1‑151 (MQ=255)
                                                         cATCCTTATGGTTTTTACCCACGAAAGGGCTTTGACCAACCGACGCCTTTTTTTTATATTTTTTTCTTATTAAATTGTGTTAAACCGG‑ACAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGc                                                                                           <  2:191082/150‑1 (MQ=255)
                                                            ccttttCGGTGTTACCCCAAAAACGGTGTGGAGAAATCGCCGCCTTTTTTTTATAGCTTATTCTTATAAAATTGTCTAAACCGGG‑AAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGa                                                                                       <  1:295624/146‑1 (MQ=255)
                                                            ccttATGGGTGTTACCCACGAACCGGTGTTAAGCAACCTACGCCGTTTTTTTATACCTTATTCTTATAAAATTTTCTCAACCCGA‑ACAAAAAAAATTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCtt                                                                                         <  2:1112402/146‑1 (MQ=255)
                                                            ccttATCGGTTTTTCCCACAAAACGGCGTTGAGAAATCGACGCCGTTTTTTTATAGTTTTTTCTTTTAAATTTGTCTAAACCCGC‑AAAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCtt                                                                                         <  1:280777/146‑1 (MQ=255)
                                                              gTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGAt                                                                                      >  1:571611/1‑150 (MQ=255)
                                                               tttcgttttttaccACGAAACGGCGTTAAGCAATCGACGCTTTTTTTTTATAGCTTTTTCTTATAAAATTTTCTTAACCGGA‑AAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGAt                                                                                      <  1:781138/137‑1 (MQ=255)
                                                               taggggTTTTATCCACGAAAGGGCTTTGACCAATCGACGCGTTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGA‑AAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTg                                                                                        <  1:315335/143‑1 (MQ=255)
                                                                  cGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAgg                                                                                  >  1:310568/1‑150 (MQ=255)
                                                                   gttgtTATCCACGAAAGGGCGTTGACCAATCCACGCTTTTTTTTTATAGTTTATTTTTAATAATTTTTCAAAACCCGA‑AAAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAgg                                                                                  <  1:1111458/147‑1 (MQ=255)
                                                                   gttgtTACCCACAAAAGGGCTTTGAGCAATCGCCGCGGTTTTTTTATAGTTTATTCTTATTAAATTGTCTTAACCCGA‑AAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGa                                                                                 <  2:12990/148‑1 (MQ=255)
                                                                   ggTTTTACCCAGGAAACGGGGTTGAGCAATCACCGCGTTTTTTTTATAGCTTATTTTTATTAATTTGTCTTAACCCGG‑AAAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGaa                                                                                <  2:678406/151‑1 (MQ=255)
                                                                     tgtTATCCCCGAAACGGCTTGGAGAAAGCTACGCGGTTTTTTTATAGCTTTTTCTTATTAATTTTTCTAAACCCGA‑ACAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGaaa                                                                               <  2:1239600/150‑1 (MQ=255)
                                                                          tCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgtcg                                                                         >  2:366224/1‑151 (MQ=255)
                                                                           ccACAAAAGGGGGTTGAAAAAGCGAGGCTTTTTTTTTATATCTTTTTTTTATTAAATTGACAAAAACCGA‑AAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                        <  1:815746/151‑1 (MQ=255)
                                                                           ccACAAAACGGCGTTGAGCAATCGACGCCTTTTTTTTATAGCTTATTTTTATTAATTTGTCAAAAACCGA‑AAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgtc                                                                          <  1:994705/149‑1 (MQ=255)
                                                                                aaagggggtggaaaaaTCGACGCCGTTTTTTTAGAGCTTATTTTTATTAATTTGTATAAAACGGA‑AAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAggg                                                                     <  1:574819/135‑1 (MQ=255)
                                                                                    ggCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTa                                                                >  1:879008/1‑150 (MQ=255)
                                                                                    ggCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTAc                                                               >  1:318824/1‑151 (MQ=255)
                                                                                    ggCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGGGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTa                                                                >  1:545472/1‑150 (MQ=255)
                                                                                    ggCGTGGAAAAAGCGCCGCTGTTTTTTAATAGTTTATTCTTTTAAAATTGTCTTAACCCGC‑AAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGAttt                                                                 <  2:706901/149‑1 (MQ=255)
                                                                                       gTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCtt                                                            >  1:424424/1‑151 (MQ=255)
                                                                                         tgagaactggcggCCGTTTTTTAACAGCTTATTCTTATTAATTTGTCTAAAACCGA‑AAAAAAAAAATTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTg                                                           <  1:300200/138‑1 (MQ=255)
                                                                                          gAGCAATCGACGCCTTTTTTTTATAGTTTATTTTTATAAATTTGTCTTAAACCGA‑AAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGc                                                          <  1:290350/150‑1 (MQ=255)
                                                                                           aGCACTCGTCGCCGTTTTTTAAAACATTATTCTTATAAATTTTTCTAAAAGCGG‑AAAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCtt                                                            <  2:562778/147‑1 (MQ=255)
                                                                                           aGCAATCGCCGCCTTTTTTTTATAGTTTTTTTTTATAAATTTGTCTAAAACCGA‑AAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGc                                                          <  1:1243498/149‑1 (MQ=255)
                                                                                             cAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCAt                                                      >  1:256470/1‑151 (MQ=255)
                                                                                              aaTCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATc                                                     >  1:605599/1‑151 (MQ=255)
                                                                                               accggcGCAGTTTTTTTATAGCTTATTCTTATAAAATTTTCTTAACCCGA‑ACAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATcc                                                    <  1:484148/146‑1 (MQ=255)
                                                                                                      ccGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtt                                             >  1:808822/1‑151 (MQ=255)
                                                                                                       cGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGTAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCCTTGAGGGGTTTACCTTGCTCATAAATTGCttt                                            >  1:1016584/1‑151 (MQ=255)
                                                                                                         tttttttATAGTTTATTTTTTTTAATTTGTTTAAACCCGG‑AAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtttt                                           <  1:688601/150‑1 (MQ=255)
                                                                                                          ttttttATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTTATCCATTGCtttttt                                         >  2:844131/1‑151 (MQ=255)
                                                                                                          ttttttATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtttttt                                         >  1:67402/1‑151 (MQ=255)
                                                                                                           tttttATATCTTATTCTTATTAATTTGTCTTAACCCGG‑ACAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtttttt                                         <  1:156101/150‑1 (MQ=255)
                                                                                                           tttttATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGTAATCTACAAAAGATTATTTCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtttttt                                         >  1:983591/1‑150 (MQ=255)
                                                                                                              ttATAGCTTATTCTTATTAAATTTTCTTAACCCGG‑ACAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGAt                                      <  2:451001/150‑1 (MQ=255)
                                                                                                                 taGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTa                                    >  2:1141761/1‑149 (MQ=255)
                                                                                                                 taGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTAc                                   >  2:520692/1‑150 (MQ=255)
                                                                                                                    cttattcttattAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCt                               >  1:471934/1‑151 (MQ=255)
                                                                                                                    cttattcttattAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAAATTATTGCTTGATCAGGAAATCCTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAg                                 >  1:439877/1‑149 (MQ=255)
                                                                                                                     ttatttttataAAATTGTCTAAAACCGA‑AAAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCt                               <  1:941432/150‑1 (MQ=255)
                                                                                                                           ttattaATTTGTCTAAAACCGA‑AAAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTa                         <  1:787025/150‑1 (MQ=255)
                                                                                                                           ttattaAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGACTTTAAGCAAGTGCAATCTACCAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCCTTGCTTTTTTGATTACAGCTGGAg                           >  1:352894/1‑148 (MQ=255)
                                                                                                                               taAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGGGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGcc                    >  1:275578/1‑151 (MQ=255)
                                                                                                                                aaaTTGTCTAAACCCGA‑AAAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGc                     <  1:314457/149‑1 (MQ=255)
                                                                                                                                  aTTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCCTCCATTGCTTTTTTGATTACAGCTGGGGTACGGCCTTg                 >  1:1005330/1‑151 (MQ=255)
                                                                                                                                   ttGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTgg                >  1:4435/1‑151 (MQ=255)
                                                                                                                                    tGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAAATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCAACCATTGCTTTTTTGATTACAGCCGGAGGACGGCCt                   >  1:158859/1‑147 (MQ=255)
                                                                                                                                     ttCTAAAACCGA‑AAAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTgg                <  1:859442/148‑1 (MQ=255)
                                                                                                                                      tCTTAACCCGG‑ACAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGcc              <  2:355045/150‑1 (MQ=255)
                                                                                                                                         tAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTc          >  2:457868/1‑151 (MQ=255)
                                                                                                                                          aacccGG‑ACAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTc          <  2:499591/147‑1 (MQ=255)
                                                                                                                                                    aaTAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGtttt  >  1:933133/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
ACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGG‑ACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTT  >  NC_000913/1292309‑1292604

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