Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 816,752 T→A 11.6% intergenic (‑105/‑292) ybhK ← / → moaA putative CofD superfamily transferase/molybdopterin biosynthesis protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913816,7520TA11.6% 276.3 / 4.5 86intergenic (‑105/‑292)ybhK/moaAputative CofD superfamily transferase/molybdopterin biosynthesis protein A
Reads supporting (aligned to +/- strand):  ref base T (51/24);  new base A (2/9);  total (53/33)
Fisher's exact test for biased strand distribution p-value = 2.38e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAAT  >  NC_000913/816517‑816976
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                
ttATCGGTGGTTGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGTCCCAGGCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTaaa                                                                                                                                                                                                                               <  1:1120287/239‑1 (MQ=255)
  ttCGGTGGGTGTGAAGATCCCCGTAAAACGCAACCCCAAAGACGAAAGGTAGAGAAGAACGCCCCCCCGCCCGTCCCCTCCCCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAgg                                                                                                                                                                                                                             <  2:714384/238‑1 (MQ=255)
   tCGGTGGTGGTGACGATACCCGTTAAACGAACACCCAAAGGCGAAAGTGATGAGAGAACGCGTCCCAGTCCATGCCCCCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAACCAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGa                                                                                                                                                                                                                            <  1:646442/239‑1 (MQ=255)
   gcGGTGGTGGTGACGATCCCCTTTAAACGAAACCCCAAAGAGGAAAGGGAGGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGa                                                                                                                                                                                                                            <  2:118053/238‑1 (MQ=255)
   gcGGTGGGGGTAACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGGCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAAATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGa                                                                                                                                                                                                                            <  1:452799/238‑1 (MQ=255)
       tggtggtgACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTTCCAATATTGAGGATACAAACTATTTTATTCTACCTCAAAAAGACGaa                                                                                                                                                                                                                        >  2:1305099/1‑238 (MQ=255)
           ggtgACGATCCCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCa                                                                                                                                                                                                                     <  1:264690/238‑1 (MQ=255)
            gtgACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCAc                                                                                                                                                                                                                    >  1:118053/1‑238 (MQ=255)
               aCGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATAAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGCGGATACAAACTATTTTATTTTACCTCAAAAGAAAGAAGAACGcag                                                                                                                                                                                                                >  2:716824/1‑237 (MQ=255)
                cGATACCCGTTAAACGAGAACCCAAAGCCGAAAGGGATGAGAGAACGCGTCCCAGCCCAGGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTa                                                                                                                                                                                                                <  1:917270/238‑1 (MQ=255)
                 gATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGTTACAAACCATTTTTTTCTACCTTTAAAGGACAATgaaaccaat                                                                                                                                                                                                               >  2:1202001/1‑230 (MQ=255)
                 gATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGGACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCAACCTCTAAAGGACGATGCAAGCTAt                                                                                                                                                                                                               >  2:134557/1‑238 (MQ=255)
                       ccttAAACGGGACCCCAAAGACGAAAGTGAGAGGGAAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGACTccc                                                                                                                                                                                                        <  2:1137345/237‑1 (MQ=255)
                         ttAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCtg                                                                                                                                                                                                      >  1:305465/1‑239 (MQ=255)
                         ttAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGAAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGAAACAAACAATTTTCTTCTACCTCTAAAGAACGATGAACGCTATGCCTCCCtg                                                                                                                                                                                                      >  2:1273357/1‑239 (MQ=255)
                           aaaCGAGAACCCAAAGACGAAAGTGATGAGCGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCATCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGGGGATACAAACTATTTTTTTCTACCTCTAAAggc                                                                                                                                                                                                                            >  1:1352414/1‑214 (MQ=255)
                            aaCGGGCACCCAAAGGCGAAAGTGATGAGAAACGCCGTCCCAGCCCAGGCCCCGCGCCGAGAGAAACGACCCGAACAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATGATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCtgatg                                                                                                                                                                                                   <  2:290560/239‑1 (MQ=255)
                            aaCCAGAACCCAAAGACGAAAGGGATGAGAACACCGCTCCCAGTCCAGGCCCTCCGCCGGAGCCAACGAAACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCtgatg                                                                                                                                                                                                   <  2:58926/239‑1 (MQ=255)
                              agagaACCACAAGCAGAACGTGATGAGAAAACGCGCCCCAGTCAATGCACTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTgatgat                                                                                                                                                                                                 <  1:1012652/238‑1 (MQ=255)
                                  aaCCCAAAGACGAAAGTGATGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGACCGCTATGCCTCCCTGATGATGtat                                                                                                                                                                                             >  1:335042/1‑239 (MQ=255)
                                       aaaGACGAAAGTGATGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCAGCTGATGTATATCaa                                                                                                                                                                                        >  1:807001/1‑239 (MQ=255)
                                               aaGTGATGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTaaaa                                                                                                                                                                                 >  1:731026/1‑238 (MQ=255)
                                                              aCGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATTAATATCAAAGAAAAACCGCCATATTCCCt                                                                                                                                                                   >  2:18570/1‑239 (MQ=255)
                                                                   gtccaacgcaaTGCCCTCCGCCGAGACAACGCACCGGTCAAAGATCATCCAGCAAAGGATGCCTCATATAGTTCCCGGGAGTCAGATTACCCGCTATACAGGAGCACAATTTGGGGGAAAAAGAAATTAACCCTGCAAATATGGAGGAAACAAATTATTTTCTTCTACCTCAAAGGGACGAGGCAGGCAATCCCTCCGTGATGGTGTATATAAAAGTAAAACCGCCATTTTCCCTTAtt                                                                                                                                                               <  1:434395/229‑1 (MQ=255)
                                                                    tCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCATTATTGCGGATCCAAACTATTTTCTCCTACATCTAAAGGACGATGCACGCTATGCACACCTAATGATGTATATCAAATTAAAACCCCCCATTTCCCTTAttt                                                                                                                                                              >  2:434395/1‑238 (MQ=255)
                                                                      ccAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCGCATGATGTAAATCAAAGTAAAACACCCATTTTCCCTTAATCTg                                                                                                                                                            >  2:673936/1‑239 (MQ=255)
                                                                           ccATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCATTTTCCCTTATTCTGTAg                                                                                                                                                         <  1:1305099/236‑1 (MQ=255)
                                                                               gCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATg                                                                                                                                                                                                              >  1:378921/1‑180 (MQ=255)
                                                                               gCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATg                                                                                                                                                                                                              <  2:378921/180‑1 (MQ=255)
                                                                               gCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCCTTTTCCCTTTTTTTTTAGCGaaaa                                                                                                                                                   >  2:391070/1‑238 (MQ=255)
                                                                                 ccaccgccgAGAGAAACGTCAGGATAAATACCAGCCAGCGTCGGTTTGCCCATATATTTCCCTGGAGCCAGATTACCCGCGCTACAGAAGCGCAATTTGCGGGAAAACAGCAATAAACCTCCAAATATGGAGGATACAACTTATTTTCTTCTACCTCAAAAGGACGAGGCACGTTTTGCCTCCCTGATGATGTATATAAAAGAAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAg                                                                                                                                                 <  1:18570/236‑1 (MQ=255)
                                                                                    ccgccgAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCACTATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGAACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc                                                                                                                                              >  1:61445/1‑239 (MQ=255)
                                                                                    ccgccgAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc                                                                                                                                              >  1:157806/1‑239 (MQ=255)
                                                                                    ccgccgAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc                                                                                                                                              >  1:90305/1‑239 (MQ=255)
                                                                                    ccgccgAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc                                                                                                                                              >  1:1070663/1‑239 (MQ=255)
                                                                                    ccgccgAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc                                                                                                                                              >  1:1101244/1‑239 (MQ=255)
                                                                                    ccgccgAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc                                                                                                                                              >  1:422662/1‑239 (MQ=255)
                                                                                    ccgccgAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTAATTCTTTAGCGTAATAGCAc                                                                                                                                              >  1:372002/1‑239 (MQ=255)
                                                                                            gaaaaaCCACACGAAAAAGATCAGCCAGCGTACGATTGCTCATAGAGTTGCCGGGAGTCAAATTACCCGCGCAACGGTAGCGAAAATTGGGGGAAAACAGAAATTACCCTGCAAATATTGAGGATACAAACTTTTTTCTTCTACCTCTAAAGGACGATGCACGTTATGCCTCCCTGATGATGTATATCAAAGAAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCAt                                                                                                                                        <  1:673936/234‑1 (MQ=255)
                                                                                              caaaCGACAAGCACAAGCTCACCAAGCGTACTATTGCGCATAGTGTTTACGGGAGTCAGATTACCCGCGACACAGTAGCCCAAATTGGGGGAAAAAAGAAATTACCCTGCCAATATTGAGGATACAAACTTTTTTCTTCTACCTCTAAAGGACGTTGCACGCTAGGCCTCCCTGAGGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGa                                                                                                                                      <  1:716824/236‑1 (MQ=255)
                                                                                               cAAAGACACGAAAAAGATCGGCCAGCGAACAATTGCCAAAATATTTTCCGGGAGCAAGATTACCCGGGCTAAAGTAGCGAAATTTGGGGGAAAACAGAAATAAACCTGCCAATATTGAGGTTACAAACTTTTTTTTTCTCCCTCAAAAGGACGATGCACGTTATGCCTCCCTGATGATGTATATAAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACg                                                                                                                                    <  1:1202001/239‑1 (MQ=255)
                                                                                                  cGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCta                                                                                                                                 >  1:762644/1‑239 (MQ=255)
                                                                                                  cGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACAATCATGACCCta                                                                                                                                 >  2:802798/1‑239 (MQ=255)
                                                                                                       cGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGAAAAACCGCCATTTTCCCTTATTCTGTAGCGAACAAGCCACATCATGACACTATATAc                                                                                                                            >  2:86458/1‑239 (MQ=255)
                                                                                                                   agccagcGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATa                                                                                                                 >  1:1221095/1‑239 (MQ=255)
                                                                                                                     ccagcGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGGAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAAAAGCACAATCATGACCCAATAAAAATGATAACAAAgg                                                                                                               >  2:552114/1‑238 (MQ=255)
                                                                                                                     ccagcGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAAAAGAACGACCATGAAGCAATATACAAAATTaaaaaga                                                                                                               >  2:487154/1‑233 (MQ=255)
                                                                                                                         cGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGcaa                                                                                                           >  1:354431/1‑239 (MQ=255)
                                                                                                                                             ttCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACg                                                                                                                                                                                                                   <  1:1291433/115‑1 (MQ=255)
                                                                                                                                             ttCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACg                                                                                                                                                                                                                   >  2:1291433/1‑115 (MQ=255)
                                                                                                                                              tCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCa                                                                                      >  1:1040948/1‑239 (MQ=255)
                                                                                                                                               ccgggatccAGTTTATCCGGGCTAAAGTAGCAAAAATTGCGGGGAAACAGAAATTAACCGGCAAATTGTGGGGATAAAACCTATTTCTTTCTACCTCAAAAGGGCGGTGCAAGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCAt                                                                                     <  1:309411/231‑1 (MQ=255)
                                                                                                                                                  ggttgcattttACCCGCGTTACAGTGGCAAAATTTGGGGGAAAACAGAAATTACCCTGCCAATATTGATGATAAAAATTTTTTTCTTCTCCCTCTAAAGGACGATGCACGTTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTa                                                                                  <  2:1040948/230‑1 (MQ=255)
                                                                                                                                                   gAGTCAGATTACCCGCGCTACAGTAGCGAAATTTGGGGGAAAACAGAATTGACCCTGCAAATATGGAGGATACAAACTTTTTTTTTCTACCTCAAAAGGACGATGCACGCTATGCCTCCCGGATGATGTATATCAAAGAAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTa                                                                                  <  1:1140507/238‑1 (MQ=255)
                                                                                                                                                      tCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGAACATGACGCTATTTACATGATTAAATAGCGAAAGTGGGGAAGGGAAAAATACCATTTaaac                                                                               >  2:449659/1‑238 (MQ=255)
                                                                                                                                                       gaatttccccGCTCTAGAGTACGAAAATTGGGGGAAAAACAAAAATTACCCGGCAAATTTGGGGGAACAAAACTTTTTTCTTCTACCACAAAAGGAAGCTCCCCTCTATCCCCCCCTGATGATGTAAAAAAAAATAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGc                                                                             <  2:157806/231‑1 (MQ=255)
                                                                                                                                                        aGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGCAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGGGTGGATGGTAAAAATCTCATTTACACGCt                                                                            >  1:1112611/1‑239 (MQ=255)
                                                                                                                                                        aGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTCCATAGCGAAAGGGTGGATGGGAAAAATCTCATTTACACGCt                                                                            >  1:305793/1‑239 (MQ=255)
                                                                                                                                                          aTTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTa                                                                                                                     >  1:696649/1‑196 (MQ=255)
                                                                                                                                                                      aCAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACAt                                                                                                                          >  1:709427/1‑179 (MQ=255)
                                                                                                                                                                        aGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAACTCTCATTTACACGCTAGTATCGGCATAACCa                                                            >  1:661419/1‑239 (MQ=255)
                                                                                                                                                                           agccaaaaTTGCGGAAAAACAAAATTTAACCCACCAAGTTGGGGGATAAAATCTATTTTCTCCAACCTCAAAAGGCAGACGCCCGTGATCCCTCCCTTATGATGTATAAAAAATTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACt                                                          <  1:86458/233‑1 (MQ=255)
                                                                                                                                                                             cgcAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAGAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATCGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGCTGGTAAAAATCTCCTTTACACGCTAGAATCGGCATAACCACTaaa                                                       >  1:1291944/1‑239 (MQ=255)
                                                                                                                                                                             cgcAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTcc                                                                                                                                                                                                         <  1:160945/93‑1 (MQ=255)
                                                                                                                                                                             cgcAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTcc                                                                                                                                                                                                         >  2:160945/1‑93 (MQ=255)
                                                                                                                                                                             cgcAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTatat                                                                                                                                                                                           >  2:484619/1‑107 (MQ=255)
                                                                                                                                                                             cgcAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTatat                                                                                                                                                                                           <  1:484619/107‑1 (MQ=255)
                                                                                                                                                                             cgcAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCtatata                                                                                                                             >  1:1193559/1‑169 (MQ=255)
                                                                                                                                                                             cgcAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACAt                                                                                                                          >  1:1191696/1‑172 (MQ=255)
                                                                                                                                                                             cgcAAAGTGGGGGAAAACGACAATTAACCGCCAAATTTAGGAGATCAAAATTATTTTCTTCTACCTTAAAAGGACGATGCACGTTATGCCTCCTTGATGAGGTATAACAAAGAAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTaa                                                        <  1:134557/238‑1 (MQ=255)
                                                                                                                                                                               cAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGAACGATCATGACGCTATATACATTATTACATAGCGAACGCGTGGCTTGTAAAACTCTCGTTTACACGAAAGAATCGGCAGAACCACAAAaca                                                     >  2:1025781/1‑239 (MQ=255)
                                                                                                                                                                                      cGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGGGTGGATGGTAAAAATCTCATTTACACGCTAGAATCGGCATAACCACTAAAAACCCTaa                                               >  2:1071928/1‑237 (MQ=255)
                                                                                                                                                                                           aaaaaAGAAATTACCCTGCCAATATGGAGGATACAAACTATTTCCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCctc                                           <  1:1071928/237‑1 (MQ=255)
                                                                                                                                                                                               cAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCtga                                                                                                                                                                                                     >  2:935118/1‑80 (MQ=255)
                                                                                                                                                                                               cAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCtga                                                                                                                                                                                                     <  1:935118/80‑1 (MQ=255)
                                                                                                                                                                                                    tttACCCGGCAAATATTGAGAAAAAAAACTTTTTCCTTCTCCCCCTAAGGGACGATGCCCGTTATGCCTCCTTGATGATGTATACCAAAGTAAATCCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTggg                                 <  1:578142/238‑1 (MQ=255)
                                                                                                                                                                                                             gCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACAt                                                                                                                          >  1:1301237/1‑142 (MQ=255)
                                                                                                                                                                                                             gCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCCTAACCAATAAAAACCCTATCCTCTCAACCCGGGGCAACGGa                         >  2:160198/1‑239 (MQ=255)
                                                                                                                                                                                                               cAATATTGGGGATAAAACTTATTTTCTTCTACCTCAAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATAAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATa                       <  1:449659/239‑1 (MQ=255)
                                                                                                                                                                                                                        ggATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCa                                                                                                                                                                                                                     >  2:560289/1‑40 (MQ=255)
                                                                                                                                                                                                                        ggATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCa                                                                                                                                                                                                                     <  1:560289/40‑1 (MQ=255)
                                                                                                                                                                                                                            aCAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCtga                                                                                                                                                                                                     >  2:1035672/1‑52 (MQ=255)
                                                                                                                                                                                                                            aCAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCtga                                                                                                                                                                                                     <  1:1035672/52‑1 (MQ=255)
                                                                                                                                                                                                                            aCAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGc                                              >  1:846077/1‑203 (MQ=255)
                                                                                                                                                                                                                            aCAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTTGcc          >  1:659264/1‑239 (MQ=255)
                                                                                                                                                                                                                                     tttCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAAt  >  1:105222/1‑238 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                
TTATCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAAT  >  NC_000913/816517‑816976

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