Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,873,746 C→A 6.2% intergenic (‑127/+31) aqpZ ← / ← PP_4283 aquaporin Z/GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,873,7460CA 6.2% 260.1 / 5.8 81intergenic (‑127/+31)aqpZ/PP_4283aquaporin Z/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (40/35);  new base A (3/2);  total (44/37)
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.81e-01

CCATGCGTTTACCCAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTGGTCCATGTGAT  >  NC_002947/4873598‑4873892
                                                                                                                                                     |                                                                                                                                                  
caaTGCGTTTAACCAGGGAAATGCTCATGAATCCATTCCTAAAGAAGTTGAAAGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGGAGTGGATGGAAGGGCTGcc                                                                                                                                                    >  6:48395/3‑149 (MQ=255)
  aTGCGTTTACCCAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGATGCCtt                                                                                                                                                  >  3:320479/1‑149 (MQ=255)
      gTTTACCCAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAg                                                                                                                                              <  8:38872/149‑1 (MQ=255)
          aCCCAGGGACATGCTCATGACTCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCAAGGCAGccc                                                                                                                                           <  7:285342/148‑1 (MQ=255)
           cccAGGGACATGCTCATGAATCAATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGcccc                                                                                                                                          >  4:24165/1‑148 (MQ=255)
            ccAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCAAGGCAGCCCTTc                                                                                                                                        <  8:174764/149‑1 (MQ=255)
              aGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGa                                                                                                                                      >  5:167726/1‑149 (MQ=255)
                       cTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAg                                                                                                                              <  7:316556/148‑1 (MQ=255)
                        tCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAgg                                                                                                                             >  7:177811/1‑148 (MQ=255)
                          aTGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCa                                                                                                                           >  6:164712/1‑148 (MQ=255)
                             aaTCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAAc                                                                                                                        <  8:233087/148‑1 (MQ=255)
                               tCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGa                                                                                                                      >  8:36033/1‑148 (MQ=255)
                                  aTTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCaa                                                                                                                  >  5:161362/1‑149 (MQ=255)
                                  aTTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCaa                                                                                                                  <  2:297759/149‑1 (MQ=255)
                                      cTCAAGAAGTTGAACTCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAAccc                                                                                                               <  3:24165/148‑1 (MQ=255)
                                       tCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCtt                                                                                                             >  5:103994/1‑149 (MQ=255)
                                       tCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCtt                                                                                                             >  3:303168/1‑149 (MQ=255)
                                       tCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCtt                                                                                                             >  4:8049/1‑149 (MQ=255)
                                       tCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCtt                                                                                                             >  1:87810/1‑149 (MQ=255)
                                         aagaagTTGAACGCAACGAATGCGCAAGCCTAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGACAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGa                                                                                                           >  2:46147/1‑149 (MQ=255)
                                           gaagTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACg                                                                                                         >  1:146534/1‑149 (MQ=255)
                                              gTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGa                                                                                                       >  1:166586/1‑148 (MQ=255)
                                                  aaCGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCt                                                                                                  >  3:79713/1‑149 (MQ=255)
                                                    cGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTg                                                                                                <  7:36033/149‑1 (MQ=255)
                                                     gCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTgg                                                                                               <  2:166586/149‑1 (MQ=255)
                                                       aaCGAATGCGCATGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCTAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGt                                                                                             >  5:110778/1‑149 (MQ=255)
                                                       aaCGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGt                                                                                             >  2:172098/1‑149 (MQ=255)
                                                       aaCGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGt                                                                                             >  5:118982/1‑149 (MQ=255)
                                                            aTGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGtttttt                                                                                        <  6:167726/149‑1 (MQ=255)
                                                            aTGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGttttt                                                                                         >  8:237853/1‑148 (MQ=255)
                                                             tGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTc                                                                                       <  8:177811/149‑1 (MQ=255)
                                                               cgcAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCaa                                                                                     <  4:153286/149‑1 (MQ=255)
                                                                gccatcccAAAGGCTTTATTAATGTTTTCTCAAAACAAAATTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCaa                                                                                     <  1:114739/144‑1 (MQ=255)
                                                                gcAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCaa                                                                                     >  1:136359/1‑148 (MQ=255)
                                                                gcAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCaa                                                                                     >  7:101998/1‑148 (MQ=255)
                                                                gcAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCaa                                                                                     >  3:236975/1‑148 (MQ=255)
                                                                gcAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAAc                                                                                    <  1:224034/149‑1 (MQ=255)
                                                                gcAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAAc                                                                                    <  7:15652/149‑1 (MQ=255)
                                                                        aaCTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTATTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCga                                                                            <  4:79713/149‑1 (MQ=255)
                                                                        aaCTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGc                                                                              >  3:143422/1‑147 (MQ=255)
                                                                                     aGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGc                                                               <  2:156593/149‑1 (MQ=255)
                                                                                            cTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGt                                                        <  2:146534/149‑1 (MQ=255)
                                                                                             tCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTc                                                       >  3:90182/1‑149 (MQ=255)
                                                                                             tCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTc                                                       >  5:92550/1‑149 (MQ=255)
                                                                                               acaAGAATTTAGCAGACCTTTGGGTTATGGCTAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTcc                                                      <  6:110778/148‑1 (MQ=255)
                                                                                               acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTcc                                                      >  8:82258/1‑148 (MQ=255)
                                                                                               acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTcc                                                      >  8:142754/1‑148 (MQ=255)
                                                                                               acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTcc                                                      >  7:294383/1‑148 (MQ=255)
                                                                                               acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTcc                                                      >  3:131663/1‑148 (MQ=255)
                                                                                               acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTcc                                                      <  5:228014/148‑1 (MQ=255)
                                                                                               acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGCGCTGCCTTGCAGCACCTCGATCACGCAGGCAAACGATCACCCCTTGACCGAAACCTTGGCTTTTTTCAACAGATGCGAGAACACTGCATGCAAATCGTcc                                                      >  5:74318/1‑148 (MQ=255)
                                                                                                caAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCGTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTcc                                                      >  1:255947/1‑147 (MQ=255)
                                                                                                caAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTcc                                                      >  4:139730/1‑147 (MQ=255)
                                                                                                 aaGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCATCAGATGCGAGAACACGGCATGCAAATCGTCCGAg                                                   >  6:223042/1‑149 (MQ=255)
                                                                                                 aaGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAg                                                   >  8:5331/1‑149 (MQ=255)
                                                                                                 aaGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAg                                                   <  6:327719/149‑1 (MQ=255)
                                                                                                 aaGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAg                                                   >  5:327719/1‑149 (MQ=255)
                                                                                                 aaGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAg                                                   <  7:82258/149‑1 (MQ=255)
                                                                                                           gCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTctcc                                          <  1:172098/149‑1 (MQ=255)
                                                                                                            cAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTctcct                                         <  3:8049/149‑1 (MQ=255)
                                                                                                             aGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTctcctc                                        >  1:191586/1‑149 (MQ=255)
                                                                                                              gACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTc                                                                                       <  2:199873/101‑1 (MQ=255)
                                                                                                              gACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCg                                       <  4:143422/149‑1 (MQ=255)
                                                                                                                 cTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTcc                                    <  8:101998/149‑1 (MQ=255)
                                                                                                                  tttGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTc                                     >  2:33843/1‑147 (MQ=255)
                                                                                                                     gggTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGa                                 <  6:118982/148‑1 (MQ=255)
                                                                                                                      ggTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGAt                                >  5:100420/1‑148 (MQ=255)
                                                                                                                       gTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGAtt                               >  3:307210/1‑148 (MQ=255)
                                                                                                                        ttATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCa                             <  5:87568/149‑1 (MQ=255)
                                                                                                                             gCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTg                        <  3:139730/149‑1 (MQ=255)
                                                                                                                                 gtagtGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGt                    <  4:303168/149‑1 (MQ=255)
                                                                                                                                  tagtGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTc                   <  6:161362/149‑1 (MQ=255)
                                                                                                                                    gtGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGa                 <  4:131663/149‑1 (MQ=255)
                                                                                                                                     tGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGAt                <  6:100420/149‑1 (MQ=255)
                                                                                                                                       gatcgaggggCTGCAAGGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGAtgt              >  7:267036/8‑149 (MQ=255)
                                                                                                                                        cTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGAtgtg             <  6:103994/149‑1 (MQ=255)
                                                                                                                                        atcgaggggCTGCAAGGCAGCCCTTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGAtgtg             >  4:137741/7‑149 (MQ=255)
                                                                                                                                        atcgaggggCTGCAAGGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGAtgtg             >  7:190590/7‑149 (MQ=255)
                                                                                                                                         tGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTgg            >  7:192717/1‑149 (MQ=255)
                                                                                                                                           gAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTGGTc          <  5:164712/149‑1 (MQ=255)
                                                                                                                                             aGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATAAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTGGTCCa        >  6:321879/1‑149 (MQ=255)
                                                                                                                                                gCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTGGTCCAtgt     <  7:237853/149‑1 (MQ=255)
                                                                                                                                                    ccTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTGGTCCATGTGAt  >  8:117179/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCATGCGTTTACCCAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTGGTCCATGTGAT  >  NC_002947/4873598‑4873892

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