Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 997,221 C→T 6.7% L126L (CTG→CTA)  PP_0860 ← sulfite reductase flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947997,2210CT 6.7% 210.0 / 3.1 70L126L (CTG→CTA) PP_0860sulfite reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base T (3/2);  total (37/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 = 1.00e+00

TGGCCCAGTCGAAGGTCAGCCAGGTGCGCCAGTCCAGCGCTTTGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGGGCGGGCATGAAGAAGACCCGTGCAGCGTTGCCTTCGCGGACATCGA  >  NC_002947/997073‑997369
                                                                                                                                                    |                                                                                                                                                    
tGGCCCAGTCGAAGGTCAGCCAGGTGCGCCAGTCCAGCGCTTTGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCa                                                                                                                                                     <  2:729937/150‑1 (MQ=255)
    ccAGTCGAAGGTCAGCCAGGTGCGCCAGTCCAGCGCTTTGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCtt                                                                                                                                                 <  2:303378/150‑1 (MQ=255)
       gTCGAAGGTCAGCCAGGTGCGCCAGTCCAGCGCTTTGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGTAGCTTGAg                                                                                                                                              >  1:818453/1‑150 (MQ=255)
        tCGAAGGTCAGCCAGGTGCGCCAGTCCAGCGCTTTGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGc                                                                                                                                             <  1:905807/150‑1 (MQ=255)
         cGAAGGTCAGCCAGGTGCGCCAGTCCAGCGCTTTGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCa                                                                                                                                            >  2:480463/1‑150 (MQ=255)
          gAAGGTCAGCCAGGTGCGCCAGTCCAGCGCTTTGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCAt                                                                                                                                           >  2:1006754/1‑150 (MQ=255)
          gAAGGTCAGCCAGGTGCGCCAGTCCAGCGCTTTGCGCAGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCAt                                                                                                                                           >  1:1115461/1‑150 (MQ=255)
            aGGTCAGCCAGGTGCGCCAGTCCAGCGCTTTGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGa                                                                                                                                         <  1:186361/150‑1 (MQ=255)
               tCAGCCAGGTGCGCCAGTCCAGCGCTTTGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGt                                                                                                                                      >  2:478151/1‑150 (MQ=255)
                cAGCCAGGTGCGCCAGTCCAGCGCTTTGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGtt                                                                                                                                     <  1:1117220/150‑1 (MQ=255)
                       gTGCGCCAGTCCAGCGCTTTGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGaagaag                                                                                                                              <  1:998925/150‑1 (MQ=255)
                        tGCGCCAGTCCAGCGCTTTGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGc                                                                                                                             <  1:469947/150‑1 (MQ=255)
                           gcCAGTCCAGCGCTTTGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCtt                                                                                                                          >  1:1134984/1‑150 (MQ=255)
                            ccAGTCCAGCGCTTTGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTg                                                                                                                         >  1:278444/1‑150 (MQ=255)
                                 ccAGCGCTTTGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCgg                                                                                                                    <  2:685778/150‑1 (MQ=255)
                                        tttGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGtt                                                                                                                                     >  2:623375/1‑126 (MQ=255)
                                        tttGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGtt                                                                                                                                     <  1:623375/126‑1 (MQ=255)
                                          tGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCAc                                                                                                           <  2:278444/150‑1 (MQ=255)
                                           gcgcGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCAcc                                                                                                          <  1:1012633/150‑1 (MQ=255)
                                            cgcagCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCt                                                                                                         <  2:1115461/146‑1 (MQ=255)
                                                        aGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGccc                                                                                             <  1:617903/150‑1 (MQ=255)
                                                        aGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCAGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGccc                                                                                             >  2:231267/1‑150 (MQ=255)
                                                          gTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCttt                                                                                                       <  1:433196/138‑1 (MQ=255)
                                                          gTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCttt                                                                                                       >  2:433196/1‑138 (MQ=255)
                                                                gCCGGACTGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACgg                                                                                     >  1:501978/1‑150 (MQ=255)
                                                                   tgACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCATCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCAGAGTGACCTTTGAGTTCGCCCGTGTACGGGTc                                                                                  <  1:635363/149‑1 (MQ=255)
                                                                     aCAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCATCTTGAGCATGAGGTTGGGGAAGCCGTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCga                                                                                >  2:635363/1‑149 (MQ=255)
                                                                         gCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGtt                                                                                                  <  1:682575/128‑1 (MQ=255)
                                                                         gCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGtt                                                                                                  >  2:682575/1‑128 (MQ=255)
                                                                         gCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTAGGCATa                                                                            >  1:851444/1‑150 (MQ=255)
                                                                         gCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTACCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATa                                                                            >  2:718293/1‑150 (MQ=255)
                                                                          cAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGTAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAg                                                                           <  2:818453/150‑1 (MQ=255)
                                                                           agaagaACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGTAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAgc                                                                          >  2:944594/1‑150 (MQ=255)
                                                                           agaagaACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCAAGTATCGCACCCGAAGAGGAAACGGCATTGCAGGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAgc                                                                          >  1:1052418/1‑150 (MQ=255)
                                                                                  cacCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTc                                                                   >  2:843658/1‑150 (MQ=255)
                                                                                   acCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCCTGTACGGGTCGGCATAGCGCAGCTCg                                                                  <  1:351776/150‑1 (MQ=255)
                                                                                    ccAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCAc                                                                                                           <  1:695854/108‑1 (MQ=255)
                                                                                    ccAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCAc                                                                                                           >  2:695854/1‑108 (MQ=255)
                                                                                     cAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGGGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGTTCGGCATAGCGCAGCTCGcc                                                                >  2:109324/1‑150 (MQ=255)
                                                                                     cAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGcc                                                                >  1:392895/1‑150 (MQ=255)
                                                                                     cAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGcc                                                                >  1:372911/1‑150 (MQ=255)
                                                                                              ggTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGTAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTc                                                       <  1:944594/150‑1 (MQ=255)
                                                                                               gTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCCCgggcgaac                                                      <  2:1052418/150‑7 (MQ=255)
                                                                                                    ggCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGGGTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCCTTGAGTTCGCCCGGGTACGGATCGGCATAACGCAGCTCGACACGGCGGTCGCACgg                                                 >  2:64276/1‑150 (MQ=255)
                                                                                                       tCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGATCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGgggc                                              >  1:861220/1‑150 (MQ=255)
                                                                                                         cGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGt                                                         >  2:153662/1‑137 (MQ=255)
                                                                                                         cGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGt                                                         <  1:153662/137‑1 (MQ=255)
                                                                                                         cGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGgggcgg                                            <  2:1063374/150‑1 (MQ=255)
                                                                                                          agTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGTTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGgggcggg                                           >  2:830829/2‑150 (MQ=255)
                                                                                                                   ccGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGGGCGGGCATgaagaa                                  <  1:1199446/150‑1 (MQ=255)
                                                                                                                     ggCCAGTATCGCCCATGGCGAGGAAACGGTGCAGc                                                                                                                                                   <  1:155798/35‑1 (MQ=255)
                                                                                                                     ggCCAGTATCGCCCATGGCGAGGAAACGGTGCAGc                                                                                                                                                   >  2:155798/1‑35 (MQ=255)
                                                                                                                     ggCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGGGCGGGCATGAagaaga                                <  1:718293/150‑1 (MQ=255)
                                                                                                                          gTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGGGCGGGCATGAAGAAGACCCGt                           <  1:480463/150‑1 (MQ=255)
                                                                                                                             tCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGACTTCGGCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGGGCGGGCATGAAGAAGACCCGTGCa                        >  1:85777/1‑150 (MQ=255)
                                                                                                                               gCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGGGCGGGCATGAAGAAGACCCGTGCAGc                      <  2:596957/150‑1 (MQ=255)
                                                                                                                                  cATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGGGCGGGCATGAAGAAGACCCGTGCAGCGtt                   <  1:583707/150‑1 (MQ=255)
                                                                                                                                   aTGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGGGCGGGCATGAAGAAGACCCGt                           <  2:396939/141‑1 (MQ=255)
                                                                                                                                   aTGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGGGCGGGCATGAAGAAGACCCGt                           >  1:396939/1‑141 (MQ=255)
                                                                                                                                    tGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGTGCGGGCATGAAGAAGACCCGTGCAGCGTTGc                 <  1:231267/150‑1 (MQ=255)
                                                                                                                                           gAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGc                                                                     >  2:829607/1‑91 (MQ=255)
                                                                                                                                           gAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGc                                                                     <  1:829607/91‑1 (MQ=255)
                                                                                                                                           gAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGGGCGGGCATGAAGAAGACCCGTGCAGCGTTGCCTTcgcg          >  1:806689/1‑150 (MQ=255)
                                                                                                                                             aaCGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGGGCGGGCATGAAGAAGACCCGTGCAGCGTTGCCTTCGCGGa        >  2:701056/1‑150 (MQ=255)
                                                                                                                                             aaCGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGGGCGGGCATGAAGAAGACCCGTGCAGCGTTGCCTTCGCGGa        >  1:1006091/1‑150 (MQ=255)
                                                                                                                                              aCGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGGGCGGGCATGAAGAAGACCCGTGCAGCGTTGCCTTCGCGGAc       >  2:431023/1‑150 (MQ=255)
                                                                                                                                              aCGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGGGCGGGCATGAAGAAGACCCGTGCAGCGTTGCCTTCGCGGAc       >  2:431003/1‑150 (MQ=255)
                                                                                                                                                ggTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGgggcgg                                            >  2:593819/1‑111 (MQ=255)
                                                                                                                                                ggTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGgggcgg                                            <  1:593819/111‑1 (MQ=255)
                                                                                                                                                   gCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGGGCGGGCATGAAGAAGACCCGTGCAGCGTTGCCTTCGCGGACATCGa  >  1:1184179/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGCCCAGTCGAAGGTCAGCCAGGTGCGCCAGTCCAGCGCTTTGCGCGGCCAGCGCAGGTACAGGCCGGACAGGCAGAAGAACACCAGCATCAGGGTGCAGGCTCCGGTGATTTGCCGGCCAGTATCGCCCATGGCGAGGAAACGGTGCAGCTTGAGCATGAGGTTGAAGAAGCCTTGCCCGGCCGCTTCACCTTTGAGTTCGCCCGTGTACGGGTCGGCATAGCGCAGCTCGCCACGGCGTTCGCCCGGGGCGGGCATGAAGAAGACCCGTGCAGCGTTGCCTTCGCGGACATCGA  >  NC_002947/997073‑997369

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: ‑