Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,903,526 G→C 9.0% intergenic (+51/‑220) PP_1703 → / → nirB assimilatory nitrate reductase/sulfite reductase/nitrite reductase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,903,5260GC 9.0% 159.6 / 5.1 55intergenic (+51/‑220)PP_1703/nirBassimilatory nitrate reductase/sulfite reductase/nitrite reductase large subunit
Reads supporting (aligned to +/- strand):  ref base G (27/23);  new base C (3/2);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

TGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCC  >  NC_002947/1903378‑1903668
                                                                                                                                                     |                                                                                                                                                
tGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTgggcgg                                                                                                                                                  <  2:38258/149‑1 (MQ=255)
       cGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAgc                                                                                                                                          >  1:108036/1‑149 (MQ=255)
        gCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCTATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAgcg                                                                                                                                         <  3:343280/149‑1 (MQ=255)
        gCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAgcg                                                                                                                                         >  7:323939/1‑149 (MQ=255)
        gCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCCCGC‑CCAgc                                                                                                                                          >  6:13383/1‑148 (MQ=255)
         ccGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAgcg                                                                                                                                         >  2:9490/1‑148 (MQ=255)
           gAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCaa                                                                                                                                      <  6:146025/149‑1 (MQ=255)
           gAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCaa                                                                                                                                      <  2:172971/149‑1 (MQ=255)
           gAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCCCGC‑CCAgccaat                                                                                                                                      <  5:13383/149‑5 (MQ=255)
            aGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGACGGCGTGAGCAAGGCCAGGCGTTATCGGGGCGATGTGTACTGAAGTGCTCCAAGATGGGGGGTTGTGCACCTTGAATGTGC‑ATTGGCGGGGCCCGC‑CCAgccaatg                                                                                                                                     <  1:376165/149‑6 (MQ=255)
              cATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAggg                                                                                                                                   >  5:334648/1‑147 (MQ=255)
                  gCCGGTTGGATGCGGCGGCGGGCGCGGGTTATTTCGAAGGCTTGCTCAAGGCCAGGCGGTATCGGCGCGATGTGGACTGCACGGCCCCAAGTTGGGGTGTTGTGCACCTTGTATGTGC‑ATTGGCTGGGCGGGC‑CCATCGGAAGTGTAgg                                                                                                                               <  2:178218/149‑1 (MQ=255)
                   ccGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGGTGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAg                                                                                                                                >  8:119000/1‑147 (MQ=255)
                     ggatggatGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CAAGCGCAAGTGTAGGAGc                                                                                                                            >  5:177736/1‑149 (MQ=255)
                          gatgCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGt                                                                                                                        <  5:287354/148‑1 (MQ=255)
                          gatgCGGCGGCGGCTGCGGCCTATCGCGCACGCGTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGGTGGGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCgg                                                                                                                          <  8:277648/146‑1 (MQ=255)
                           atgCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTgggcgg                                                                                                                                                  >  1:209954/1‑122 (MQ=255)
                           atgCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTgggcgg                                                                                                                                                  <  2:209954/122‑1 (MQ=255)
                             gcggcggcgGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGGTGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTa                                                                                                                     <  7:119000/148‑1 (MQ=255)
                                   gcggctgcggctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGa                                                                                                              <  4:134033/149‑1 (MQ=255)
                                       ctgcggctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAagag                                                                                                          >  7:196404/1‑149 (MQ=255)
                                         gcggctTATGTCGAAGGCTTGAGCATCGGCACGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                        <  6:177736/149‑1 (MQ=255)
                                          cggctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAgg                                                                                                         <  1:268233/147‑1 (MQ=255)
                                                aTGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTAATGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGaca                                                                                                   <  6:334648/147‑1 (MQ=255)
                                                 tGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcag                                                                                                <  7:7480/149‑1 (MQ=255)
                                                   tCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACACGCGAAGAGGCCGACAcagg                                                                                               >  2:306051/1‑148 (MQ=255)
                                                         gCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcagg‑caaccca                                                                                       <  1:265718/149‑1 (MQ=255)
                                                         gCTTGAGCAAGGCCAGGCGTTATCGGAGCGATGTGTACTGAAGTGCTACAAGATGGGGTGTTGTGCCCCTGGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACACGCGAAGAGGCCGACAcagg‑caacc                                                                                         >  2:354987/1‑147 (MQ=255)
                                                           ttGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCACCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcagg‑caacccagg                                                                                     >  7:296227/1‑149 (MQ=255)
                                                              aGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcagg‑caacccaggcaa                                                                                  >  6:190122/1‑149 (MQ=255)
                                                                 aaGGCCAGGGGTTATCGGCGCGTTGTGTACTGAAGTGCTACAAGTTGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGCTCTAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcagg‑caacccaggcaac                                                                                 <  1:306051/148‑1 (MQ=38)
                                                                    gCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACg                                                                            <  8:296227/149‑1 (MQ=255)
                                                                     ccAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACg                                                                            >  2:20103/1‑148 (MQ=255)
                                                                         gCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAAcc                                                                       <  7:91108/149‑1 (MQ=255)
                                                                            ttATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGa                                                                    >  7:214668/1‑149 (MQ=255)
                                                                               tCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTg                                                                 <  2:108036/149‑1 (MQ=255)
                                                                                  gcgcGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCa                                                              >  7:213776/1‑149 (MQ=255)
                                                                                   cgcgATGAGTACTGAACTGCTCCAAGCTGGCGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCAc                                                             <  8:323939/149‑1 (MQ=255)
                                                                                            tACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcagg‑caacccaggc                                                                                    <  8:138129/117‑1 (MQ=255)
                                                                                            tACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcagg‑caacccaggc                                                                                    >  7:138129/1‑117 (MQ=255)
                                                                                               tGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCtt                                                  >  2:94627/1‑148 (MQ=255)
                                                                                                gAAGTGCTCCAAGATGGGGTGTTGTGCCCCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGCCACAGG‑CTACCCAGGCAACCCAAGCAACCCGATTGGCACGCTGCCTGCTTaa                                                >  7:263405/1‑149 (MQ=39)
                                                                                                  aGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGCCACCGG‑CAACCCAGGCAACCCACGCAGCCCGATTGGCACGCTCCCTGCTTAAcc                                              >  5:300/1‑149 (MQ=255)
                                                                                                    tGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCAGATTGGCACGCTCCCTGCTTAACCTc                                            <  1:9490/149‑1 (MQ=255)
                                                                                                      cTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTACt                                                   <  6:255376/140‑1 (MQ=255)
                                                                                                      cTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTACt                                                   >  5:255376/1‑140 (MQ=255)
                                                                                                       tCCAAGATGGGGTGTTGTGCACCTTGAATGTGCAATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCg                                           >  1:249216/1‑148 (MQ=255)
                                                                                                       tCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGa                                         >  7:292777/1‑149 (MQ=255)
                                                                                                       tCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGa                                         >  7:337973/1‑149 (MQ=255)
                                                                                                         cAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACCCAGG‑CACCCCCGGCACCCCACGCCACCCGCTTGGCACGCTCCCTGCGAAACCTCGGaga                                       >  8:284150/1‑147 (MQ=38)
                                                                                                           aGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATaca                                     <  5:190122/149‑1 (MQ=255)
                                                                                                             aTGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATacaac                                   <  8:213776/149‑1 (MQ=255)
                                                                                                                       gtgCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACt                         >  8:226540/1‑149 (MQ=255)
                                                                                                                                     tgC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGTGGGTTTACCCGCGAAGAGGCCGACACAGGAGAGCCCAGGGAACCCACGCGACCCGATTTGCACGCTCCCGGCTGTAACGCGGGGACCACAAAGCCCACTGATAAACGGCGAt            >  1:315153/1‑149 (MQ=25)
                                                                                                                                       acttgcgCTGGGCCCGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATc           >  2:147580/7‑148 (MQ=255)
                                                                                                                                          ttGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAAcc       >  6:369203/1‑149 (MQ=255)
                                                                                                                                           gcgCTGGGCCCGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACGACAAAGCCCACTGATCAACGGCGATCAACCg      >  8:281433/3‑149 (MQ=255)
                                                                                                                                               tGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGcc  >  8:114006/1‑149 (MQ=255)
                                                                                                                                                 ggcgggc‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGa             <  8:163785/136‑1 (MQ=255)
                                                                                                                                                 ggcgggc‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGa             >  7:163785/1‑136 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGC‑ATTGGCTGGGCGGGC‑CCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGG‑CAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCC  >  NC_002947/1903378‑1903668

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