Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,260,8880CA23.2% 90.7 / 10.2 57intergenic (‑85/+40)dauA/prsaerobic C4‑dicarboxylate transporter DauA/ribose‑phosphate diphosphokinase
Reads supporting (aligned to +/- strand):  ref base C (35/8);  new base A (0/13);  total (35/22)
Fisher's exact test for biased strand distribution p-value = 1.08e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCACTACCAATAGCCAACGCCATCGCCAGCGGGATAGCAATAATCCCGACGGTTATCCCGGCAATCAGGTCACGGGTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGTATCGCAGACAACGACTTCATCAATTACAGAGTTACGCAGGTTGTTCGCCGCGTTGCCAGAGAAGATCGGGTGAGTCGCGTA  >  NC_000913/1260640‑1261137
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                         
ccacTACCAATAGCCAACGCCATCGCCAGCGGGAGAGCAATAATCCCGACGGTTATCCCGGCAATCAGGTCACGGGTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACaa                                                                                                                                                                                                                                                         <  2:369817/251‑1 (MQ=255)
  acTACCAATAGCCAACGCCATCGCCAGCGGGATAGCAATAATCCCGACGGTTATCCCGGCAATCAGGTCACGGGTAAACCGTGCGGCAGTATCTTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACaaa                                                                                                                                                                                                                                                        >  1:682702/1‑250 (MQ=255)
     aCCAATAGCCAACGCCATCGCCAGCGGGATAGCAATAATCCCGACGGTTATCCCGGCAATCAGGTCACGGGTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACaaaaa                                                                                                                                                                                                                                                      <  1:455369/249‑1 (MQ=255)
        aaTAGCCAACGCCATCGCCAGCGGGATAGCAATAATCCCGACGGTTATCCCGGCAATCAGGTCACGGGTAAACCGTGCGGCAGTATATTTTTCTTGCCAGCAAGCGTCGAGCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGGCAACTAAATGGCAGGTGAAGGAGGCATAGGGCATACAAATGGATAGTACATGCTAAAAAAccc                                                                                                                                                                                                                                                 <  2:36318/251‑1 (MQ=255)
         aTAGCCAACGCCATCGCCAGCGGGATAGCAATAATCCCGACGGTTATCCCGGCAATCAGGTCACGGGTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACccg                                                                                                                                                                                                                                                >  1:612736/1‑251 (MQ=255)
                aCGCCATCGCCAGCGGGATAGCAATAATCCAGACGGTTATCCCGGCCATCCGGTCACGGGTAAACCGTGCGGCCGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGACAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGACATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGc                                                                                                                                                                                                                                         >  2:631118/1‑251 (MQ=255)
                     aTCGCCAGCGGGATAGCAATAATCCCGACGGTTCTCCCGGCACTCAGGTCACGGGTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCACGCGCCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTCCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGATGGCGGCATAGGTCATACAAATGGAAATTACAGACAAAACACACCGCCGCAGCGGGTc                                                                                                                                                                                                                                    >  2:386626/1‑251 (MQ=255)
                                gATAGCAATAATCCCGACGGTTATCCCGGCAATCAGGTCACGGGTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCggg                                                                                                                                                                                                                         >  1:461485/1‑251 (MQ=255)
                                  aaGCAATAATCCCGACGTTTATCCCGGCAATCAGGTCACGGGTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGAAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATTTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGAATAGGTCATAAAAATGGATATTACGGAAAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGtt                                                                                                                                                                                                                       <  1:467041/250‑1 (MQ=255)
                                      aataatCCCGACGGTTATCCCGGCAATCAGGTCACGGGTAAACCGTGCGGCAGAATATTTTTCTTTCCAGCAAGAGCCGATCAGAGCGCGGAAAGCCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAATTAAATGGCAGGTGAAGAGGGCAAAGGTCATAAAAAGGGATATTACGGACAAAAAAACCGGCGGCACCGGGTCTTTGAGCCGGGTTCgct                                                                                                                                                                                                                   >  1:240461/1‑249 (MQ=255)
                                       ataatCCCGACGGATATCCCGGCAATCAGGTCACGGGTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGAtt                                                                                                                                                                                                                  >  1:577043/1‑251 (MQ=255)
                                                                  agtcaaCGGGTAAACCTGGCGGCAGTATTTTTGTCTTTCAAGCGAGCGTGTGCCAGGGCAGAGAAAGGCTTCAAATGGGAGAAAAATTTTTTGTTCACAATAATGTTTCTGCCGGGGGCGCATCGTGTTTCAACAAAATGGCGGGGGAGGGAGGCAGAGGTCATAAAATTGGATTTTACAAAAAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAtt                                                                                                                                                                                       <  2:521789/246‑1 (MQ=255)
                                                                  aGGTCACGGGTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAtt                                                                                                                                                                                       >  2:655569/1‑251 (MQ=255)
                                                                      caggggTAATCCGTGCGGCATAATATTTGTCTTCCCGGCAAGCGTCGTGCAGGGCGCGGAAGGGCATAACATGTGAGAAAAATTTTTTGTAAACAAAAATTTTTACTCCGTGAGCGCATCGTCTGTGAGCTAAGTGGGAGGGAAATGAGGCATGGGCAGTAAAAGTGGATATTACAGATAAAAATTGCCGCCGCAGCGGGTCTTTGAGCCGGTTTCGATTAGTGTTCGAACATGGCAGAGATCGAttctt                                                                                                                                                                                    <  1:691316/247‑1 (MQ=255)
                                                                              aaCCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAgtgt                                                                                                                                                                                                             >  1:490408/1‑217 (MQ=255)
                                                                                 cGTGCGGCCGTCTATTTTTCTTTACAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCCTAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCAGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATac                                                                                                                                                                        >  2:432825/1‑251 (MQ=255)
                                                                                         aGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTc                                                                                                                                                                                                                                    >  1:155691/1‑183 (MQ=255)
                                                                                          gtaatgttttcgccccgtgaggCGTCGATCAGCGCGCTAAAAGGCATCACAGGTAGGAAAAATTTTTTGTACAAAATAATGTTTAATCCGTGGGTGCAAGCTATGTAAAAAAAAGTGCAGGAGAAGGCGACATAGGAAATACAAATGGAATGTACAGAAAAAAAAGCCCGCCGCAGCGGGTCTGTGAGCCGGGTTCGAGTAGTGTTCGAACATGGCAGAGATCGATTCTTCGGTGCTTATACGGCGAATCg                                                                                                                                                               <  2:335952/230‑1 (MQ=255)
                                                                                              ttttttCTTCCCCGCAAGGGTCGAACAGAGCGCGGAAGGGCATCACATGGGAGAAAAATTTTTTGTCCAAAAAAATTTTTCATCCGGGGGCGCAACAGCTGTCAAATAAATGGCAGGTAAAGGAGGCATGGGTAATAAAAAGGTATATTAAAGAAAAAAAACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTGCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCtt                                                                                                                                                            <  2:597024/249‑1 (MQ=255)
                                                                                                   tcttccCAGCAACCGGCGATCAGAGCGCGAAAAGGCAGCTCTTGTGAGGAAATTTTTTTTTTCATAATATTGTTTCCCCCGTGGGGGATTCACCTTCCAACTAAATGGCGGTTGAAGGAGGAATGGTGCATGCAAGGGGAGATGCCAGCCAAAAAAACCCGCCGCAGCGGGGCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCTGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGcc                                                                                                                                                       <  1:633930/250‑1 (MQ=255)
                                                                                                      ttcCAGCAAGCGTCGATCAGAGCGCGGAAAGGAATCACATGTGAGAAAAATATTTTGTTCACAATAATGTTTCACCCGTGAGCGCATCATCTGTAAACAAAAGGGCAGGGAAAGGAGGCTTAGGTCATAAAAAGGGATATTACAGACAAAAAACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAg                                                                                                                                                     <  1:15691/249‑1 (MQ=255)
                                                                                                      taccAGCAAGCGTCGATCAGAGCGCGAAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGAAAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGc                                                                                                                                                    <  1:591531/248‑1 (MQ=255)
                                                                                                      ccccAGCAAGCGGCGATCAGAGCGCGAAAAGGCACAAAATGGGAGAAAAATTTTTTGTTCACAATAATGTTTCATCGGGGAGCGCATCATCTGTCAACAAAAGGGAAGGTGAAGGGGGCATAGGTAAAAAAAATGAATATACAAGAAAAAAAACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCa                                                                                                                                                   <  2:329200/249‑1 (MQ=255)
                                                                                                             aaGCGTCGATCAGAGCGCGGAAAGGCAGCACATGTGAGAAAAATATTTTGTTCAAAATAATGTTTCATCGGTGAGCGCATCATCTGTAAAATAAATGGCAGGTGAAGGGGGCATAGGTCATACAAATGGATATTACAGAAAAAAAACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGa                                                                                                                                            <  2:127471/251‑1 (MQ=255)
                                                                                                              aGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTGATTAGTGTTCGAACATGGCag                                                                                                                                                                                                >  1:547800/1‑197 (MQ=255)
                                                                                                              aGCGTCGATCAGAGCGCGGAAAGGCATCACAGGTGAGAAAAATATTTTGTTCAAAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACAAAAGGGCAGGTGAAGGAGGCATAGGTCTTAAAAATGAATATTACAGAAAAAAAACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTGATTAGTGTTCGAACATGGCag                                                                                                                                                                                                <  2:547800/197‑1 (MQ=255)
                                                                                                                 gTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACACGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGCATTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGCCCAGCATACCTGACAgg                                                                                                                                        >  1:414757/1‑251 (MQ=255)
                                                                                                                               ggAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGAAACGACGAATCGCTTCGGCCAGCCTACCTTACAGGGTCAGAGTAAGCAc                                                                                                                          >  2:391218/1‑251 (MQ=255)
                                                                                                                                        tCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGc                                                                                                                   >  1:566058/1‑249 (MQ=255)
                                                                                                                                         cacaTGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGGGAGCGCATCATCTGTCAACTAAATGGAAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAgg                                                                                                                >  1:30521/1‑250 (MQ=255)
                                                                                                                                         cacaTGTGAGAAAAATTTTTTGTTCACAATAATTTTGCACCCGGGATCTCATCATCTGTAAACAAAGTGGCAGGTGAAGGAGGCATAGGTCATACAAATGAATATTACAGCAAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTGGCTGATACGACGAATCGCTTTGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAg                                                                                                                 <  1:101015/250‑1 (MQ=255)
                                                                                                                                           caTGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGtgg                                                                                                              >  1:308933/1‑250 (MQ=255)
                                                                                                                                               gtgaggacaccatTTTGTTCACAATAATGTTTCATCCGTGAGAGCATCATCTGCCAACTAAATAGCAGGTGAAGGAGGCATAGGTCATCCAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGCTACGACGCATCGCTTGGGCCAGAATACCTGTCCGGGTCAGAGTACGCCCGTTCGGCAGtggttt                                                                                                           >  2:280913/12‑251 (MQ=255)
                                                                                                                                                   gAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTAAACTAAATGGCAGGTGAAGGAGGCATAGGTAATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTAGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTgatttgattt                                                                                                      >  1:478810/1‑251 (MQ=255)
                                                                                                                                                     aaaTATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTc                                                                                                     >  2:229933/1‑250 (MQ=255)
                                                                                                                                                                ttCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACAAAATGGCAGGTGAAGGGGGCATAGGTCATACAAATGGATATTACAGAAAAAAAACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGc                                                                                           <  1:7255/250‑1 (MQ=255)
                                                                                                                                                                      ataatGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAt                                                                                                                                                                                        <  1:201783/151‑1 (MQ=255)
                                                                                                                                                                      ataatGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAt                                                                                                                                                                                        >  2:201783/1‑151 (MQ=255)
                                                                                                                                                                       aaatGTTTCACCCGGGAGCGCATCATCGTACAACAAAATGGCAGGGAAGGGAGGCATAGGGCATAAAAAGGAAAATCACAGACAAAAAACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGt                                                                                   <  2:612736/250‑1 (MQ=255)
                                                                                                                                                                         atGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCCTAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTAGTTGCTGATACGACGAATCGCTTCGACCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGTAt                                                                                 >  2:236042/1‑251 (MQ=255)
                                                                                                                                                                                aTCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGTATCGCAGAc                                                                          >  1:155885/1‑251 (MQ=255)
                                                                                                                                                                                     ggAGCGAATAATCTGTAAACAAAATGGAAGGGGAAGGAGGAAAGGGTCAAACAAAGGGATATAAAAGAAAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGTATCGCAGACAACGa                                                                     <  2:228874/250‑1 (MQ=255)
                                                                                                                                                                                                              tGGCAGGCGACGGCGGCATAGGTCATAAAACTGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGCACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGTATCGCAGACAACGACTTCATCAATTACAGAGTTCCGc                                              >  2:466781/1‑250 (MQ=255)
                                                                                                                                                                                                                 cAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCt                                                                                                                                              >  2:250731/1‑151 (MQ=255)
                                                                                                                                                                                                                 cAGGTGAAGGAGGCATAGGTCATAAAAAGGGATATAACAGAAAAAAAACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCt                                                                                                                                              <  1:250731/151‑1 (MQ=255)
                                                                                                                                                                                                                  aGGTGAAGGAGGCATAGGTCATACAAATGGATATTAAAGACAAAAAAAACCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTATTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGTATCGCAGACAACGACTTCATCAATTACAGAGTTACGCAGgtt                                         >  1:643836/1‑251 (MQ=255)
                                                                                                                                                                                                                    gTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACCACACACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGAATCGCAGACAACGACTTCATCAATTACAGAGTTACGCAGgttgt                                       >  1:119322/1‑251 (MQ=255)
                                                                                                                                                                                                                            aggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGTATCGCAGACAACGACTTCATCAATTACATAGTTACGCAGGTTGTTcgccgc                                >  1:563628/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGTATCGCAGACAACGACTTCATCAATTACAGAGTTACGCAGGTTGTTCGCCGCGTTGcc                          >  1:483542/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aTACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGTATCGCAGACAACGACTTCATCAATTACAGAGTTACGCAGGTTGTTCGCCGCGTTGCCag                        >  1:64432/1‑248 (MQ=255)
                                                                                                                                                                                                                                           aCAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGACTGGTATCGCAGACAACGACTTCATCAATTACAGAGTTACGCAGGTTGTTCGCCGCGTTGCgagaga                     >  2:208295/1‑249 (MQ=255)
                                                                                                                                                                                                                                           aCAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGTATCGCAGACAACGACTTCATCAATTACAGAGTTACGCAGGTTGTTCGCCGCGTTGCCAgaga                     >  1:142156/1‑249 (MQ=255)
                                                                                                                                                                                                                                                tGGATATTAAAGACAAAAAACCCCGACGCAGCGGGTCTTTGAGCCGGGGTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTACGTTGCTGATACGACGACTCGCTTCGGCCAGCATAACTGACAGGGTCAGAGTAAGCAAGTTCGGCAGTGATTTGATTTCATCGCTCAGAGGAATGGTATCGCAGACAACGACTTCATCAATTACAGAGTTACGCAGGTTGTTCGCCGCGGTGCCagagc                     >  1:485412/1‑243 (MQ=255)
                                                                                                                                                                                                                                                gggttatTACACAAAAAAACACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGTATCGCAGACAACGACTTCATCAATTACAGAGTTACGCAGGTTGTTCGCCGCGTTGCCAGAGAAGATCgg              <  2:131037/247‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGTATCGCAGACAACGACTTCATCAATTACAGAGTTACGCAGGTTGTTCGCCGCGTTGCCAGAGAAGATCggg             >  2:215769/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     atTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGTATCGCAGACAACGACTTCATCAATTACAGAGTTACGCAGGTTGTTCGCCGCGTTGCCAGAGAAGATCGGGTg           >  1:429931/1‑249 (MQ=255)
                                                                                                                                                                                                                                                              aCAAAAAACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGTATCGCAGACAACGACTTCATCAATTACAGAGTTACGCAGGTTGTTCGCCGCGTTGCCAGAGAAGATCGGGTGAGTCGCGTa  <  1:580497/251‑1 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                         
CCACTACCAATAGCCAACGCCATCGCCAGCGGGATAGCAATAATCCCGACGGTTATCCCGGCAATCAGGTCACGGGTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAATGGTATCGCAGACAACGACTTCATCAATTACAGAGTTACGCAGGTTGTTCGCCGCGTTGCCAGAGAAGATCGGGTGAGTCGCGTA  >  NC_000913/1260640‑1261137

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