Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,085 G→T 6.1% intergenic (+94/+313) fdnI → / ← yddM formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,552,0850GT 6.1% 209.2 / 3.6 65intergenic (+94/+313)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (26/35);  new base T (2/2);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

TATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTA  >  NC_000913/1551938‑1552227
                                                                                                                                                   |                                                                                                                                              
tATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCTGATGCGGCGTg                                                                                                                                                >  6:14641/1‑148 (MQ=255)
  tCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                                                                                              <  3:243257/148‑1 (MQ=255)
  tCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                                                                                              <  4:322222/148‑1 (MQ=255)
     tGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                          <  7:385487/149‑1 (MQ=255)
            gagaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                   <  8:251533/149‑1 (MQ=255)
            gagaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                   >  3:411378/1‑149 (MQ=25)
               aaGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                                                                                <  1:362459/149‑1 (MQ=255)
               aaGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                <  4:411378/149‑1 (MQ=25)
                 ggCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                 >  7:212586/1‑146 (MQ=255)
                 ggCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                 <  8:212586/146‑1 (MQ=255)
                    aGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaa                                                                                                                            >  4:85231/1‑148 (MQ=255)
                     gAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGATAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaaaa                                                                                                                          >  1:184724/1‑149 (MQ=255)
                        gCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAccc                                                                                                                       <  1:424027/149‑1 (MQ=255)
                           aaaaaaGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCtt                                                                                                                     >  3:141451/1‑148 (MQ=255)
                           aaaaaaGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCtt                                                                                                                     >  3:147611/1‑148 (MQ=255)
                             aaaaGAGAGTGAAGAAGGTATATAATCTCTTTTGAACTTTAAGCTGATAATGGCGCTGTAAATGTCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                    <  6:423924/131‑1 (MQ=255)
                             aaaaGAGAGTGAAGAAGGGATATAAACTCTTTTGACCTATAAGCTGAAAATGGCGCAGTAAAAGGCGCCATTTTCATATTGTAGACAACGAAGGATTTGTTCATGCCGGATGAGGCGTGAACGCCTTATcc                                                                                                                                    >  5:423924/1‑131 (MQ=21)
                                 gagagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATcc                                                                                                              <  7:319504/149‑1 (MQ=255)
                                   gagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCaa                                                                                                            <  3:85231/149‑1 (MQ=255)
                                            agGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCa                                                                                                   >  3:222677/1‑149 (MQ=255)
                                            agGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCa                                                                                                   <  3:223730/149‑1 (MQ=255)
                                              ggATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                    >  5:23390/1‑114 (MQ=255)
                                              ggATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                    <  6:23390/114‑1 (MQ=255)
                                              ggATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGt                                                                                                 >  4:146114/1‑149 (MQ=255)
                                                atatAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTg                                                                                                >  4:174060/1‑148 (MQ=255)
                                                    aaTCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACgt                                                                                           >  4:111296/1‑149 (MQ=255)
                                                       ctctTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACgtg                                                                                          >  1:358305/1‑147 (MQ=255)
                                                        tctTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAg                                                                                       >  2:200123/1‑149 (MQ=255)
                                                          ttttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTa                                                                                        <  7:98527/146‑1 (MQ=255)
                                                          ttttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGc                                                                                     >  7:270068/1‑149 (MQ=255)
                                                           tttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCaa                                                                                                                  >  2:209066/1‑119 (MQ=255)
                                                              gAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCa                                                                                                   >  4:198126/1‑131 (MQ=255)
                                                              gAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCa                                                                                                   <  3:198126/131‑1 (MQ=255)
                                                               aaCTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGa                                                                                 <  2:358305/148‑1 (MQ=255)
                                                               aaCTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGAt                                                                                <  3:397457/149‑1 (MQ=255)
                                                                  tttAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAg                                                                             >  1:204677/1‑149 (MQ=40)
                                                                      aGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTa                                                                         <  5:14641/149‑1 (MQ=37)
                                                                        cTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc                                                                       <  1:313148/149‑1 (MQ=35)
                                                                        cTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc                                                                       <  3:174060/149‑1 (MQ=35)
                                                                         tGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcg                                                                      >  7:371510/1‑149 (MQ=34)
                                                                          gAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCa                                                                                                   <  1:131774/119‑1 (MQ=34)
                                                                          gAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCa                                                                                                   >  2:131774/1‑119 (MQ=34)
                                                                          gAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcg                                                                      >  3:297774/1‑148 (MQ=33)
                                                                           aaaaTGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCa                                                                    <  1:200123/149‑1 (MQ=32)
                                                                               tGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGATAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                <  2:184724/149‑1 (MQ=17)
                                                                                  ctcTGTAAAAGGCGCCATTTTCATGTTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGCCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                             <  4:315752/147‑1 (MQ=12)
                                                                                           aGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGt                                                    <  3:146114/149‑1 (MQ=32)
                                                                                             gcgcCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGttt                                                  <  4:222677/149‑1 (MQ=32)
                                                                                              cgcCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTg                                                 <  4:147611/149‑1 (MQ=32)
                                                                                              cgcCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTg                                                 <  4:141451/149‑1 (MQ=32)
                                                                                                ccATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGtc                                               <  4:297774/149‑1 (MQ=11)
                                                                                                              tAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATg                                 >  4:142711/1‑149 (MQ=34)
                                                                                                              tAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATg                                 <  4:328606/149‑1 (MQ=34)
                                                                                                                 aCAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcg                              >  7:362764/1‑149 (MQ=34)
                                                                                                                          gCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAg                     <  8:362764/149‑1 (MQ=40)
                                                                                                                           cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAgg                    <  7:435122/149‑1 (MQ=40)
                                                                                                                           cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAgg                    >  1:54245/1‑149 (MQ=40)
                                                                                                                                     gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAtttt           >  4:396704/1‑148 (MQ=255)
                                                                                                                                     gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAtttt           >  3:236113/1‑148 (MQ=255)
                                                                                                                                     gCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTaaaa                      <  5:254496/137‑1 (MQ=14)
                                                                                                                                     gCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTaaaa                      >  6:254496/1‑137 (MQ=14)
                                                                                                                                      ccGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCTAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc          <  2:204677/148‑1 (MQ=255)
                                                                                                                                       cGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCat        <  3:142711/149‑1 (MQ=255)
                                                                                                                                             cggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa  <  4:236113/149‑1 (MQ=255)
                                                                                                                                                 gTGAACGCCTTGTCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa  >  2:89759/1‑145 (MQ=255)
                                                                                                                                                 gTGAACGCCTTGTCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa  <  1:89759/145‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTA  >  NC_000913/1551938‑1552227

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