Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,852,543 A→G 5.1% P136P (CCA→CCG hypD → hydrogenase maturation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,852,5430AG 5.1% 248.7 / 2.9 78P136P (CCA→CCGhypDhydrogenase maturation protein
Reads supporting (aligned to +/- strand):  ref base A (40/34);  new base G (2/2);  total (42/36)
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.80e-01

GGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCCAGCACATTACGCTTATCCCGACGTTGCGCAGTTTG  >  NC_000913/2852402‑2852690
                                                                                                                                             |                                                                                                                                                   
ggAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCg                                                                                                                                                >  4:375102/1‑147 (MQ=255)
  aaGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCTGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGATAATCCAACCCGCaa                                                                                                                                             <  5:457532/148‑1 (MQ=255)
      cATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAgtg                                                                                                                                         >  8:441276/1‑148 (MQ=255)
        taatatGTACCTTTGGCGAAGCCATGCCCCAGCCGGGGAAACAGGGATCGCAGTTGCAGGCAAAAGCACGCCGAGCCGATGTACCCATCGATTAATCTACGATGGAAGCGTTCAAACTTGAGCAGGAGAATCCAACACGCAAAgtgct                                                                                                                                       >  8:367652/6‑146 (MQ=255)
           tctGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGgtgt                                                                                                                                     >  6:203837/1‑147 (MQ=255)
             tGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGttct                                                                                                                                  <  1:179502/148‑1 (MQ=255)
                     tGGCGACGCCATGCGCGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCtt                                                                                                                           <  3:142860/147‑1 (MQ=255)
                        cGACGCCATGCGCGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGt                                                                                                                       >  8:457748/1‑148 (MQ=255)
                         gACGCCATGCGCGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGtt                                                                                                                      <  2:160649/148‑1 (MQ=255)
                          aCGCCATGCGCGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGttt                                                                                                                     <  4:171757/148‑1 (MQ=255)
                            gCCATGCGCGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTg                                                                                                                   <  7:141402/148‑1 (MQ=255)
                              cATGCGCGTGACTGGGAAACTGGGATCGCTGTTGCAGTCAAATGCACGCGGTGCCGATGTGCGCATCGTTTACTCGTCGATGGATGCGTTGAATCTGGTGCAGGAGAATCCAACGCGCAAATTGGTGTTCTTCGTCTTAGGTTTTGaa                                                                                                                 <  4:438180/148‑1 (MQ=255)
                               aTGCGCGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGtt                                                                                                                      <  4:188670/142‑1 (MQ=255)
                               aTGCGCGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGtt                                                                                                                      >  3:188670/1‑142 (MQ=255)
                                 gcgcGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGaa                                                                                                                 >  2:209426/1‑145 (MQ=255)
                                 gcgcGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGaa                                                                                                                 <  1:209426/145‑1 (MQ=255)
                                 gcgcGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAAc                                                                                                               >  5:413513/1‑147 (MQ=255)
                                 gcgcGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAAc                                                                                                               >  8:324116/1‑147 (MQ=255)
                                 gcgcGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGAGGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGAAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGATAc                                                                                                               >  5:404443/1‑147 (MQ=255)
                                 gcgcGTGCCGGGGAAAAAGGGATCGCTGTTGCAGGCAAAAGGACGCGGTGCCGATGTGCGCATCGTGTCCTCGCCGATGGATGCGGTGAAACTGGCGCACGAGAATCCAACCCGCAAAGTCGTGTTCTTCGGCTTAGGTTATGAAAc                                                                                                               >  6:104992/1‑147 (MQ=255)
                                    cgTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGtt                                                                                                                                    >  5:184008/1‑123 (MQ=255)
                                    cgTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGtt                                                                                                                                    <  6:184008/123‑1 (MQ=255)
                                    cgTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACt                                                                                                           <  3:189495/148‑1 (MQ=255)
                                    cgTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACt                                                                                                           <  7:441276/148‑1 (MQ=255)
                                    cgTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACt                                                                                                           <  3:375102/148‑1 (MQ=255)
                                    cgTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACt                                                                                                           <  1:186194/148‑1 (MQ=255)
                                    cgTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACt                                                                                                           <  5:203837/148‑1 (MQ=255)
                                         cGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGcc                                                                                                      >  5:154412/1‑148 (MQ=255)
                                          ggggAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACt                                                                                                           >  2:325054/1‑142 (MQ=255)
                                          ggggAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACt                                                                                                           <  1:325054/142‑1 (MQ=255)
                                          ggggAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGcc                                                                                                      >  3:140605/1‑147 (MQ=255)
                                               aaCAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCgg                                                                                                                              <  1:226759/118‑1 (MQ=255)
                                               aaCAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCgg                                                                                                                              >  2:226759/1‑118 (MQ=255)
                                                aCAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGaccacc                                                                                               <  2:190909/148‑1 (MQ=255)
                                                   gggATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGTTGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCAAACCCGCAAAGTGGTGTTCTTCGGCTTAGGttt                                                                                                                     <  7:35417/123‑1 (MQ=255)
                                                    ggATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCAAACCCGCAAAGTCCTGTTCTTCGGCTTAGGtatt                                                                                                                    >  8:35417/1‑123 (MQ=255)
                                                     gATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTAt                                                                                          >  6:114591/1‑148 (MQ=255)
                                                          cTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCGACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCCctc                                                                                     >  7:24956/1‑148 (MQ=255)
                                                          cTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCGACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCActc                                                                                     >  7:32207/1‑148 (MQ=255)
                                                            gTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCAGATGGATGCGTTGAAACTGGCGAAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTg                                                                                   >  3:311693/1‑148 (MQ=255)
                                                             ttGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTg                                                                                   <  6:381154/147‑1 (MQ=255)
                                                                  ggCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTAGTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTTCAACAg                                                                             <  2:162237/148‑1 (MQ=255)
                                                                  ggCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCGACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCCCTCTGCAACAg                                                                             <  8:24956/148‑1 (MQ=255)
                                                                  ggCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCGACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAg                                                                             <  8:32207/148‑1 (MQ=255)
                                                                      aaaGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGaccac                                                                                                <  3:192673/125‑1 (MQ=255)
                                                                       aaGAACGAGCTGCCGATGTGCGCATCGTTTACTAGAATATGGATGCGTTGAAACTGGCGCAGGAGAATCAAACACGCAAAGTGGAGTTCTTCGGCTTATGTATTGAAACCACTATGACGacaacc                                                                                               >  4:192673/1‑125 (MQ=255)
                                                                        aGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGa                                                                         >  3:2354/1‑146 (MQ=255)
                                                                            cgcgGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAgcgc                                                                   >  5:368640/1‑148 (MQ=255)
                                                                                          gcgcATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAAt                                                     >  2:231540/1‑148 (MQ=255)
                                                                                                cGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGa                                                               <  7:392993/132‑1 (MQ=255)
                                                                                                cGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGa                                                               >  8:392993/1‑132 (MQ=255)
                                                                                                  tttACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAAttttt                                                 <  7:253724/144‑1 (MQ=255)
                                                                                                   ttACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTActtc                                            >  2:400890/1‑148 (MQ=255)
                                                                                                     aCTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTActtct                                           <  2:225136/147‑1 (MQ=255)
                                                                                                      cTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACttcttc                                         >  7:71821/1‑148 (MQ=255)
                                                                                                         gCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGaa                                                                        >  1:175886/1‑114 (MQ=255)
                                                                                                         gCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGaa                                                                        <  2:175886/114‑1 (MQ=255)
                                                                                                         gCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTg                                       <  5:114591/147‑1 (MQ=255)
                                                                                                         gCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCAGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCAATCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGc                                      >  3:428810/1‑148 (MQ=255)
                                                                                                              tggatgCGTTTAAACTGGCGCAGGAGAATACAACCAGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCCCTCTGCAACAGGCGAAAGCGCGTGATGAGAAGAATTTTTACTTCTTCTGCCAgaa                                 >  8:350461/1‑146 (MQ=255)
                                                                                                              tggatgCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCCAGca                                 >  2:246539/1‑148 (MQ=255)
                                                                                                                gatgCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGACACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCAAGcaca                               >  6:31706/1‑148 (MQ=255)
                                                                                                                gatgCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCCAGcac                                <  4:140605/147‑1 (MQ=255)
                                                                                                                  tgCGTTGAAACTGGCGCAGGAGAATCCAACCCGCaa                                                                                                                                             <  3:181206/36‑1 (MQ=255)
                                                                                                                  tgCGTTGAAACTGGCGCAGGAGAATCCAACCCGCaa                                                                                                                                             >  4:181206/1‑36 (MQ=255)
                                                                                                                  tgCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCCAGCACAt                              >  2:380409/1‑147 (MQ=255)
                                                                                                                       tGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCCAGCACATTAc                           <  1:311284/145‑1 (MQ=255)
                                                                                                                       tGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCCAGCACATTAc                           >  2:311284/1‑145 (MQ=255)
                                                                                                                       tGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCCAGCACATTACGCt                        >  6:311498/1‑148 (MQ=255)
                                                                                                                                cgcAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCCAGCACATTACGCTTATCCCGAc               <  7:367652/148‑1 (MQ=255)
                                                                                                                                  cAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCCAGCACATTACGCTTATCCCGACGt             >  6:245188/1‑148 (MQ=255)
                                                                                                                                    ggAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCCAGCACAt                              >  6:404235/1‑129 (MQ=255)
                                                                                                                                    ggAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCCAGCACAt                              <  5:404235/129‑1 (MQ=255)
                                                                                                                                       gaATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCCAGCACATTACGCTTATCCCGACGTTgcgc        <  5:311498/148‑1 (MQ=255)
                                                                                                                                        aaTCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCCAGCACATTACGCTTATCCCGACGTTgcg         >  4:183217/1‑146 (MQ=255)
                                                                                                                                           ccAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAAttttt                                                 >  7:404928/1‑103 (MQ=255)
                                                                                                                                           ccAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAAttttt                                                 <  8:404928/103‑1 (MQ=255)
                                                                                                                                           ccAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCCAGCACATTACGCTTATCCCGACGTTGCGCAGtt    >  4:28009/1‑148 (MQ=255)
                                                                                                                                             aaCCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCCAGCACATTACGCTTATCCCGACGTTGCGCAGTTTg  <  6:154412/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACAGGGATCGCTGTTGCAGGCAAAAGCACGCGGTGCCGATGTGCGCATCGTTTACTCGCCGATGGATGCGTTGAAACTGGCGCAGGAGAATCCAACCCGCAAAGTGGTGTTCTTCGGCTTAGGTTTTGAAACCACTATGCCGACCACCGCTATCACTCTGCAACAGGCGAAAGCGCGTGATGTGCAGAATTTTTACTTCTTCTGCCAGCACATTACGCTTATCCCGACGTTGCGCAGTTTG  >  NC_000913/2852402‑2852690

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: