Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,229,2040TC44.4% 62.5 / 25.7 63intergenic (+139/+234)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base T (14/21);  new base C (28/0);  total (42/21)
Fisher's exact test for biased strand distribution p-value = 8.80e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAC  >  NC_000913/2229056‑2229351
                                                                                                                                                     |                                                                                                                                                    
gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGGGGGCAAACACCCCCCCCCCCCGc                                                                                                                                                    >  1:334667/1‑151 (MQ=255)
  tatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGCCACAAGGCGCCGCCCCCTCCAc                                                                                                                                                  >  1:793176/1‑151 (MQ=255)
   atTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCCGCAcc                                                                                                                                                 >  2:830135/1‑151 (MQ=255)
     tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCCCCCCCCGCAcccc                                                                                                                                               >  1:386379/1‑151 (MQ=255)
            ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGGCAATCGCCGCCGCCCCCCCCCCCCCCCGGctc                                                                                                                                        >  1:1271226/1‑151 (MQ=255)
                ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCATCGCCGCCGCCCCCCGCACCCCCGGGCTCTgg                                                                                                                                     >  2:318378/1‑150 (MQ=255)
                 taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGGCAAACGCCCCCCCCCCCCGCCCCCCCCGGCTCCGGCg                                                                                                                                   >  1:910918/1‑151 (MQ=37)
                   atGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCCCCGCCCCCCCCACCCCCGGGCCCTGGcca                                                                                                                                  >  2:807469/1‑148 (MQ=255)
                    tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGGCAAATGCCCCCGGCCCCCGCACCCCCGGGCCCTTGCg                                                                                                                                   >  1:451155/1‑148 (MQ=37)
                     gTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGGCAAACGCCGCCGCCCCCCCCACCCCCCGGCCCCGGCgga                                                                                                                                 >  1:524685/1‑147 (MQ=37)
                      tcatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCCCCGCCCCCCGCACCCCCCGGCTCTGGCGAAc                                                                                                                                >  2:731685/1‑149 (MQ=255)
                       catcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCAATCGCCCCCGCCCCCCGCACCCCCCGGCTCTGGCGAACaa                                                                                                                              >  1:282926/1‑150 (MQ=255)
                          caGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATACAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATc                                                                                                                          <  2:894534/151‑1 (MQ=255)
                           aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCCCCCCCCCCTGCACCCCCCGGCCCTGGCGAACaaaa                                                                                                                            >  2:711670/1‑148 (MQ=255)
                               ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgccgc                                                                                                                     <  1:1183688/151‑1 (MQ=255)
                                cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGAGGGGGCGGCCAATGGCCGCCGCCCCCTGCACCCCCGGGCTCCGGCGAACAAAATcggc                                                                                                                       >  2:392250/1‑146 (MQ=255)
                                 gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCCGCACCCCCCGGCTCCGGCGAACAAAAAcgccgc                                                                                                                     >  2:1171838/1‑149 (MQ=255)
                                  aaaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGGCAATCGCCGCCGCCCCCCGCACCCCCCCGCTCTTGCGAACAAAAACGCCGCt                                                                                                                    >  1:1003851/1‑149 (MQ=25)
                                  aaaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTgc                                                                                                                  <  1:1243891/151‑1 (MQ=255)
                                   aaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCACACGCCGCCGCCCCCCTCACCCCCCGGCGCTCGCGAACaaaaaaaccgccgcg                                                                                                                 >  1:1032605/1‑140 (MQ=21)
                                    aaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAAGGGGCGGGCCAACCCCGCCGCCCCCCCCACCCCCCCGCCCTGGCCGACaaaaacccccc                                                                                                                     >  1:1026894/1‑139 (MQ=255)
                                    aaTAAAGAAAATGTTTTAGTCTCCCATCCTGCTGGGCTCCTTTCAAGTGGATTGCTTCATATAATCACTGGCGCAAGGAGGGCGCCAGGCGGCGGCCAACCACCCCCCCCCCCGCGCCCCCCGGCACTGCCGCACaaa                                                                                                                             >  2:601198/1‑138 (MQ=255)
                                        aCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGCCAAACGCCGCCGCCCCCCCCACCCCCCGGCCCCGGCGAACAAAAAAGCCCCCGCGCgg                                                                                                              >  1:806936/1‑149 (MQ=255)
                                        aCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTg                                                                                                            <  2:386365/151‑1 (MQ=255)
                                         cGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGCGGCAAACCCCGCCGCCCCCCCCCCCCCCCGGCCCCGGCGCACAAAATCGCCGCCGCGCgg                                                                                                              >  1:562190/1‑148 (MQ=21)
                                              atatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGCGGGCAATCGCCGCCGCCCCCTGCACCCCCCGGCTCTGGCGGACAAAAAAGCCGCTGCGCGGTGCCCTc                                                                                                       >  2:1117004/1‑150 (MQ=21)
                                               tatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCAAACGCCCCCCCCCCCCCCACCCCCCGGCTCTGGCGGACAAAAACCCCGCCGCGCGGGGccccc                                                                                                       >  1:344274/1‑147 (MQ=18)
                                                    cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGCCAAACGCCCCCGCCCCCCCCACCCCCCGGCCCTTGCGAAAAAAAAAACCCCCGCGCGGGGCCCCCCGCCTAt                                                                                                >  1:1153834/1‑151 (MQ=255)
                                                      gTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGGCAATCGCCGCCGCCCCCCGCACCCCCCGGCTCTGGCGAACAAAAACGCCCCCGCGCGGGGCCCCCCGCCTATcc                                                                                              >  2:950899/1‑151 (MQ=18)
                                                      gTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCCCCACCCCCCGGCTCCGGCGAAAAAAAACGCCGCCGCGCGGGGGCCTCCGCCTTTcc                                                                                              >  2:856477/1‑151 (MQ=18)
                                                           tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCAGGGGCGGGCCATCGCCGCCCCCCCCCGCACCCCCCGGGTCTGGCGGACAAAATTTCCCCTGCGCGGGGccccccgcgcattccct                                                                                           >  2:1012537/1‑135 (MQ=255)
                                                            cATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                                                           <  1:778394/148‑1 (MQ=255)
                                                              tCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACgg                                                                                       <  2:137671/150‑1 (MQ=255)
                                                                     gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                                               <  2:505968/150‑1 (MQ=255)
                                                                     gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                                               <  1:759729/150‑1 (MQ=255)
                                                                      cgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCCGGCTCTGGCGAACAAAAACGCCGCTGCGCGGGGCCCCCGGCTTATCCCTTTCGGCTaaaccgg                                                                              >  2:903884/1‑145 (MQ=17)
                                                                          ccTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCgg                                                                          <  2:578785/150‑1 (MQ=255)
                                                                            tGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCCCCGCCCCCCGCACCCCCCCGCCCTGGCGGACAAAAACCCCGCTGCGCCGGGCCCCCCGCCTATCCCCTTCGGCTAccccggccgggg                                                                       >  1:1015551/1‑141 (MQ=17)
                                                                             gTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgcg                                                                       <  2:644657/150‑1 (MQ=255)
                                                                                  gtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                                                  <  1:988925/150‑1 (MQ=255)
                                                                                   tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCATCGCCGCCGCCCCCCCCACCCCCGGGCTCTGGCGGACAAAAACGCCGCTGCGCGGTGCCCCCCGCTTATCCCCTCCCGCTACCCCGGCGGGcgcggg                                                                   >  2:822938/1‑146 (MQ=14)
                                                                                              atatAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTa                                                     <  1:3440/151‑1 (MQ=37)
                                                                                               tataATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCATCCCCCCCCCCCCCCGCACCCCCCGGCTCTGGCGGACAAAATCGCCGCTGCGCGGGGCCCTCCGCCCATCCCTTACGGCTACCCGGGCCGGCGCGCGGGAAAAccccccta                                                     >  2:749654/1‑146 (MQ=2)
                                                                                                     cACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCCGCACCCCCGGGCTCTTGCGAACAAAAACGCCGCTGCGCGGTGCCCTCCACTTATCCCTTACGGCTAACCGGGCGGGGGCGCGGGAAAAAcccccttaaaa                                                 >  2:134091/1‑140 (MQ=11)
                                                                                                     cACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCCCCGCCCCCCCCACCCCCCGGCCCTTGCGAAAAAAAACCCCCCCGCGCGGGGCCCTCCGCCTATCCCCTTCCGCTACCCGGGTGGGGGCGGGGGAAAAATCCCTTTaaaa                                                  >  1:517402/1‑148 (MQ=2)
                                                                                                            gcAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg                                       <  2:815854/151‑1 (MQ=31)
                                                                                                             cAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc                                      <  2:339032/151‑1 (MQ=31)
                                                                                                                          aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt                          <  2:15524/151‑1 (MQ=31)
                                                                                                                             gggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTggg                       >  2:667405/1‑151 (MQ=31)
                                                                                                                              ggcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGc                      >  1:195824/1‑151 (MQ=31)
                                                                                                                               gcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCt                     <  2:102684/151‑1 (MQ=31)
                                                                                                                               gcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCt                     <  2:1062278/151‑1 (MQ=31)
                                                                                                                                   ccAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGccc                 <  1:419166/151‑1 (MQ=31)
                                                                                                                                    cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc                >  2:1025940/1‑151 (MQ=31)
                                                                                                                                    cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGGAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc                >  1:605742/1‑151 (MQ=14)
                                                                                                                                     aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGGCGcccc                >  2:1035437/1‑150 (MQ=14)
                                                                                                                                      aTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCgg              <  2:334653/151‑1 (MQ=31)
                                                                                                                                              cgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGGGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAgg      >  2:188629/1‑151 (MQ=14)
                                                                                                                                               gccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTc         >  1:1100368/1‑147 (MQ=31)
                                                                                                                                               gccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAggg     <  2:414517/151‑1 (MQ=31)
                                                                                                                                               gccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAggg     <  1:1238583/151‑1 (MQ=31)
                                                                                                                                               gccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGGAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAggg     >  2:742814/1‑151 (MQ=14)
                                                                                                                                                  cccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAc  >  1:980300/1‑151 (MQ=31)
                                                                                                                                                     |                                                                                                                                                    
GATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAC  >  NC_000913/2229056‑2229351

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