Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,873,747 T→A 7.9% intergenic (‑128/+30) 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,7470TA 7.9% 183.2 / 5.5 63intergenic (‑128/+30)aqpZ/PP_4283aquaporin Z/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base A (2/3);  total (27/36)
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.75e-01

ATGCGTTTACCCAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTGGTCCATGTGATGCA  >  NC_002947/4873600‑4873895
                                                                                                                                                   |                                                                                                                                                      
aTGCGTTTACCCAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCtg                                                                                                                                                       >  8:5265/1‑148 (MQ=255)
          ccAGGGCCATGCTCATGAATCCATTCCTCAAGAAGATGAACGCAACGAATGCGCAAGCCCAAATGATTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGTGCTGCCTTGCAGCCCCTc                                                                                                                                             >  4:182196/1‑149 (MQ=255)
           cAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCtt                                                                                                                                                       >  3:170029/1‑138 (MQ=255)
           cAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCtt                                                                                                                                                       <  4:170029/138‑1 (MQ=255)
             gggACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCAAGGCAGCCCTTccag                                                                                                                                          <  3:343384/149‑4 (MQ=255)
               gACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCAAGGCAGccct                                                                                                                                               >  5:120865/1‑141 (MQ=255)
               gACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCAAGGCAGccct                                                                                                                                               <  6:120865/142‑2 (MQ=255)
                     cTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAg                                                                                                                                 <  5:85943/148‑1 (MQ=255)
                     cTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAg                                                                                                                                 <  8:322928/148‑1 (MQ=255)
                               cATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATc                                                                                                                       >  5:78830/1‑148 (MQ=255)
                               cATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATc                                                                                                                       >  2:135553/1‑148 (MQ=255)
                               cATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATc                                                                                                                       <  6:78830/148‑1 (MQ=255)
                                aTTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGa                                                                                                                                          >  4:36653/1‑129 (MQ=255)
                                aTTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGa                                                                                                                                          <  3:36653/129‑1 (MQ=255)
                                aTTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCaa                                                                                                                     <  7:173663/149‑1 (MQ=255)
                                aTTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCaa                                                                                                                     <  1:186639/149‑1 (MQ=255)
                                    cTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAAccc                                                                                                                  <  1:91292/148‑1 (MQ=255)
                                     tCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCAAGGCAGCCCTTc                                                                                                                                             >  4:62705/1‑122 (MQ=255)
                                     tCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCAAGGCAGCCCTTc                                                                                                                                             <  3:62705/122‑1 (MQ=255)
                                       aagaagTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCtt                                                                                                                <  5:65388/147‑1 (MQ=255)
                                       aagaagTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGa                                                                                                              >  3:183143/1‑149 (MQ=255)
                                            gTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCa                                                                                                                                      >  1:153462/1‑120 (MQ=255)
                                            gTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCa                                                                                                                                      <  2:153462/120‑1 (MQ=255)
                                                 aCGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCtt                                                                                                    <  5:33006/149‑1 (MQ=255)
                                                  cGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTg                                                                                                   <  6:233310/149‑1 (MQ=255)
                                                   gCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCaa                                                                                                                     <  2:222458/130‑1 (MQ=255)
                                                   gCAACGAATGCGCAAGCCCAACTGCTTGATGAAAGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCaa                                                                                                                     >  1:222458/1‑130 (MQ=255)
                                                     aaCGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGt                                                                                                >  2:289705/1‑149 (MQ=255)
                                                     aaCGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGt                                                                                                >  2:213832/1‑149 (MQ=255)
                                                          aTGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATc                                                                                                                       >  2:89010/1‑121 (MQ=255)
                                                          aTGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATc                                                                                                                       <  1:89010/121‑1 (MQ=255)
                                                             cgcAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCaa                                                                                        <  4:183143/149‑1 (MQ=255)
                                                                      aaCTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCGAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCTATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTgg                                                                                                  <  6:286220/130‑1 (MQ=255)
                                                                      aaCTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTgg                                                                                                  >  5:286220/1‑130 (MQ=255)
                                                                      aaCTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCga                                                                               <  5:326139/149‑1 (MQ=255)
                                                                      aaCTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCga                                                                               <  7:104354/149‑1 (MQ=255)
                                                                      aaCTGCTTGATGAAGGTTTGCTCAAAAGAATTTAGCAGACCTTTGGGTTATTGCAATTAGTGGCTGGAATGGCTTCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCga                                                                               <  3:182196/149‑1 (MQ=255)
                                                                      aaCGGCTTGATGAAGGTTTGCTCACAAGATTTGAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCGTGCAGCCCCTCGATCCCGCAGGCATACGATCCACCCTTGACGGAAACCTTGGGTTTTTTCAAGAGATGCga                                                                               <  4:138770/149‑1 (MQ=255)
                                                                         tGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGaa                                                                            >  7:291440/1‑149 (MQ=255)
                                                                            ttGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAacac                                                                         >  6:313691/1‑149 (MQ=255)
                                                                                 gaAGGTTTGCTCACAAGAATTAAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGCAAGGGCTGCCTTGCAGACACTCGATCAAGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTATTTTCAACAGATGCGAGAACACGGCAt                                                                    >  1:263015/1‑149 (MQ=255)
                                                                                           tCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAAc                                                                                       <  2:124816/120‑1 (MQ=255)
                                                                                             acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGaa                                                                                                         <  7:245180/100‑1 (MQ=255)
                                                                                             acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGaa                                                                                                         >  8:245180/1‑100 (MQ=255)
                                                                                             acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTcc                                                         >  5:26372/1‑148 (MQ=255)
                                                                                             acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTcc                                                         >  6:244688/1‑148 (MQ=255)
                                                                                               aaGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAg                                                      <  1:213832/149‑1 (MQ=255)
                                                                                                        gCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTctcc                                             <  8:291440/149‑1 (MQ=255)
                                                                                                         cAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCga                                                                               <  7:318269/114‑1 (MQ=255)
                                                                                                         cAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCga                                                                               >  8:318269/1‑114 (MQ=255)
                                                                                                         cAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTctcct                                            >  3:262689/1‑149 (MQ=255)
                                                                                                          aGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTctcctc                                           >  8:310479/1‑149 (MQ=255)
                                                                                                                  gggTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGa                                    <  5:244688/148‑1 (MQ=255)
                                                                                                                          ggcAGTTGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACTCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTg                           <  2:263015/147‑1 (MQ=255)
                                                                                                                          gCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTg                           <  1:135553/149‑1 (MQ=255)
                                                                                                                           ccAGTATTGGCTGGAATGGCTGCCTTTGAGCCCCTCGCACTCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCtt                            <  1:289705/149‑1 (MQ=255)
                                                                                                                              gtagtGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGt                       >  7:324276/1‑149 (MQ=255)
                                                                                                                                   ggCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGAtg                  >  1:316016/1‑149 (MQ=255)
                                                                                                                                    gCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGAtg                  <  2:316016/148‑1 (MQ=255)
                                                                                                                                     cTGGAAGGGCTGCCTTGCAGCCCCTCGATCAGGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGAtgtg                >  1:213012/1‑149 (MQ=255)
                                                                                                                                     cTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGAtgtg                <  7:310479/149‑1 (MQ=255)
                                                                                                                                      tGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTgg               >  6:120358/1‑149 (MQ=255)
                                                                                                                                          aGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTGGTCCa           <  6:26372/149‑1 (MQ=255)
                                                                                                                                                 ccTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTGGTCCATGTGATg    <  5:313691/149‑1 (MQ=255)
                                                                                                                                                  cTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAgg                                                     <  8:97787/99‑1 (MQ=255)
                                                                                                                                                  cTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAgg                                                     >  7:97787/1‑99 (MQ=255)
                                                                                                                                                   ttGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAg                                                      >  4:223236/1‑97 (MQ=255)
                                                                                                                                                   ttGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAg                                                      <  3:223236/97‑1 (MQ=255)
                                                                                                                                                   ttGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGCTTATTTCAACAGCTGCGAGAACACCGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGAAGTGGGCCATGTGATGca  >  7:257376/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
ATGCGTTTACCCAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCGATCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTGGTCCATGTGATGCA  >  NC_002947/4873600‑4873895

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: