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

AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGC  >  NC_000913/2001769‑2002066
                                                                                                                                                        |                                                                                                                                                    
aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAc                                                                                                                                                     >  2:679450/1‑151 (MQ=255)
 gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcc                                                                                                                                                    <  2:1013116/151‑1 (MQ=255)
   ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcccc                                                                                                                                                  >  1:383473/1‑151 (MQ=255)
       taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgccg                                                                                                                                              >  1:1056767/1‑151 (MQ=255)
             gAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCg                                                                                                                                        >  1:593707/1‑151 (MQ=255)
                  cACGATAAACAGCCCTGCGTTATAGGAGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAACCGGACGATTAGGGGGTAAAATGAGGGGTTATTGGGGGGTAACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGaa                                                                                                                                   <  1:27149/151‑1 (MQ=255)
                  cACGAAAACCAGCCTTGGTTTATAGGTGTTATCGGCATGATTATCCGTTTTTGCAGGTTTTTAAATCGGACGATTAGTGGGTAAAAGGGGGGGTTTTTTGGGGGTGAAAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGaa                                                                                                                                   <  2:366865/151‑1 (MQ=255)
                      aaaaaCAGCCCTGCTTTATATGAGTTATCGGCATGTTTATCCGTTTCTGCGGGGTTTTAAAGCGGACGATTAGTGGGAAAAAGGGGGGTTTATTGGGGGGTTAAAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                               <  1:662296/149‑1 (MQ=255)
                      aTAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACg                                                                                                                              >  2:620188/1‑151 (MQ=255)
                         aaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                            >  1:776543/1‑151 (MQ=255)
                          accgCCCTGTTTTATAAGGTTTATCGGCATAATTATCCGTTTAGGCAGGTTTTTTAATGGGACAATAAGGGGGTGAAAGGAGGGTTTATTGGGGGTTTACAGGAAATCTCCAGCAAAAAAAAACCCCCGCCGGTGGCGGGGAAg                                                                                                                                  <  1:687780/141‑1 (MQ=255)
                          aCAGCCCTGCTTTATATGAGTTAGCGGCATGATTATCCGTTTCTGCGGGTTTTTAAACCGGACGATTGGGGGGTAAAAGGGGGGTTTTTTGGGGGGTTACAGAAAAATCCCAGGAAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                            <  1:1020642/150‑1 (MQ=255)
                              cccTGCTTTATATGAGTTATCGGCATGATTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGGGGGTAAAATGAGGGTTTATTGGGGGTTTACAGTTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGa                                                                                                                       <  2:891235/151‑1 (MQ=255)
                                cTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                                      >  2:572170/1‑150 (MQ=255)
                                    tttATATGGTTAATGGGAATGTTTATCCGTTTCGGGAGGTTTTTAAAGCGGAGTTTTAGTGGTTAAAATGGGGGTTTATTGGGGGTTTACAGTAAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAt                                                                                                                  <  2:699408/149‑1 (MQ=255)
                                     ttATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGTTTAGTGGGTGAAATGGGGGTTTTTTTGGGGGTTACGGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                                 <  2:890071/150‑1 (MQ=255)
                                      tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                                >  2:167909/1‑150 (MQ=255)
                                           gatttATGGGCATGATTATCCGTTTCTGCAGGGTTTTTAATGGGAGGATTAGTGGGTAAAAGGGGGGGTTATTGGGGGGTAAAAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                           <  1:347122/147‑1 (MQ=255)
                                              ttATCGGCATGATTATCCGTTTCTGCAGGTTTTTAAATCGGACGATTAGGGGGTAAAAGGGGGGGTTATTGGGGGGTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                              <  1:20301/144‑1 (MQ=255)
                                                   ggCATGATTATCCTTTTCTGCAGGTTTTTTAATCGGACGATTGGGGGGAAAAAGGGGGGGTTATTGGGGGGTAACAGGAAAATTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                                    <  2:963076/149‑1 (MQ=255)
                                                    gCAGGATTATCCGTTTCTGCGGGTTTTTTAAGCGGACGATGAGTGGGTGAAATGAGGGTTTTTTTGGGGGTTACAGATAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGcc                                                                                                  <  1:617565/150‑1 (MQ=255)
                                                     cATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                                >  1:420250/1‑151 (MQ=255)
                                                        gATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                             >  1:228404/1‑151 (MQ=255)
                                                         aTTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                            >  2:1108746/1‑151 (MQ=255)
                                                           tATCCGTTTCTGCAGGGTTTTTAATCGGACGATTGGGGGGTGAAATGAGGGTTTTTTGGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                           <  1:325642/150‑1 (MQ=255)
                                                           tATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAAgcgg                                                                                          >  1:1146230/1‑151 (MQ=255)
                                                              ccGTTTCGGAGGGTTTTTTAATCGGACGATGAGGGGGAAAAAGGAGGGTTTTTTGGGGGTTAAAAGTAAAATCCCAGAAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                              <  1:1136189/144‑1 (MQ=255)
                                                               cGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                        >  1:91636/1‑149 (MQ=255)
                                                                 tttCTGCGGGTTTTTTAATCGGACGATTAGTGGGTGAAATGGGGGTTTTTTTGGGGGTTACAGGTAATTTCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                                     <  1:13598/150‑1 (MQ=255)
                                                                 tttCTGCAGGTTTTTTAATCGGACGATTGGGGGGTAAAATGGGGGTTTATTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                          <  2:876993/145‑1 (MQ=255)
                                                                  ttcggtgggttttttAAGCGGACTATTAGTGGGAAAAAGGGGGGTTTTTTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaaa                                                                                    <  1:936260/140‑1 (MQ=255)
                                                                   tCTGCGGGTTTTTTAATCGGACTATTAGTGGGTGAAAGGGGGGTTTTTTTGGGGGTTACAGTAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                   <  2:650968/150‑1 (MQ=255)
                                                                         gggTTTTTAATGGGACGATTAGGGGGTGAAATGGGGGGTTATTTGGGGGTTACGGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                                  <  2:184019/145‑1 (MQ=255)
                                                                          gttttttAATCGGACGATTAGGGGGTGAAATGGGGGGTTTTTGGGGGGTTACAGTTAAATTCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                           <  2:808479/149‑1 (MQ=255)
                                                                           gTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                          >  1:313971/1‑151 (MQ=255)
                                                                             ttttAATGGGCGGATTAGGGGGTAAAAGGAGGTTTTATTGGGGGGTTAAAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                         <  2:671938/150‑1 (MQ=37)
                                                                              tttAATCGGACGATTAGGGGGTAAAAGGGGGGTTTTTTGGGGGTTTACGGGTAATTCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                        <  2:542709/150‑1 (MQ=37)
                                                                                       aggATTGGGGGGTAAAATGGGGGTTTATTGGGGGGTTAAAGGTAATTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                               <  1:115843/148‑1 (MQ=25)
                                                                                       aCGATTAGTGGGTAAAATGGGGGTTTATTTGGGGTTTACAGGAAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                               <  2:689257/150‑1 (MQ=37)
                                                                                        cGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                             >  1:361614/1‑151 (MQ=255)
                                                                                           ttAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                          >  1:1123641/1‑151 (MQ=255)
                                                                                              gggggTAAAATGGGGGTTTTTTGGGGGTTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                         <  1:870164/147‑1 (MQ=37)
                                                                                                gggTAAAATGGGGGGTTATTGGGGGGTTAAAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                         <  1:227681/147‑1 (MQ=37)
                                                                                                      aaagggggtgtttttGGGGTTTTACAGGTAAATCCCAGACAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                                <  1:984608/139‑1 (MQ=25)
                                                                                                      aaaggggggTTTTTTGGGGGGTTAAAGTAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                  <  1:832586/144‑1 (MQ=25)
                                                                                                         gggggtTTATTTGGGGTTTACAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaaa                                             <  2:585191/148‑1 (MQ=37)
                                                                                                         gAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                            >  2:219297/1‑151 (MQ=255)
                                                                                                           gggtTTATTGGGGGGTTACAGGAAAATTCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                          <  1:974368/151‑1 (MQ=37)
                                                                                                             ggTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                         >  2:112353/1‑150 (MQ=255)
                                                                                                                tATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                     >  2:363039/1‑151 (MQ=255)
                                                                                                                tATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                     >  2:1032061/1‑151 (MQ=255)
                                                                                                                       ggggTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                               >  2:305836/1‑151 (MQ=255)
                                                                                                                       ggggTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                               >  2:1041776/1‑151 (MQ=255)
                                                                                                                                 ggTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt                       >  1:320492/1‑150 (MQ=255)
                                                                                                                                 ggTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg                      >  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)
                                                                                                                                                        |                                                                                                                                                    
AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGC  >  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: