Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,487,141 C→T 14.8% intergenic (+77/‑70) purH → / → purD bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,487,1410CT14.8% 421.2 / 60.9 185intergenic (+77/‑70)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base C (79/76);  new base T (16/11);  total (95/90)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-02

Only 100 of 201 total aligned reads displayed.

TGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAA  >  NC_002947/5486992‑5487290
                                                                                                                                                       |                                                                                                                                                     
      cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                                 <  2:3316131/150‑1 (MQ=37)
      cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                                 <  2:6860532/150‑1 (MQ=37)
      cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                                                 <  1:3832214/150‑1 (MQ=255)
       aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCTTGCGCCGAACACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                                 <  1:44091/150‑1 (MQ=14)
       aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgcg                                                                                                                                                <  1:5745670/150‑1 (MQ=37)
       aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                                                <  1:25157/150‑1 (MQ=255)
       aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                                                >  2:8364020/1‑150 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGTCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGGGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCTTTCGCgg                                                                                                                                               <  1:1968724/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                               <  1:149632/150‑1 (MQ=37)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                               <  1:3802355/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                               <  2:5322921/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                               <  1:3809658/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                               <  2:7540239/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                               <  2:8340206/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                               <  1:7049690/150‑1 (MQ=255)
          gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGCCCGCGCAAAACATGGGGGATTCACCGCCGCCGTCGCGGGc                                                                                                                                             >  2:8688192/1‑150 (MQ=255)
            tgctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCaa                                                                                                                                           <  2:4499665/150‑1 (MQ=255)
                   gacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTTCCCCCGCAATAGATGGGGGAATGACACCCGCATTAGCGGGCAAGCCCGcg                                                                                                                                    >  2:8589231/1‑149 (MQ=255)
                   gacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGGCCCCGCATAGGGTGGTGGTGTGACCGCCCCCGGAGCGGGCCAGcccccg                                                                                                                                    >  1:686031/1‑147 (MQ=255)
                    acgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCCCGGAAACCAAGGTGGGGTCACCGCCCCAATAGCGGgcccgcccccgc                                                                                                                                   >  2:4790203/1‑146 (MQ=255)
                                tggcgTTGGTATTCACTGGCAGGCGCCACTTCCCCAAATATTTACGCGTATCCCGAGCAAAGTCTAGATCCGGTGGGAGCGGGCTTGCCCGCTAATACTATGGTGGATACCCCGCCGCATTCGCGGGAAAGCCCGCTCCCACAGTGTTCg                                                                                                                       <  2:5465089/148‑1 (MQ=255)
                                            tCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACg                                                                                                           <  2:5944258/150‑1 (MQ=255)
                                             cACTGGCTGGCGCCACTCCCCCCACTATTTCCGCGGATCCGGTGCAATGTCGAGATCCTGTGGGAGGGGGCTTGCCCGCGAATACGTTGGTGGATGCCCCGCCGCTTTCGCGGGAATGCCCGCTCCCACAGTGTTCGGCGCAAGCCACgg                                                                                                          <  2:2761080/150‑1 (MQ=255)
                                                  gcttGCGCCACTTCCCCCACTATTTACGCGGAGCCCGCGCCACGTCGAGATCCTGTGGGCGCGGGCTTGCCCGCGAATACTTTGGTGGATCCACCGCCGCATTCGGGGGAACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCa                                                                                                     <  2:2423364/147‑1 (MQ=255)
                                                          cACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCGCGGATACGGTGGTGGGTTGAACCCCCCCGTCGCGGGGAAGCCCCCGCCCACAGGGTTGGGCGCAACCCCCCGGGTCACGTggagc                                                                                             >  2:2187077/1‑146 (MQ=255)
                                                            ctgccGCAACTATTTCGGCGGAGCCGGCGGCAGGTCGTGCTCGTGTGGGCGGGGGGTGGCCGGCTAATCCTTTGGTGTATACCCCGCCGCTTGGGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGa                                                                                           <  2:1045567/147‑1 (MQ=255)
                                                             ttCCGCAAATAATTACGCGTATCCCGAGGAAACTCCATATCATGTGGGAGCGGGCTTGCCCGCGAATACTATGTTCTATTAACCGCCGAATTCGCTGGCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                                          <  2:1410216/150‑1 (MQ=255)
                                                             tcccGCCACTAATTACGGGGACCCCGAGCAATGTCTATATCCTGTGGGAGCGGGCTTGCCCGCCATTACGTTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                                          <  2:8623658/148‑1 (MQ=37)
                                                             tcccGCCACTAATCCCGCGGACCCCGCGCCATGTCTCTACCGTGGGGGTGCGGGCGGGCCCCCGAACACTTGGGTGTATTCCCCGCCGCATGCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                                          <  2:830038/148‑1 (MQ=255)
                                                                 gCCACTAATTCCGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGTTGGTGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGttt                                                                                      <  2:2071378/150‑1 (MQ=255)
                                                                  ccACTATTTCCGGGTTTCCCGTGGCATGTCGTGATCCTGGGGGTGGGGGATTGCCCGCTATTACGTGGGTGTATTCACCGCCGCTTTCGCGGCCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt                                                                                     <  1:7601324/150‑1 (MQ=255)
                                                                  ccACTAATTCCGCGTATCCCGCGGAAAGTCTAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGTTGGTGGATAAACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt                                                                                     <  2:3974726/150‑1 (MQ=37)
                                                                  ccACTAATTACGCGGACCCCGCGGCACGGCGAGATCCTGTGGGCGGGGGCTTGCCCGCGAATACGAGGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt                                                                                     <  2:2864952/150‑1 (MQ=38)
                                                                    aCTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcg                                                                                                                                                  >  1:5056169/1‑87 (MQ=255)
                                                                    aCTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcg                                                                                                                                                  <  2:5056169/87‑1 (MQ=255)
                                                                         ttACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGCCCGCGCATAAGGTGGTGGATTGACAGCCGCAGTAGCGGGGAAGCCCGCACCCCCACTGTTGGGCGCAACCCCCCGGGGAACATGAATAGAGGgtttttaaaga                                                                              >  2:2155733/1‑143 (MQ=255)
                                                                           aCGGGGAGCCCGCGGCAAGTCGAGTGCCTGTGGGAGGGGGCTGGCCGGCGAATACGTTGGTTGTTCCCCCGCCGCATGCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                                            <  2:8501605/150‑1 (MQ=255)
                                                                            cgcgTATCCCGAGGCAAGTCGAGAGCCTGTGGGAGCGGGCTTGCCCGCGAATACGAGGGTGGATCCCCCGCCGCATTCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAg                                                                           <  2:7111869/150‑1 (MQ=38)
                                                                                taTCCCGAGGCACGCGTAGATCTTGTGGGTGCGGGCTGGCCGGCTATTAGTTTGGTGTATACACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                                       <  2:7786258/149‑1 (MQ=255)
                                                                                gtgcccGCGGCAAGTCGAGATCGTGTGGGTGCGGGCTTGCCCGCGAATAGGTTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                                       <  2:8767718/147‑1 (MQ=38)
                                                                                gagcccGCGGCACGGCGAGATCCGGGGGGAGCGGGCTGGCCCGCGAATCCGATGGTGGATGCACCGCCGCATTGGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                                       <  2:7187565/147‑1 (MQ=255)
                                                                                gATCCGGCGGCAAGTCGATATCGTGTGGGAGCGGGCTTGCCCGCTAATAGTTTGGTGGATTCACCGCCGCTTTCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                                       <  1:8389237/150‑1 (MQ=255)
                                                                                gATCCCGAGTCAAGTCGATTTCCTGTGTGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGTCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                                       <  2:5311025/150‑1 (MQ=255)
                                                                                  tCCCGAGCCAAGCCGAGATCCTGTGGGAGCGGGCTCCCCCGCTAATACGATGTTGCATACACCCCCCCATACGCGGCCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAACTTTTGa                                                                     <  2:1336967/150‑1 (MQ=255)
                                                                                     cGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCGCGAAAACGATTGTGGTTTAAACGCCCCCTTCGCGGGCAAGCCCCCTCCCACCGTGTTGGGCGCAACCCACGGGCTCACGTGAATCGAGGGTTTGACATGAAAGgttttgtca                                                                  >  2:6699595/1‑145 (MQ=255)
                                                                                       cggCATGGCGATATCCTGTGGGCGCGGGCTTGCCCGCGAATACGATGGGGGCTCCCCCGCCGCATTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcatc                                                                <  2:2695577/149‑1 (MQ=255)
                                                                                       cggCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGGGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcatc                                                                <  2:8242247/149‑1 (MQ=255)
                                                                                        ggCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCg                                                               <  2:8996042/150‑1 (MQ=255)
                                                                                            aGTCGATATCGTGTGGGCGGGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAg                                                           <  1:9040120/150‑1 (MQ=255)
                                                                                                     tcgtgtgGGCGGGGGCTGGCCCGCTAACCGGCGGGTGTTTCCCCCGCCGCCTGCGGGGGAACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGcc                                                   <  2:769732/147‑1 (MQ=255)
                                                                                                         gtgGGCGCGGGCTGGCCCGCGAATACGTTGGTGGATACACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                               <  2:4340556/150‑1 (MQ=255)
                                                                                                                  ggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcct                                      >  1:8115148/1‑150 (MQ=255)
                                                                                                                   gCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg                                     >  1:5255702/1‑150 (MQ=255)
                                                                                                                   gCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg                                     >  2:4753249/1‑150 (MQ=255)
                                                                                                                     ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAGGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                                   >  1:217508/1‑150 (MQ=255)
                                                                                                                        cccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACCTGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg                                >  1:6251987/1‑150 (MQ=255)
                                                                                                                        cccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg                                >  1:6081923/1‑150 (MQ=255)
                                                                                                                        cccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg                                >  2:6725530/1‑150 (MQ=255)
                                                                                                                         ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCAGTTGAATCGAGGTTTTGAGATGTAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcgtgg                               >  2:8852832/1‑150 (MQ=25)
                                                                                                                         ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  1:5449636/1‑150 (MQ=37)
                                                                                                                         ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  2:1481106/1‑150 (MQ=37)
                                                                                                                         ccGCGAATGCGGCGGTGAATCCACCATCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  2:1177548/1‑150 (MQ=37)
                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  2:1349288/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  1:421727/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  2:2673622/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  1:874351/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  1:6683563/1‑150 (MQ=255)
                                                                                                                          cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa                              >  2:6915759/1‑150 (MQ=37)
                                                                                                                           gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                             >  2:7873371/1‑150 (MQ=255)
                                                                                                                           gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                             >  1:4915124/1‑150 (MQ=255)
                                                                                                                           gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                             >  2:765616/1‑150 (MQ=255)
                                                                                                                           gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                             >  2:765555/1‑150 (MQ=255)
                                                                                                                            cgAATACGATGGTGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                            <  2:566649/150‑1 (MQ=255)
                                                                                                                             gAACACGTTGGTGCTTTCCCCGCCGATTTGGGGGGCATGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCTAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGATCACGCCCTGGCCTGGAAAg                           <  1:8364020/150‑1 (MQ=255)
                                                                                                                             ctattACGTTGTTGTATCCCCCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                            <  1:2808808/146‑1 (MQ=255)
                                                                                                                                aTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc                         >  1:6481258/1‑150 (MQ=255)
                                                                                                                                 tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCg                        >  1:34875/1‑150 (MQ=255)
                                                                                                                                  gcggcggtgaATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGc                       >  2:3764540/6‑150 (MQ=38)
                                                                                                                                  aCGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTTTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGc                       >  1:326751/1‑150 (MQ=255)
                                                                                                                                    gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc                     >  1:1451992/1‑150 (MQ=255)
                                                                                                                                    gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc                     >  2:5130159/1‑150 (MQ=255)
                                                                                                                                    gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc                     >  1:1698933/1‑150 (MQ=255)
                                                                                                                                        gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt                                                                                     <  1:3072986/82‑1 (MQ=255)
                                                                                                                                        gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt                                                                                     >  2:3072986/1‑82 (MQ=255)
                                                                                                                                        gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc                                           <  1:7971700/124‑1 (MQ=255)
                                                                                                                                        gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGa                 >  2:499284/1‑150 (MQ=255)
                                                                                                                                        gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGa                 >  2:5073003/1‑150 (MQ=255)
                                                                                                                                         tggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAc                >  2:4186433/1‑150 (MQ=255)
                                                                                                                                           gATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAccc              >  1:2404380/1‑150 (MQ=255)
                                                                                                                                           aatccACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCATGAGCCCGCCCTGGCCTGGACAGTCGCCCAGGAccc              >  2:2579004/5‑150 (MQ=255)
                                                                                                                                            aTTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCa             >  2:5890776/1‑150 (MQ=255)
                                                                                                                                             ttCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAc            >  2:6997564/1‑150 (MQ=255)
                                                                                                                                                 ccgccgCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGt        <  1:6062760/150‑1 (MQ=255)
                                                                                                                                                  cgccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTc       >  2:6756803/1‑150 (MQ=255)
                                                                                                                                                  cgccgcATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCActttag       <  2:6314272/150‑6 (MQ=255)
                                                                                                                                                   gccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCg      >  2:417258/1‑150 (MQ=255)
                                                                                                                                                   gccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCg      >  2:530883/1‑150 (MQ=255)
                                                                                                                                                    ccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCga     >  2:1905048/1‑150 (MQ=255)
                                                                                                                                                      gcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCgaga   <  2:7316367/150‑1 (MQ=255)
                                                                                                                                                       cATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGaa  >  2:6062436/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
TGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAA  >  NC_002947/5486992‑5487290

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