Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,258,051 A→C 6.4% intergenic (+21/‑41) dcd → / → PP_1101 deoxycytidine triphosphate deaminase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,258,0510AC 6.4% 625.2 / 24.0 188intergenic (+21/‑41)dcd/PP_1101deoxycytidine triphosphate deaminase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (91/85);  new base C (4/8);  total (95/93)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 189 total aligned reads displayed.

AGATCTACGCCAACGAAGGCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTGGGTTGTTGCCCTGGTCAC  >  NC_002947/1257903‑1258200
                                                                                                                                                    |                                                                                                                                                     
aGATCTACGCCAACGAAGGCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTAc                                                                                                                                                      <  2:449685/150‑1 (MQ=255)
 gATCTACGCCAACGAAGGCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTAcc                                                                                                                                                     >  1:2616024/1‑150 (MQ=255)
 gATCTACGCCAACGAAGGCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTAcc                                                                                                                                                     >  2:54560/1‑150 (MQ=255)
      aCGCCAACGAAGGCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGTAATTCCCCgcgc                                                                                                                                                <  1:1670731/150‑1 (MQ=255)
        gCCAACGAAGGCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGTAATTCCCCgcgctt                                                                                                                                              <  1:4706668/150‑3 (MQ=255)
        gCCAACGAAGGCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGTAATTCCCCgcgctt                                                                                                                                              <  2:5386698/150‑3 (MQ=255)
          cAACGAAGGCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTTCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCttt                                                                                                                                            >  2:154383/1‑150 (MQ=255)
          cAACGAAGGCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCttt                                                                                                                                            >  2:2043716/1‑150 (MQ=255)
          cAACGAAGGCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCttt                                                                                                                                            <  2:1316195/150‑1 (MQ=255)
          cAACGAAGGCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCttt                                                                                                                                            <  2:590542/150‑1 (MQ=255)
           aaCGAAGGCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTg                                                                                                                                           <  1:1438165/150‑1 (MQ=255)
           aaCGAAGGCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTg                                                                                                                                           <  1:4842924/150‑1 (MQ=255)
           aaCGAAGGCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTg                                                                                                                                           <  1:2470693/150‑1 (MQ=255)
            aCGAAGGCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTgg                                                                                                                                          <  2:4289021/150‑1 (MQ=255)
                  gCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTcc                                                                                                                                    >  1:1590448/1‑150 (MQ=255)
                    gTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCaa                                                                                                                                  >  1:2851030/1‑150 (MQ=255)
                     tGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAAc                                                                                                                                 <  2:1693108/150‑1 (MQ=255)
                     tGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAAc                                                                                                                                 <  1:3802791/150‑1 (MQ=255)
                       gcgcAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGa                                                                                                                               >  2:4594852/1‑150 (MQ=255)
                        cgcAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGaa                                                                                                                              <  2:2365854/150‑1 (MQ=255)
                        cgcAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGaa                                                                                                                              <  1:174177/150‑1 (MQ=255)
                         gcAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAg                                                                                                                             >  1:5044022/1‑150 (MQ=255)
                         gcAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAg                                                                                                                             >  2:1136762/1‑150 (MQ=255)
                           aGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTa                                                                                                                           <  1:3871409/150‑1 (MQ=255)
                           aGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTa                                                                                                                           <  1:3871403/150‑1 (MQ=255)
                              tGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACg                                                                                                                        >  1:5438215/1‑150 (MQ=255)
                               gCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACTTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGc                                                                                                                       <  2:784502/150‑1 (MQ=255)
                               gCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGc                                                                                                                       >  1:3677863/1‑150 (MQ=255)
                                cTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGcc                                                                                                                      <  2:3326474/150‑1 (MQ=255)
                                   ttCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACCAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGgt                                                                                                                   <  2:2820669/150‑1 (MQ=255)
                                      ttGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGtgtg                                                                                                                >  1:3274420/1‑150 (MQ=255)
                                           gTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCatga                                                                                                           >  1:2369107/1‑150 (MQ=255)
                                           gTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCatga                                                                                                           >  1:2574050/1‑150 (MQ=255)
                                              cgacgaAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGAtgt                                                                                                        <  2:215071/150‑1 (MQ=255)
                                              cgacgaAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGAtgt                                                                                                        >  2:4625469/1‑150 (MQ=255)
                                              cgacgaAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCACGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGAtgt                                                                                                        >  1:4315918/1‑150 (MQ=255)
                                               gacgaAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGAtgtg                                                                                                       >  1:4984176/1‑150 (MQ=255)
                                                   aaGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATc                                                                                                   <  1:5298236/150‑1 (MQ=255)
                                                     gagTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCgg                                                                                                 <  1:2658234/150‑1 (MQ=255)
                                                         gCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTc                                                                                                                                     <  2:3561355/110‑1 (MQ=255)
                                                         gCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTc                                                                                                                                     >  1:3561355/1‑110 (MQ=255)
                                                            taGTGTCGTACACGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGt                                                                                          <  2:3894293/149‑1 (MQ=255)
                                                              gtgtCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCa                                                                                        <  1:1916717/150‑1 (MQ=255)
                                                               tgtCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAg                                                                                       <  2:1691708/150‑1 (MQ=255)
                                                               tgtCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAg                                                                                       >  2:4472995/1‑150 (MQ=255)
                                                               tgtCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAg                                                                                       <  1:510148/150‑1 (MQ=255)
                                                               tgtCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAg                                                                                       >  1:1500221/1‑150 (MQ=255)
                                                               tgtCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAg                                                                                       >  1:114245/1‑150 (MQ=255)
                                                                 tCGTACACGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGTGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGa                                                                                     <  1:1023320/150‑1 (MQ=255)
                                                                 tCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGa                                                                                     <  1:1023067/150‑1 (MQ=255)
                                                                 tCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGa                                                                                     >  1:4430776/1‑150 (MQ=255)
                                                                 tCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGa                                                                                     <  2:4295111/150‑1 (MQ=255)
                                                                      cAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATc                                                                                                   >  2:3371217/1‑131 (MQ=255)
                                                                       aaGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGcc                                                                               >  1:1807610/1‑150 (MQ=255)
                                                                       aaGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGcc                                                                               >  1:5253582/1‑150 (MQ=255)
                                                                       aaGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGcc                                                                               >  2:1407215/1‑150 (MQ=255)
                                                                         ggACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCaa                                                                             <  1:2283978/150‑1 (MQ=255)
                                                                         cgACCGTGGCGGCAAATACCAGCGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCAATGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCATTCGTCAGGAGCAGCCaa                                                                             <  2:1117830/149‑1 (MQ=255)
                                                                          gACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAATTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAAt                                                                            <  1:4097150/150‑1 (MQ=255)
                                                                          gACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAAt                                                                            >  1:4446922/1‑150 (MQ=255)
                                                                          gACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAAt                                                                            >  1:3346535/1‑150 (MQ=255)
                                                                              gTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGaaa                                                                        >  2:1346782/1‑150 (MQ=255)
                                                                              gTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGaaa                                                                        >  1:1524113/1‑150 (MQ=255)
                                                                                      aaaTACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATcc                                                                >  1:961391/1‑150 (MQ=255)
                                                                                      aaaTACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACACCATcc                                                                >  1:2065996/1‑150 (MQ=255)
                                                                                        aTACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCg                                                              >  2:4577434/1‑150 (MQ=255)
                                                                                            cAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAAt                                                          >  1:1816155/1‑150 (MQ=255)
                                                                                             aGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATc                                                         >  1:692285/1‑150 (MQ=255)
                                                                                              gggCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCa                                                        >  2:3830346/1‑150 (MQ=255)
                                                                                              gggCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCa                                                        <  2:4712993/150‑1 (MQ=255)
                                                                                                    aCGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCa                                                  >  1:3811596/1‑150 (MQ=255)
                                                                                                    aCGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCa                                                  >  2:2702298/1‑150 (MQ=255)
                                                                                                        ggcgTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCt                                              <  2:670266/150‑1 (MQ=255)
                                                                                                        ggcgTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCt                                              <  1:3352823/150‑1 (MQ=255)
                                                                                                            tCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCa                                          >  1:4558190/1‑150 (MQ=255)
                                                                                                                 cTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCt                                     >  1:4237796/1‑150 (MQ=255)
                                                                                                                 cTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCt                                     >  2:1939433/1‑150 (MQ=255)
                                                                                                                  tGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTg                                    <  2:2574050/150‑1 (MQ=255)
                                                                                                                  tGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTg                                    >  1:5365559/1‑150 (MQ=255)
                                                                                                                       cgTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGcc                               <  2:2031173/150‑1 (MQ=255)
                                                                                                                            cTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATc                          >  2:3519327/1‑150 (MQ=255)
                                                                                                                            cTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATc                          <  1:4123142/150‑1 (MQ=255)
                                                                                                                             tGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGAGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCa                         >  1:2099948/1‑150 (MQ=255)
                                                                                                                                  cGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTg                    <  2:4273681/150‑1 (MQ=255)
                                                                                                                                    aaaggcgcGGGTAATTCCCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTggg                  >  2:2932567/5‑150 (MQ=255)
                                                                                                                                    aCAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTggg                  <  1:2028893/150‑1 (MQ=255)
                                                                                                                                     cAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTGGgt                 <  1:1958961/150‑1 (MQ=255)
                                                                                                                                     cAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTGGgt                 <  2:4331357/150‑1 (MQ=255)
                                                                                                                                      aaGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTGGgtt                <  2:4430776/150‑1 (MQ=255)
                                                                                                                                       ggcgcGGGTAATTCCCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTGGgttg               >  1:477330/2‑150 (MQ=255)
                                                                                                                                        gcgcGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTGGgttgt              <  2:3669474/150‑1 (MQ=255)
                                                                                                                                         cgcgGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTGGgttgtt             <  1:2043712/150‑1 (MQ=255)
                                                                                                                                          gcgGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCACTCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTGGGttgttg            >  1:2778286/1‑150 (MQ=255)
                                                                                                                                           cgGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTGGGTTGTTGc           >  1:3378371/1‑150 (MQ=255)
                                                                                                                                               gAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTGGGTTGTTGCCCTg       <  1:5049748/150‑1 (MQ=255)
                                                                                                                                               gAATTACCCGCGCCTATGGCACGCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGGGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTGGGTTGTTGCCCTg       <  1:1253900/150‑1 (MQ=255)
                                                                                                                                                aaTTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTGGGTTGTTGCCCTgg      <  2:5017031/150‑1 (MQ=255)
                                                                                                                                                aaTTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTGGGTTGTTGCCCTgg      >  1:3206279/1‑150 (MQ=255)
                                                                                                                                                 aTTACCCGCGCCATTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTGGGTTGTTGCCCTGGt     <  2:2616024/150‑1 (MQ=255)
                                                                                                                                                    aCCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTGGGTTGTTGCCCTGGTCAc  >  1:3202435/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGATCTACGCCAACGAAGGCGTGGCGCAGATGCTGTTCTTGCAGTCCGACGAAGAGTGCGAAGTGTCGTACAAGGACCGTGGCGGCAAATACCAGGGCCAACGCGGCGTCACCCTGCCGCGTACCTGAGCCGACAAGCGCGGGGAATTACCCGCGCCTTTGGCACTCCAACGAAGTAACGCCGGTGTGCATGATGTGCATCGGCTTTCGTCAGGAGCAGCCAATGAAACTCCATCCCGCAATCAGTGCCAAGCTGGCAGGCCTGCAGCCCAATCACGTTGGGTTGTTGCCCTGGTCAC  >  NC_002947/1257903‑1258200

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