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

GATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAC  >  NC_000913/2229056‑2229351
                                                                                                                                                     |                                                                                                                                                    
gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CGGGGGGGGGGCAAACACCCCCCCCCCCCGc                                                                                                                                                    >  1:334667/1‑151 (MQ=255)
  tatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CCGGGGGGGGCCACAAGGCGCCGCCCCCTCCAc                                                                                                                                                  >  1:793176/1‑151 (MQ=255)
   atTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCCGCAcc                                                                                                                                                 >  2:830135/1‑151 (MQ=255)
     tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCCCCCCCCGCAcccc                                                                                                                                               >  1:386379/1‑151 (MQ=255)
            ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGGCAATCGCCGCCGCCCCCCCCCCCCCCCGGctc                                                                                                                                        >  1:1271226/1‑151 (MQ=255)
                ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCCATCGCCGCCGCCCCCCGCACCCCCGGGCTCTgg                                                                                                                                     >  2:318378/1‑150 (MQ=255)
                 taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGGGGGCAAACGCCCCCCCCCCCCGCCCCCCCCGGCTCCGGCg                                                                                                                                   >  1:910918/1‑151 (MQ=37)
                   atGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCCCCGCCCCCCCCACCCCCGGGCCCTGGcca                                                                                                                                  >  2:807469/1‑148 (MQ=255)
                    tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGGGGGCAAATGCCCCCGGCCCCCGCACCCCCGGGCCCTTGCg                                                                                                                                   >  1:451155/1‑148 (MQ=37)
                     gTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGGGGGCAAACGCCGCCGCCCCCCCCACCCCCCGGCCCCGGCgga                                                                                                                                 >  1:524685/1‑147 (MQ=37)
                      tcatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCCCCGCCCCCCGCACCCCCCGGCTCTGGCGAAc                                                                                                                                >  2:731685/1‑149 (MQ=255)
                       catcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGGGGCCAATCGCCCCCGCCCCCCGCACCCCCCGGCTCTGGCGAACaa                                                                                                                              >  1:282926/1‑150 (MQ=255)
                          caGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATACAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATc                                                                                                                          <  2:894534/151‑1 (MQ=255)
                           aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCCCCCCCCCCTGCACCCCCCGGCCCTGGCGAACaaaa                                                                                                                            >  2:711670/1‑148 (MQ=255)
                               ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgccgc                                                                                                                     <  1:1183688/151‑1 (MQ=255)
                                cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑GAGGGGGCGGCCAATGGCCGCCGCCCCCTGCACCCCCGGGCTCCGGCGAACAAAATcggc                                                                                                                       >  2:392250/1‑146 (MQ=255)
                                 gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCCGCACCCCCCGGCTCCGGCGAACAAAAAcgccgc                                                                                                                     >  2:1171838/1‑149 (MQ=255)
                                  aaaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGGCAATCGCCGCCGCCCCCCGCACCCCCCCGCTCTTGCGAACAAAAACGCCGCt                                                                                                                    >  1:1003851/1‑149 (MQ=25)
                                  aaaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTgc                                                                                                                  <  1:1243891/151‑1 (MQ=255)
                                   aaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGGGGCCACACGCCGCCGCCCCCCTCACCCCCCGGCGCTCGCGAACaaaaaaaccgccgcg                                                                                                                 >  1:1032605/1‑140 (MQ=21)
                                    aaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAAGGGGCGGGCCAACCCCGCCGCCCCCCCCACCCCCCCGCCCTGGCCGACaaaaacccccc                                                                                                                     >  1:1026894/1‑139 (MQ=255)
                                    aaTAAAGAAAATGTTTTAGTCTCCCATCCTGCTGGGCTCCTTTCAAGTGGATTGCTTCATATAATCACTGGCGCAAGGAGGGCGCCAGGCGGCGGCCAACCACC‑CCCCCCCCCGCGCCCCCCGGCACTGCCGCACaaa                                                                                                                             >  2:601198/1‑138 (MQ=255)
                                        aCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CCGGGGGGGGCCAAACGCCGCCGCCCCCCCCACCCCCCGGCCCCGGCGAACAAAAAAGCCCCCGCGCgg                                                                                                              >  1:806936/1‑149 (MQ=255)
                                        aCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTg                                                                                                            <  2:386365/151‑1 (MQ=255)
                                         cGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGGCGGCAAACCCCGCCGCCCCCCCCCCCCCCCGGCCCCGGCGCACAAAATCGCCGCCGCGCgg                                                                                                              >  1:562190/1‑148 (MQ=21)
                                              atatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CGGGGGGCGGGCAATCGCCGCCGCCCCCTGCACCCCCCGGCTCTGGCGGACAAAAAAGCCGCTGCGCGGTGCCCTc                                                                                                       >  2:1117004/1‑150 (MQ=21)
                                               tatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGGGGCCAAACGCCCCCCCCCCCCCCACCCCCCGGCTCTGGCGGACAAAAACCCCGCCGCGCGGGGccccc                                                                                                       >  1:344274/1‑147 (MQ=18)
                                                    cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CCGGGGGGGGCCAAACGCCCCCGCCCCCCCCACCCCCCGGCCCTTGCGAAAAAAAAAACCCCCGCGCGGGGCCCCCCGCCTAt                                                                                                >  1:1153834/1‑151 (MQ=255)
                                                      gTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGGCAATCGCCGCCGCCCCCCGCACCCCCCGGCTCTGGCGAACAAAAACGCCCCCGCGCGGGGCCCCCCGCCTATcc                                                                                              >  2:950899/1‑151 (MQ=18)
                                                      gTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCCCCACCCCCCGGCTCCGGCGAAAAAAAACGCCGCCGCGCGGGGGCCTCCGCCTTTcc                                                                                              >  2:856477/1‑151 (MQ=18)
                                                           tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CCAGGGGCGGGCCATCGCCGCCCCCCCCCGCACCCCCCGGGTCTGGCGGACAAAATTTCCCCTGCGCGGGGccccccgcgcattccct                                                                                           >  2:1012537/1‑135 (MQ=255)
                                                            cATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                                                           <  1:778394/148‑1 (MQ=255)
                                                              tCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACgg                                                                                       <  2:137671/150‑1 (MQ=255)
                                                                     gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCgg                                                                               <  2:505968/150‑1 (MQ=255)
                                                                     gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCgg                                                                               <  1:759729/150‑1 (MQ=255)
                                                                      cgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCCGGCTCTGGCGAACAAAAACGCCGCTGCGCGGGGCCCCCGGCTTATCCCTTTCGGCTaaaccgg                                                                              >  2:903884/1‑145 (MQ=17)
                                                                          ccTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCgg                                                                          <  2:578785/150‑1 (MQ=255)
                                                                            tGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCCCCGCCCCCCGCACCCCCCCGCCCTGGCGGACAAAAACCCCGCTGCGCCGGGCCCCCCGCCTATCCCCTTCGGCTA‑ccccggccgggg                                                                       >  1:1015551/1‑141 (MQ=17)
                                                                             gTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGgcg                                                                       <  2:644657/150‑1 (MQ=255)
                                                                                  gtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAgg                                                                  <  1:988925/150‑1 (MQ=255)
                                                                                   tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCCATCGCCGCCGCCCCCCCCACCCCCGGGCTCTGGCGGACAAAAACGCCGCTGCGCGGTGCCCCCCGCTTATCCCCTCCCGCTA‑CCCCGGCGGGcgcggg                                                                   >  2:822938/1‑146 (MQ=14)
                                                                                              atatAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTa                                                     <  1:3440/151‑1 (MQ=37)
                                                                                               tataATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCCATCCCCCCCCCCCCCCGCACCCCCCGGCTCTGGCGGACAAAATCGCCGCTGCGCGGGGCCCTCCGCCCATCCCTTACGGCTA‑CCCGGGCCGGCGCGCGGGAAAAccccccta                                                     >  2:749654/1‑146 (MQ=2)
                                                                                                     cACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCCGCACCCCCGGGCTCTTGCGAACAAAAACGCCGCTGCGCGGTGCCCTCCACTTATCCCTTACGGCTA‑ACCGGGCGGGGGCGCGGGAAAAAcccccttaaaa                                                 >  2:134091/1‑140 (MQ=11)
                                                                                                     cACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCCCCGCCCCCCCCACCCCCCGGCCCTTGCGAAAAAAAACCCCCCCGCGCGGGGCCCTCCGCCTATCCCCTTCCGCTACCCGGGTGGGGGCGGGGGAAAAATCCCTTTaaaa                                                  >  1:517402/1‑148 (MQ=2)
                                                                                                            gcAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg                                       <  2:815854/151‑1 (MQ=31)
                                                                                                             cAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc                                      <  2:339032/151‑1 (MQ=31)
                                                                                                                          aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt                          <  2:15524/151‑1 (MQ=31)
                                                                                                                             gggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTggg                       >  2:667405/1‑151 (MQ=31)
                                                                                                                              ggcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGc                      >  1:195824/1‑151 (MQ=31)
                                                                                                                               gcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCt                     <  2:102684/151‑1 (MQ=31)
                                                                                                                               gcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCt                     <  2:1062278/151‑1 (MQ=31)
                                                                                                                                   ccAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGccc                 <  1:419166/151‑1 (MQ=31)
                                                                                                                                    cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc                >  2:1025940/1‑151 (MQ=31)
                                                                                                                                    cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGGAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc                >  1:605742/1‑151 (MQ=14)
                                                                                                                                     aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGGCGcccc                >  2:1035437/1‑150 (MQ=14)
                                                                                                                                      aTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCgg              <  2:334653/151‑1 (MQ=31)
                                                                                                                                              cgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGGGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAgg      >  2:188629/1‑151 (MQ=14)
                                                                                                                                               gccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTc         >  1:1100368/1‑147 (MQ=31)
                                                                                                                                               gccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAggg     <  2:414517/151‑1 (MQ=31)
                                                                                                                                               gccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAggg     <  1:1238583/151‑1 (MQ=31)
                                                                                                                                               gccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGGAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAggg     >  2:742814/1‑151 (MQ=14)
                                                                                                                                                  cccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAc  >  1:980300/1‑151 (MQ=31)
                                                                                                                                                     |                                                                                                                                                    
GATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTA‑CCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAC  >  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: ‑