Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,484,358 G→A 6.2% C275Y (TGC→TAC)  purH → bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,484,3580GA 6.2% 190.9 / 2.7 64C275Y (TGC→TAC) purHbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base G (29/31);  new base A (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

GTTATGGCGAGAACCCGCACCAGAGCGCGGCGTTCTATGTAGAAGCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCA  >  NC_002947/5484217‑5484507
                                                                                                                                             |                                                                                                                                                     
gTTATGGCGAGAACCCGCAACAACGCGCCGCGCTCTATGTAGAATCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCGTTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAaga                                                                                                                                               <  2:1250395/150‑1 (MQ=255)
    tGGCGAGAACCCGCACCAGAGCGCGGCGTTCTATGTAGAAGCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCtt                                                                                                                                           >  1:163430/1‑150 (MQ=255)
     ggCGAGAACCCGCACCAGAGCGCGGCGTTCTATGTAGAAGCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTc                                                                                                                                          <  2:576376/150‑1 (MQ=255)
            aCCCGCACCAGAGCGCGGCGTTCTATGTAGAAGCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGc                                                                                                                                   <  1:159187/150‑1 (MQ=255)
              ccGCACCAGAGCGCGGCGTTCTATGTAGAAGCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCg                                                                                                                                 >  1:553755/1‑150 (MQ=255)
              ccGCACCAGAGCGCGGCGTTCTATGTAGAAGCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGACGATACTGACGCGGCGCTGGAATGGGTAAAGAGCTTCGTCAAGCCg                                                                                                                                 >  2:1249615/1‑150 (MQ=255)
                gCACCAGAGCGCGGCGTTCTATGTAGAAGCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAATGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGc                                                                                                                               <  1:1605544/150‑1 (MQ=255)
                    cAGAGCGCGGCGTTCTATGTAGAAGCAAAGAAAGGCGCAGCCAGCATTACCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGc                                                                                                                           >  2:185816/1‑150 (MQ=255)
                     agagCGCGGCGTTCTATGTAGAAGCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCg                                                                                                                          >  1:195604/1‑150 (MQ=255)
                            ggcgTTCTATGTAGAAGCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGt                                                                                                                   >  2:714425/1‑150 (MQ=255)
                                    aTGTAGAAGCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCAcc                                                                                                           >  2:23468/1‑149 (MQ=255)
                                      gTAGAAGCCAAGAAAGGCGAAGCCAGCATTTACACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTc                                                                                                                        >  1:1340257/1‑135 (MQ=255)
                                      gTAGAAGCCAAGAAAGGCGAAGCCAGCATTTACACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTc                                                                                                                        <  2:1340257/135‑1 (MQ=255)
                                         gAAGCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGcc                                                                                                                                  >  2:1548077/1‑122 (MQ=255)
                                         gAAGCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGcc                                                                                                                                  <  1:1548077/122‑1 (MQ=255)
                                         gAAGCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAAc                                                                                                      >  2:1243143/1‑150 (MQ=255)
                                            gCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGc                                                                                                                               >  2:922865/1‑122 (MQ=255)
                                            gCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGc                                                                                                                               <  1:922865/122‑1 (MQ=255)
                                            gCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATACGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATc                                                                                                                     >  1:1530265/1‑132 (MQ=255)
                                            gCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATACGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATc                                                                                                                     >  1:1529847/1‑132 (MQ=255)
                                            gCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATACGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATc                                                                                                                     <  2:1530265/132‑1 (MQ=255)
                                            gCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATACGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATc                                                                                                                     <  2:1529847/132‑1 (MQ=255)
                                              cAAGGAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTg                                                                                                 <  1:22618/150‑1 (MQ=255)
                                                     ggCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTg                                                                                          <  1:1243143/150‑1 (MQ=255)
                                                      gCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTgg                                                                                         <  1:236746/150‑1 (MQ=255)
                                                       cGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGc                                                                                        <  2:938910/150‑1 (MQ=255)
                                                           gCCAGCATTTCCACCGCGGTCCAGCTGCAAGGCAAGGAACTGTTATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCTTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCagta                                                                                    <  2:69269/150‑1 (MQ=255)
                                                              aGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAgtagta                                                                                 >  1:123364/1‑150 (MQ=255)
                                                              aGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAgtagta                                                                                 >  1:328497/1‑150 (MQ=255)
                                                                     ccACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAg                                                                          >  2:1394212/1‑150 (MQ=255)
                                                                         ctcTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCt                                                                                                                             <  2:1413038/93‑1 (MQ=255)
                                                                         cGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCt                                                                                                                             >  1:1413038/1‑95 (MQ=255)
                                                                         cGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGa                                                                      >  2:1294087/1‑150 (MQ=255)
                                                                            tGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAgg                                                                   <  1:185816/150‑1 (MQ=255)
                                                                             gTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAggc                                                                  <  2:163430/150‑1 (MQ=255)
                                                                                 aGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCa                                                              <  2:195604/150‑1 (MQ=255)
                                                                                    tGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATcc                                                           >  2:1286524/1‑150 (MQ=255)
                                                                                          gcaaggAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAg                                                     <  1:714425/150‑1 (MQ=255)
                                                                                                  cTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGAc                                             >  1:1591098/1‑150 (MQ=255)
                                                                                                   tGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGTCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTa                                                 >  2:1577321/1‑145 (MQ=255)
                                                                                                   tGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTa                                                 <  1:1577321/145‑1 (MQ=255)
                                                                                                     tCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTg                                          <  1:1294087/150‑1 (MQ=255)
                                                                                                       aTTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGc                                        <  1:1309126/150‑1 (MQ=255)
                                                                                                       aTTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGc                                        >  1:307920/1‑150 (MQ=255)
                                                                                                        ttCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTg                                          >  2:97524/1‑147 (MQ=255)
                                                                                                        ttCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTg                                          <  1:97524/147‑1 (MQ=255)
                                                                                                        ttCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGcc                                       >  2:748437/1‑150 (MQ=255)
                                                                                                            acaacGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAg                                                                    >  1:58412/1‑117 (MQ=255)
                                                                                                            acaacGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAg                                                                    <  2:58412/117‑1 (MQ=255)
                                                                                                            acaacGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGcc                                       <  1:901596/146‑1 (MQ=255)
                                                                                                            acaacGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGcc                                       >  2:901596/1‑146 (MQ=255)
                                                                                                             caacGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAggcggc                                                               <  2:1223112/121‑1 (MQ=255)
                                                                                                             caacGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAggcggc                                                               >  1:1223112/1‑121 (MQ=255)
                                                                                                              aacGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTAcgcg                                 >  2:409030/1‑150 (MQ=255)
                                                                                                               acGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGcga                                >  2:1444072/1‑150 (MQ=255)
                                                                                                                cGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGcgac                               <  2:328497/150‑1 (MQ=255)
                                                                                                                   ggCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGcgaccga                            <  1:748437/150‑1 (MQ=255)
                                                                                                                   ggCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGcgaccga                            <  1:1286524/150‑1 (MQ=255)
                                                                                                                          aCTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGaa                     <  2:123364/150‑1 (MQ=255)
                                                                                                                           cTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAAt                    <  1:409030/150‑1 (MQ=255)
                                                                                                                              aCGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCgg                 >  1:566524/1‑150 (MQ=255)
                                                                                                                                gccgcGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAggcg                                                                 <  1:227654/100‑1 (MQ=255)
                                                                                                                                gccgcGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAggcg                                                                 >  2:227654/1‑100 (MQ=255)
                                                                                                                                   gcgcTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTc            <  2:553755/150‑1 (MQ=255)
                                                                                                                                             gCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTAAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGcatca  <  1:1249615/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
GTTATGGCGAGAACCCGCACCAGAGCGCGGCGTTCTATGTAGAAGCCAAGAAAGGCGAAGCCAGCATTTCCACCGCTGTCCAGCTGCAAGGCAAGGAACTGTCATTCAACAACGTGGCCGATACTGACGCCGCGCTGGAATGCGTAAAGAGCTTCGTCAAGCCGGCCTGCGTCATCGTCAAGCACGCCAACCCGTGCGGCGTGGCAGTAGTACCGGAAGACGAAGGCGGCATCCGCAAAGCCTACGACCTGGCCTACGCGACCGACACCGAATCGGCGTTCGGCGGCATCA  >  NC_002947/5484217‑5484507

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