Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,973,872 A→T 5.8% intergenic (+29/‑119) gyrA → / → serC DNA gyrase subunit A/3‑phosphoserine/phosphohydroxythreonine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,973,8720AT 5.8% 204.9 / 2.7 69intergenic (+29/‑119)gyrA/serCDNA gyrase subunit A/3‑phosphoserine/phosphohydroxythreonine aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base T (2/2);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

GCCGTCCGAGGACGAACTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTAACTTCTGCG  >  NC_002947/1973726‑1974018
                                                                                                                                                   |                                                                                                                                                  
gCCGTCCGAGGACGAACTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAgg                                                                                                                                                  <  1:214515/149‑1 (MQ=255)
 ccGTCCGAGCACGAACTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGt                                                                                                                                                 <  7:202563/149‑1 (MQ=255)
  cGTCCGAGGACGAACTCGTTGTAGTGATCGAGACGGACGAGGATGGTTTCGAGGCTGTAGCGACAGACAATGAAAATTAATATTCTGCCGGTGCAGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCTTGGCGACCAAGGTc                                                                                                                                                <  6:282206/149‑1 (MQ=255)
  cGTCCGAGGACGAACTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTc                                                                                                                                                <  8:318883/149‑1 (MQ=255)
       gAGGACGAACTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCTTGGTCAcccc                                                                                                                                           <  8:365081/149‑1 (MQ=255)
       gAGGACGAACTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGcccc                                                                                                                                           >  3:46668/1‑149 (MQ=255)
          gACGAACTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGtt                                                                                                                                        >  2:92809/1‑149 (MQ=255)
             gAACTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGa                                                                                                                                     <  4:189655/149‑1 (MQ=255)
              aaCTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGAc                                                                                                                                    <  1:8197/149‑1 (MQ=255)
                cTCGATGAAGTGATCGAGACGGACGAGGAGGATGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACtg                                                                                                                                  >  1:46607/1‑149 (MQ=255)
                   gATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACtgtg                                                                                                                                >  6:111200/1‑148 (MQ=255)
                   gATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACtgtg                                                                                                                                <  1:92809/148‑1 (MQ=255)
                            aTCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCTCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCggg                                                                                                                      >  3:266147/1‑149 (MQ=255)
                            aTCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTACACCCGAGCGGGGCGACCAAGGTAGCCCCGTTTGACTGTGTAGGAGAggg                                                                                                                      >  6:94525/1‑149 (MQ=255)
                                 gacggacgAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAACCCCGAGCGGGGCGACCACGGTCGCGCCGTTTGCCTGTGTAGGAGCGGGCTTGc                                                                                                                 >  8:200352/1‑149 (MQ=255)
                                 gacggacgAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGc                                                                                                                 <  3:212699/149‑1 (MQ=255)
                                    ggacgAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCcg                                                                                                              <  4:46668/149‑1 (MQ=255)
                                     gacgAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGct                                                                                                                        >  8:220223/1‑137 (MQ=255)
                                     gacgAGGAGGGTGTCGAGGCTGTAGAGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGGAGCCCAGTTTGACTGTGTATGAGCg                                                                                                                        >  8:214132/1‑138 (MQ=255)
                                      acgAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCcgc                                                                                                             >  5:155609/1‑148 (MQ=255)
                                          ggaggGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGtt                                                                                                                                        <  5:377086/117‑1 (MQ=255)
                                          ggaggGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGtt                                                                                                                                        >  6:377086/1‑117 (MQ=255)
                                                  tCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGAACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGcc                                                                                                 <  1:336978/149‑1 (MQ=255)
                                                    gAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCAcc                                                                                              >  3:359894/1‑149 (MQ=255)
                                                       gCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGAt                                                                                           >  7:54952/1‑149 (MQ=255)
                                                         tGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGTCTGTGTAGGAGCCGGCTTGCACGCGATGCAGGCGCCCCCGATcc                                                                                         >  3:40112/1‑149 (MQ=255)
                                                             gcgcCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCTAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGc                                                                                      <  7:280211/148‑1 (MQ=255)
                                                                ccAGACAATGAAAATGAAGAAGCTGCAGGTGCCGATGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATg                                                                                   >  6:196997/1‑148 (MQ=255)
                                                                   gACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAAt                                                                               <  5:187351/149‑1 (MQ=255)
                                                                       atgaaaatgaaGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTgg                                                                           <  3:156956/149‑1 (MQ=255)
                                                                         gaaaatgaaGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTgggc                                                                         >  7:230198/1‑147 (MQ=255)
                                                                          aaaatgaaGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCa                                                                        >  7:231704/1‑149 (MQ=255)
                                                                           aaatgaaGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCa                                                                        <  2:46607/148‑1 (MQ=255)
                                                                                  gaagCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTa                                                                >  7:21700/1‑149 (MQ=255)
                                                                                        gCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGAttt                                                          >  4:302458/1‑149 (MQ=255)
                                                                                          aGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTg                                                        <  4:170352/149‑1 (MQ=255)
                                                                                                 gacgaGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAATGTAGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCCATCCGGCATGGAATGTGGGCCCCGAATTACGATGTTGGTCTAtt                                                 >  8:407456/1‑149 (MQ=255)
                                                                                                  acgaGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTAtt                                                 <  8:374543/148‑1 (MQ=255)
                                                                                                  acgaGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGGTTGACTGTGTAGGCGCGGGCTTGCCCGCGATGCAGGCGGCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTAttt                                                >  7:374543/1‑149 (MQ=255)
                                                                                                   cgaGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTg                                               >  2:42336/1‑149 (MQ=255)
                                                                                                   cgaGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTg                                               <  4:359894/149‑1 (MQ=255)
                                                                                                       gCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCag                                           <  6:155609/149‑1 (MQ=255)
                                                                                                         cccGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCagag                                         >  2:119280/1‑149 (MQ=255)
                                                                                                         cccGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGGTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCagag                                         >  8:214072/1‑149 (MQ=255)
                                                                                                           cGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCg                                       <  1:42336/149‑1 (MQ=255)
                                                                                                            gCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCga                                      <  3:335545/149‑1 (MQ=255)
                                                                                                                     agtaagCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGGTTGACTGTGTAGGAGCGGGCTTTGCCGCGATGCAGGCGCCACCGATCCTGCATGGGATGTGGCCACCGAATTAAGATTTTGTTGTATTTGGCAGTGGGAGAGTGGAtg                             >  8:221592/1‑149 (MQ=255)
                                                                                                                      gtaagCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGAtgt                            <  4:266147/149‑1 (MQ=255)
                                                                                                                         agCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGtt                                                      >  3:336225/1‑120 (MQ=255)
                                                                                                                         agCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGtt                                                      <  4:336225/120‑1 (MQ=255)
                                                                                                                         agCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGa                          <  8:21700/148‑1 (MQ=255)
                                                                                                                          gCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAg                         >  6:395373/1‑148 (MQ=255)
                                                                                                                             tAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCaaa                     >  3:272019/1‑149 (MQ=255)
                                                                                                                               aaCCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAAc                    <  5:111200/148‑1 (MQ=255)
                                                                                                                                aCCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGACCCCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGa                  >  6:108770/1‑149 (MQ=255)
                                                                                                                                aCCCGAGCGGGGCGACCAAGGTAGCCCCGTTTGCCTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCAGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGAAGTGAGCAAAc                    >  4:329417/1‑147 (MQ=255)
                                                                                                                                  ccGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGc                >  5:422062/1‑149 (MQ=255)
                                                                                                                                  ccGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGc                >  4:47293/1‑149 (MQ=255)
                                                                                                                                  ccGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGc                <  8:54952/149‑1 (MQ=255)
                                                                                                                                    aaacGGGGCGACCTTGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCaga                                          >  4:121015/4‑121 (MQ=255)
                                                                                                                                    aaacGGGGCGACCTTGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCaga                                          <  3:121015/118‑1 (MQ=255)
                                                                                                                                      gCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCttt            <  5:395373/149‑1 (MQ=255)
                                                                                                                                         gggCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTAAc         <  4:40112/149‑1 (MQ=255)
                                                                                                                                         gggCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCCATCCTGCATCGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTAAc         >  8:115355/1‑149 (MQ=255)
                                                                                                                                          ggCGACCTTGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTAACt        >  6:379760/1‑149 (MQ=255)
                                                                                                                                           gCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTAACtt       >  8:38021/1‑149 (MQ=255)
                                                                                                                                             gACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTAACTTCt     <  4:272019/149‑1 (MQ=255)
                                                                                                                                                 cAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTAACTTCTgcg  >  2:393023/1‑149 (MQ=255)
                                                                                                                                                 cAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCAACGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTAACTTCTgag  >  2:113145/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCGTCCGAGGACGAACTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGCGACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTAACTTCTGCG  >  NC_002947/1973726‑1974018

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