Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9180AC44.3% 97.3 / 19.7 61intergenic (‑129/+192)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (29/5);  new base C (0/27);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 2.02e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGC  >  NC_000913/2001769‑2002066
                                                                                                                                                        |                                                                                                                                                    
aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAAc                                                                                                                                                     >  2:679450/1‑151 (MQ=255)
 gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAAcc                                                                                                                                                    <  2:1013116/151‑1 (MQ=255)
   ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAAcccc                                                                                                                                                  >  1:383473/1‑151 (MQ=255)
       taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCccgccg                                                                                                                                              >  1:1056767/1‑151 (MQ=255)
             gAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCg                                                                                                                                        >  1:593707/1‑151 (MQ=255)
                  cACGATAAACAGCCCTGCGTTATAGGAGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAACCGGACGATTAGGGGGTAAAATGAGGGGTTATTG‑GGGGG‑TAACAGGT‑AAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGaa                                                                                                                                   <  1:27149/151‑1 (MQ=255)
                  cACGAAAACCAGCCTTGGTTTATAGGTGTTATCGGCATGATTATCCGTTTTTGCAGGTTTTTAAATCGGACGATTAGTGGGTAAAAGGGGGGGTTTTTT‑GGGGG‑TGAAAGGT‑AAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGaa                                                                                                                                   <  2:366865/151‑1 (MQ=255)
                      aaaaaCAGCCCTGCTTTATATGAGTTATCGGCATGTTTATCCGTTTCTGCGGGGTTTTAAAGCGGACGATTAGTGGGAAAAAGGGGGGTTTATTG‑GGGGG‑TTAAAGGA‑AAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                               <  1:662296/149‑1 (MQ=255)
                      aTAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGT‑AAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACg                                                                                                                              >  2:620188/1‑151 (MQ=255)
                         aaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                            >  1:776543/1‑151 (MQ=255)
                          accgCCCTGTTTTATAAGGTTTATCGGCATAATTATCCGTTTAGGCAGGTTTTTTAATGGGACAATAAGGGGGTGAAAGGAGGGTTTATTG‑GGGGT‑TTACAGGA‑AATCTCCAGCAAAAAAAAACCCCCGCCGGTGGCGGGGAAg                                                                                                                                  <  1:687780/141‑1 (MQ=255)
                          aCAGCCCTGCTTTATATGAGTTAGCGGCATGATTATCCGTTTCTGCGGGTTTTTAAACCGGACGATTGGGGGGTAAAAGGGGGGTTTTTTG‑GGGGG‑TTACAGAA‑AAATCCCAGGAAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                            <  1:1020642/150‑1 (MQ=255)
                              cccTGCTTTATATGAGTTATCGGCATGATTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGGGGGTAAAATGAGGGTTTATTG‑GGGGT‑TTACAGTT‑AAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGa                                                                                                                       <  2:891235/151‑1 (MQ=255)
                                cTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                                      >  2:572170/1‑150 (MQ=255)
                                    tttATATGGTTAATGGGAATGTTTATCCGTTTCGGGAGGTTTTTAAAGCGGAGTTTTAGTGGTTAAAATGGGGGTTTATTG‑GGGGT‑TTACAGTA‑AAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAt                                                                                                                  <  2:699408/149‑1 (MQ=255)
                                     ttATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGTTTAGTGGGTGAAATGGGGGTTTTTTT‑GGGGG‑TTACGGGT‑AAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                                 <  2:890071/150‑1 (MQ=255)
                                      tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                                >  2:167909/1‑150 (MQ=255)
                                           gatttATGGGCATGATTATCCGTTTCTGCAGGGTTTTTAATGGGAGGATTAGTGGGTAAAAGGGGGGGTTATTG‑GGGGG‑TAAAAGGA‑AAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                           <  1:347122/147‑1 (MQ=255)
                                              ttATCGGCATGATTATCCGTTTCTGCAGGTTTTTAAATCGGACGATTAGGGGGTAAAAGGGGGGGTTATTG‑GGGGG‑TTACAGGA‑AAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                              <  1:20301/144‑1 (MQ=255)
                                                   ggCATGATTATCCTTTTCTGCAGGTTTTTTAATCGGACGATTGGGGGGAAAAAGGGGGGGTTATTG‑GGGGG‑TAACAGGA‑AAATTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                                    <  2:963076/149‑1 (MQ=255)
                                                    gCAGGATTATCCGTTTCTGCGGGTTTTTTAAGCGGACGATGAGTGGGTGAAATGAGGGTTTTTTT‑GGGGG‑TTACAGAT‑AAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGcc                                                                                                  <  1:617565/150‑1 (MQ=255)
                                                     cATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                                >  1:420250/1‑151 (MQ=255)
                                                        gATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                             >  1:228404/1‑151 (MQ=255)
                                                         aTTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                            >  2:1108746/1‑151 (MQ=255)
                                                           tATCCGTTTCTGCAGGGTTTTTAATCGGACGATTGGGGGGTGAAATGAGGGTTTTTTG‑GGGGG‑TTACAGGT‑AAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                           <  1:325642/150‑1 (MQ=255)
                                                           tATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAAgcgg                                                                                          >  1:1146230/1‑151 (MQ=255)
                                                              ccGTTTCGGAGGGTTTTTTAATCGGACGATGAGGGGGAAAAAGGA‑GGGTTTTTT‑GGGGG‑TTAAAAGTAAAATCCCAGAAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                              <  1:1136189/144‑1 (MQ=255)
                                                               cGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                        >  1:91636/1‑149 (MQ=255)
                                                                 tttCTGCGGGTTTTTTAATCGGACGATTAGTGGGTGAAATGGGGGTTTTTTT‑GGGGG‑TTACAGGT‑AATTTCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                                     <  1:13598/150‑1 (MQ=255)
                                                                 tttCTGCAGGTTTTTTAATCGGACGATTGGGGGGTAAAATGGGGGTTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                          <  2:876993/145‑1 (MQ=255)
                                                                  ttcggtgggttttttAAGCGGACTATTAGTGGGAAAAAGGGGGGTTTTTTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaaa                                                                                    <  1:936260/140‑1 (MQ=255)
                                                                   tCTGCGGGTTTTTTAATCGGACTATTAGTGGGTGAAAGGGGGGTTTTTTT‑GGGGG‑TTACAGTA‑AAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                   <  2:650968/150‑1 (MQ=255)
                                                                         gggTTTTTAATGGGACGATTAGGGGGTGAAATGGGGGGTTATTT‑GGGGG‑TTACGGGA‑AAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                                  <  2:184019/145‑1 (MQ=255)
                                                                          gttttttAATCGGACGATTAGGGGGTGAAATGGGGGGTTTTTG‑GGGGG‑TTACAGTT‑AAATTCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                           <  2:808479/149‑1 (MQ=255)
                                                                           gTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                          >  1:313971/1‑151 (MQ=255)
                                                                             ttttAATGGGCGGATTAGGGGGTAAAAGGAGGTTTTATTG‑GGGGG‑TTAAAGGT‑AAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                         <  2:671938/150‑1 (MQ=37)
                                                                              tttAATCGGACGATTAGGGGGT‑AAAAGGGGGGTTTTTT‑GGGGGTTTACGGGT‑AATTCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                        <  2:542709/150‑1 (MQ=37)
                                                                                       aggATTGGGGGGTAAAATGGGGGTTTATTG‑GGGGG‑TTAAAGGT‑AATTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                               <  1:115843/148‑1 (MQ=25)
                                                                                       aCGATTAGTGGGTAAAATGGGGGTTTATTT‑GGGGT‑TTACAGGA‑AAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                               <  2:689257/150‑1 (MQ=37)
                                                                                        cGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                             >  1:361614/1‑151 (MQ=255)
                                                                                           ttAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                          >  1:1123641/1‑151 (MQ=255)
                                                                                              gggggTAAAATG‑GGGGTTTTTT‑GGGGGTTTACAGGT‑AAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                         <  1:870164/147‑1 (MQ=37)
                                                                                                gggTAAAATGGGGGGTTATTG‑GGGGG‑TTAAAGGA‑AAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                         <  1:227681/147‑1 (MQ=37)
                                                                                                      aaagggggtgttttt‑GGGGTTTTACAGGT‑AAATCCCAGACAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                                <  1:984608/139‑1 (MQ=25)
                                                                                                      aaaggggggTTTTTTGGGGGG‑TTAAAGTA‑AAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                  <  1:832586/144‑1 (MQ=25)
                                                                                                         gggggtTTATTT‑GGGGT‑TTACAGGA‑AAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaaa                                             <  2:585191/148‑1 (MQ=37)
                                                                                                         gAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                            >  2:219297/1‑151 (MQ=255)
                                                                                                           gggtTTATTG‑GGGGG‑TTACAGGA‑AAATTCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                          <  1:974368/151‑1 (MQ=37)
                                                                                                             ggTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                         >  2:112353/1‑150 (MQ=255)
                                                                                                                tATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                     >  2:363039/1‑151 (MQ=255)
                                                                                                                tATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                     >  2:1032061/1‑151 (MQ=255)
                                                                                                                       gggg‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                               >  2:305836/1‑151 (MQ=255)
                                                                                                                       gggg‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                               >  2:1041776/1‑151 (MQ=255)
                                                                                                                                 ggT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt                       >  1:320492/1‑150 (MQ=255)
                                                                                                                                 ggT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg                      >  2:291792/1‑151 (MQ=255)
                                                                                                                                      aaTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa                  <  1:809183/151‑1 (MQ=255)
                                                                                                                                          ccAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc              >  1:963792/1‑151 (MQ=255)
                                                                                                                                           cAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc              >  2:866043/1‑150 (MQ=255)
                                                                                                                                             ggaaaaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca             <  1:963800/144‑1 (MQ=255)
                                                                                                                                               cAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc         >  2:104222/1‑151 (MQ=255)
                                                                                                                                               cAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc         >  1:699339/1‑151 (MQ=255)
                                                                                                                                                aaaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc         <  2:1070179/148‑1 (MQ=255)
                                                                                                                                                aGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc         >  2:1168781/1‑150 (MQ=255)
                                                                                                                                                 aaaaaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt       <  2:383474/150‑1 (MQ=255)
                                                                                                                                                    aaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc      >  1:348908/1‑149 (MQ=255)
                                                                                                                                                     aaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg   <  1:519721/151‑1 (MQ=255)
                                                                                                                                                       aaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc  >  2:98717/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGT‑AAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGC  >  NC_000913/2001769‑2002066

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