Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,552,1910CA25.9% 90.5 / 18.2 58intergenic (+200/+207)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (14/29);  new base A (9/6);  total (23/35)
Fisher's exact test for biased strand distribution p-value = 7.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAG  >  NC_000913/1552042‑1552335
                                                                                                                                                     |                                                                                                                                                
atatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCTTGAACGCCTTATCCGGCATGAAAAACCTTCAAATCCATTAGTTTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTc                                                                                                                                                  <  1:642451/150‑1 (MQ=9)
atatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTc                                                                                                                                                  <  2:1363512/150‑1 (MQ=31)
 tatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCa                                                                                                                                                 <  1:1289303/150‑1 (MQ=32)
 tatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCa                                                                                                                                                 <  2:584131/150‑1 (MQ=32)
 tatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCa                                                                                                                                                 <  2:1079501/150‑1 (MQ=32)
      tAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCATTGAACGTTTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATgg                                                                                                                                            <  1:195729/150‑1 (MQ=21)
         aCAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcgc                                                                                                                                         <  1:1006441/150‑1 (MQ=34)
                aGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGAATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGGGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGGTTGTCATCAGTTTCAAATGGCGCTGTaaaa                                                                                                                                  <  1:1411531/150‑1 (MQ=255)
                 ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAg                                                                                                                                 >  1:714861/1‑150 (MQ=17)
                  gCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAgg                                                                                                                                >  1:1364180/1‑150 (MQ=40)
                  gCTTTGTGCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCTAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCTCTTCAGGAAATGTTGCATTTGTCATCAGTTTCAAATGGCGCTGTAAAAgg                                                                                                                                <  2:916876/150‑1 (MQ=14)
                   gtttGTTCATGCCGGATCGGGCTGTAACG‑CATATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                                                                                                                              <  2:1208457/149‑1 (MQ=14)
                    tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                                                                                                                              <  2:1098397/150‑1 (MQ=255)
                    tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                                                                                                                              <  2:1567045/150‑1 (MQ=255)
                      tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTaaaaagagac                                                                                                                            >  2:1583155/1‑144 (MQ=255)
                       gTTCATGCCGGATGCGGAGTGAACGCCTTATCGGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGAAATGTTTCGTTTGTCATCAGTTTCAATTGGCGCTGTAAAAGGCGTCa                                                                                                                           <  2:872869/150‑1 (MQ=14)
                                tgATGAGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCATATCCAATAGATTGCAGTGAACGTGTAGTCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                                                                                                  <  1:1570028/149‑1 (MQ=255)
                                ggATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                                                                                                  <  1:407016/150‑1 (MQ=255)
                                  aTGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTg                                                                                                                <  2:1567989/150‑1 (MQ=255)
                                    gcggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGTGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa                                                                                                              <  2:415738/150‑1 (MQ=255)
                                      ggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGa                                                                                                            >  2:1028793/1‑150 (MQ=255)
                                       gcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGAc                                                                                                           <  1:1409517/150‑1 (MQ=255)
                                        cgTTAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGATTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGAct                                                                                                          <  1:492392/150‑2 (MQ=255)
                                        cgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACa                                                                                                          <  1:1345059/150‑1 (MQ=255)
                                        cgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACa                                                                                                          <  1:1351108/150‑1 (MQ=255)
                                             aCGCCTTATCCGTAATGAAAACACTTCAACTCCAATAGATTGCAGTGAACGTGTAGTCCTGATAAGCGTAGCGCATCAGGAAATGTTGCGTTTTTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTa                                                                                                     <  2:255711/150‑1 (MQ=18)
                                              cGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAACGAGACATTTTCATATTGTATACAAAGTAg                                                                                                    >  2:1108133/1‑150 (MQ=14)
                                                    aTCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTg                                                                                              <  1:76670/150‑1 (MQ=255)
                                                     tCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGt                                                                                             <  2:1336177/150‑1 (MQ=255)
                                                       cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAGGAGACATTTTGATATTGTTGACAAAGAAGGCTTTGtaa                                                                                           >  1:955297/1‑148 (MQ=14)
                                                       cGGAATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAAGGAGTCATTTTGATATTTTATACAAAGTCTGCTTTGTTc                                                                                           >  1:955300/1‑150 (MQ=14)
                                                        ggCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGAAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCa                                                                                          <  1:802904/150‑1 (MQ=255)
                                                          cATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGAGTCATTTTATTATTGTATAAAAAGAAGGCTTTATACAta                                                                                        >  1:166378/1‑149 (MQ=14)
                                                            tGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGcc                                                                                      <  1:425353/150‑1 (MQ=255)
                                                            tGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTATAAGGAGGCATTTTAATATTGAATAAAACGTAGGCTTAGTTCATGcc                                                                                      >  2:378854/1‑150 (MQ=14)
                                                                     ttCAAATCCAATAGATTGCATTGAAAGTGTAGTCCTGATAAGCGTAGCTCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggc                                                                             <  2:1265026/150‑1 (MQ=25)
                                                                        aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTg                                                                          <  1:1309762/150‑1 (MQ=255)
                                                                           tCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTTTAAAAGGCGACATTTTAATATTGTATACAAAGTAGGCTTTGTTCATGCCGGATGCGGCGTGaaa                                                                       >  1:1487417/1‑149 (MQ=21)
                                                                                      gCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                            >  2:543936/1‑150 (MQ=255)
                                                                                      gCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAAAGTAGGCTTTGTTCATGCCGGATGCGGCGTGATCGACTTATCCgg                                                            >  1:1591410/1‑150 (MQ=255)
                                                                                         gTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCAt                                                         <  2:1364180/150‑1 (MQ=255)
                                                                                              agtgtAGGTCTGATAAGCGTAGCGCATCAGGCAATGTTGAGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaaaa                                                    <  2:1515171/149‑1 (MQ=255)
                                                                                                tgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAcc                                                  >  1:1422335/1‑150 (MQ=255)
                                                                                                tgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAcc                                                  <  1:536156/150‑1 (MQ=255)
                                                                                                     gCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGTCGCTGTAAAATGCTTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCa                                             <  2:987796/150‑1 (MQ=255)
                                                                                                     gCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAAAGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATTCCGGATGCTGCGTGAACGCCTTATAATGCATGAAAACACTTCa                                             >  2:878890/1‑150 (MQ=18)
                                                                                                              gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTATTATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGGGTGAAAGACTTATACGGCATGAAAACCCTTCAAATCCAATa                                    >  1:618153/1‑150 (MQ=18)
                                                                                                                 tAGCTCATCAGGAAATGTTGGGTTTGTCATAAGATTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGAt                                 <  2:316658/150‑1 (MQ=14)
                                                                                                                    cttaTCATGCAATTTTGCGTTTCGAATTAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                              <  2:315507/147‑1 (MQ=14)
                                                                                                                     gcATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGAGGCGTTAAAGACTTATCCGGCATGAAAACCCTTAAAATCCAATAGaattca                             >  1:4728/1‑145 (MQ=17)
                                                                                                                            gCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATACGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACg                      >  2:1087206/1‑150 (MQ=255)
                                                                                                                                         tGTCATGAGATTAAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATa         <  2:1316985/150‑1 (MQ=255)
                                                                                                                                          gTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATaa        >  2:298239/1‑150 (MQ=255)
                                                                                                                                           tcatcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAg       >  2:930466/1‑150 (MQ=40)
                                                                                                                                           ccataaGTTTAAAATGGCGCGGTAAAAGGCGGCATTTTCAAATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGAGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGAATGCAGTGAACGTGTAGGCCTGATAAg       <  2:258182/149‑1 (MQ=11)
                                                                                                                                            catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc      >  2:353746/1‑150 (MQ=40)
                                                                                                                                            catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc      >  2:452873/1‑150 (MQ=40)
                                                                                                                                            catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCATTCAAATCCAATAGATTGCAGTGGACGTGTAGGCCTGATAAGc      >  2:1267503/1‑150 (MQ=255)
                                                                                                                                                aGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAg  >  2:1272220/1‑150 (MQ=37)
                                                                                                                                                aGTTGCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAg  <  1:371652/150‑1 (MQ=21)
                                                                                                                                                     |                                                                                                                                                
ATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAG  >  NC_000913/1552042‑1552335

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