Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,037 A→G G114G (GGA→GGG creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,637,0370AG96.8% 179.7 / ‑2.3 63G114G (GGA→GGGcreCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (30/31);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

AGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGC  >  NC_000913/4636801‑4637279
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
agagCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTAGCAGGGTAAAGTATTGTTCGATTAGTCAAATAAAGCCGTTGGGCAGGATTATTggcg                                                                                                                                                                                                                                       >  2:572861/1‑251 (MQ=255)
       gtgagggggagttaatcgaacccGCAACGTTGCTGGGGGAGCTGGCGCGTCCCGATGTGCTCTCTGGAGACCCAACGCATGGTCAACTGGCGCAGGCGTTTAATCAGCTACAACAGCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATAGGTCCGATTCGTAGGGTAAAGTATGGTTCGATTCGGCATATAAACCCGTTTGGCAGGATGATTCGCGCTGGAAt                                                                                                                                                                                                                                <  2:279435/230‑1 (MQ=255)
                   taaTCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGGTTTGCTCTCTGGGGACCCAACGCATGGTCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGGTTCGCGCCAATATCGGTTGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTTTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCACGCTGGAATGACGTCTGGcg                                                                                                                                                                                                                     >  1:519507/3‑249 (MQ=255)
                           cccGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTTCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACTTCTGGCTAACGTTGCg                                                                                                                                                                                                            <  2:297082/250‑1 (MQ=255)
                                  cGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAg                                                                                                                                                                                                     <  1:408000/251‑1 (MQ=255)
                                                  ggCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGCTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCgcagc                                                                                                                                                                                      >  1:511138/1‑250 (MQ=255)
                                                   gcgcGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCAc                                                                                                                                                                                    >  1:51527/1‑251 (MQ=255)
                                                                 gCTCTCTGGGGACCCACCGCATGGGCAACTGGCGAAGGCGTTTAATCAGCTCCAACACCGCCCGGTTCGCGCCAATATTGGTGGCATTAACAAAGGTCGCAATGAATATCATGTCTTTATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCTGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCTCGCAGTACGTTGCAATATCCTg                                                                                                                                                                      <  2:519507/251‑1 (MQ=255)
                                                                 gCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTg                                                                                                                                                                      >  1:590219/1‑251 (MQ=255)
                                                                 gCTCGCTGTGGACCCAACGCATGGGCAACTGGCGCAGGGTTTTAATTATCTACAACATCTCCCGTTTCGCGCCAATATCGGTGGGATTAACAAAGTGCTCAATGAATATCATGTCTTTATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCt                                                                                                                                                                       <  1:628725/250‑1 (MQ=255)
                                                                      ctGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCg                                                                                                                                                                   <  1:136620/249‑1 (MQ=255)
                                                                        ggggACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATcc                                                                                                                                                               <  1:49543/251‑1 (MQ=255)
                                                                        ggggACCCAACGCATGGGCAACTGGCCCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGTTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAATTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGAt                                                                                                                                                                 <  1:117721/249‑1 (MQ=255)
                                                                             cccAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCg                                                                                                                                                                                                            >  1:331181/1‑201 (MQ=255)
                                                                             cccAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGaaa                                                                                                                                                          <  1:242457/251‑1 (MQ=255)
                                                                             cccAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGaa                                                                                                                                                           >  1:217120/1‑250 (MQ=255)
                                                                             cccAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGGCTGGCTAACGTTGCg                                                                                                                                                                                                            <  2:331181/201‑1 (MQ=255)
                                                                                       gggCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGat                                                                                                                                                <  1:223812/251‑1 (MQ=255)
                                                                                        ggCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGatg                                                                                                                                               >  1:529102/1‑251 (MQ=255)
                                                                                                   ccggCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGAGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGTCGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCATGCCGTTCCCGAAAGTTCGGTGATGTATGTTGCCGc                                                                                                                                    >  2:4276/3‑251 (MQ=255)
                                                                                                           tAATCCGCTACAACATATCCCGTTTCGCGCCAATATCGGTGGTATTAACAAAGTTCGCATTGATATTCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATCACGTCTGCCTAACGTTGCGTGGTCAGTATGTTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGCCCGGTGATGTATGTTGCCGCTCCGAtt                                                                                                                             <  2:97660/250‑1 (MQ=255)
                                                                                                           tAATCAGCTACAACATCGCCCGTTTCGCGCCACTATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAACGCAGTTGGGCAGGATTATTCGCGCTGGAACGACGTCTGGCTAACGTTGCGTGGTCAGTTTGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCTGAAAGTTCGGTGATGTATGGTGCCGCGCCggtt                                                                                                                             >  2:270227/1‑247 (MQ=255)
                                                                                                               cAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGa                                                                                                                        <  1:134192/251‑1 (MQ=255)
                                                                                                                     caacaTCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTc                                                                                                                  <  1:304694/251‑1 (MQ=255)
                                                                                                                        caTCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTTCCGCGCCGATTATGGACGGCTcgcg                                                                                                               >  2:260222/1‑251 (MQ=255)
                                                                                                                        caTCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTcgcg                                                                                                               >  1:508434/1‑251 (MQ=255)
                                                                                                                             cgcttttCGCGCCATTTGCGGTTGCATTAACAAAGTGCGCAATGAATATCATGTCAATATGACCGATTCGCAGGTTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGGGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTa                                                                                                          <  2:656271/247‑1 (MQ=255)
                                                                                                                             ccctgttCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAAGAAAGACGTTGGGCAGGATTTTTCGCGCTGGAATGTCGGCTGGTTAACGTTGCGTGGTCCGTATGGTGCGCGCCGCACGTTGGAAAATCCTGCCGATCCCGAAAGTGCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTa                                                                                                          <  1:594315/246‑1 (MQ=255)
                                                                                                                                 tttCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTgg                                                                                                      >  1:335265/1‑251 (MQ=255)
                                                                                                                                      cgcCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGttt                                                                                                 >  1:14871/1‑251 (MQ=255)
                                                                                                                                         cAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTg                                                                                               >  1:229176/1‑250 (MQ=255)
                                                                                                                                             atCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGatg                                                                                                                                               >  1:52283/1‑198 (MQ=255)
                                                                                                                                             atCGGTGGCATTAACAAAGTGCGCAAGGAATATCATGGCTATATGACTGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATGATGCGCGCTGGAATGACGTCTGGCTACCGTTGCGTGGTCAGTAGGGTGCGCGAAGCACGGTGCAAAACCCTGCCGATCCCGAAAGTTCGGTGagg                                                                                                                                               <  2:52283/198‑3 (MQ=255)
                                                                                                                                                 gTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTc                                                                                                                                                      <  1:381595/187‑1 (MQ=255)
                                                                                                                                                 gTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTc                                                                                                                                                      >  2:381595/1‑187 (MQ=255)
                                                                                                                                                           aaaaGTGCGAAAGGAATGTCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCGGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTGGCAAACCCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAAcg                                                                            <  2:590219/250‑1 (MQ=255)
                                                                                                                                                            aaaGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCACACCGAAcgc                                                                           >  2:121264/1‑251 (MQ=255)
                                                                                                                                                                   gctaTGAATATCATGTCTATATGACCGATTCGGAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTGGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATg                                                                    <  2:255465/248‑1 (MQ=255)
                                                                                                                                                                     aaTGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGc                                                                 >  1:600317/1‑251 (MQ=255)
                                                                                                                                                                     aaTGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGc                                                                 >  1:196753/1‑251 (MQ=255)
                                                                                                                                                                            aatcctgtctatctgctcgtttcgcAGGGTAAAGTATTGCTCGATTCGGAAAAGAATGCCGTTGGGAAGGTGTTTTCGCGCTGGAATGATGACTGTCTAACGGTGCGTGGTGAGTAGGGTGCTCGCAGCATGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGGAGGTTGCCGCGCCGATTATGGACGGGGCGCGGCTGATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGTCTCCGGt                                                           <  2:602867/229‑1 (MQ=255)
                                                                                                                                                                               cATGTCTATATGACTGATTCGCAGGTTAAAGTTTTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGTCGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCa                                                         <  2:14871/250‑1 (MQ=255)
                                                                                                                                                                                             cGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTagcgagc                                          >  1:183838/1‑251 (MQ=255)
                                                                                                                                                                                               attcgcAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGCTGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGc                                                                                                     >  2:82391/4‑191 (MQ=255)
                                                                                                                                                                                               attcgcAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGc                                                                                                     <  1:82391/188‑1 (MQ=255)
                                                                                                                                                                                               agtcgcAGGGTAACGTATTGTGCGATGCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCTCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAAGCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGGAGCGAgcgg                                        <  2:126663/248‑1 (MQ=255)
                                                                                                                                                                                                   gcAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAAGCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGaa                                    <  1:620334/251‑1 (MQ=255)
                                                                                                                                                                                                          aaaGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATg                             <  1:580329/251‑1 (MQ=255)
                                                                                                                                                                                                           aaGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATgg                            <  1:260222/251‑1 (MQ=255)
                                                                                                                                                                                                             gTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTATGGCTAACGTTGAGTGGACAGTATGGTGCGAGCAGCACGTAGCAAAATCCTGCCGATCCCGAAAGTGCGGAGATGTATGTTGACGCGCAGATTATGGACGGCTCGCGGCTTATTGGAGTTTTGAGCGTAGGCAAACCGAAAGCGGCGTGGACTCCGGGCCTTAAGCGTCGCGAGAGGCGGATTTTATgggg                          >  2:203695/1‑250 (MQ=255)
                                                                                                                                                                                                              tATTGTTCGATGCGGCAAATAAAACCGTTGGGCAGGATTATTCGCGCTGGACTGACGTCTGGATAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGACAGTTCGGTGATGTATGTTGCCGAGACGATTATGGACGGCCCGCGGCTTAATGGCGTTTTGATCGTAGGCTACACGCACGCGGAGCAGGCTACGGTCATGACGCGTAGCGAGCCGCGAATTTTATTGGcc                         >  2:126654/1‑251 (MQ=255)
                                                                                                                                                                                                               aTTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGGTCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCTTTAAGCGTAGCGAGCGGCGAATTTTAGGGGcc                         >  2:337048/1‑250 (MQ=255)
                                                                                                                                                                                                                 tgttcctttCGGAAAAAAAAGCCGTGGGGCAGGATTATTCGCGCTGGAATAACTTCTGTCTAACGTTGCTTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAg                       <  2:508434/243‑1 (MQ=255)
                                                                                                                                                                                                                     cGATTCGGCAAATAAAGACGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGAAAACTCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCACGGCTTATTGGCGTTTTTAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAgcaccc                  >  2:585194/1‑247 (MQ=255)
                                                                                                                                                                                                                         tCGGCAAATAAAGCCGTTGGGCAGGATTACTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTGGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCAttt               <  1:120465/250‑1 (MQ=255)
                                                                                                                                                                                                                          cGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTttg             >  1:335950/1‑251 (MQ=255)
                                                                                                                                                                                                                          cGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTttg             <  1:229325/251‑1 (MQ=255)
                                                                                                                                                                                                                            gCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTttgct           >  1:310793/1‑249 (MQ=255)
                                                                                                                                                                                                                                 tAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTgggg       >  1:673389/1‑250 (MQ=255)
                                                                                                                                                                                                                                 tAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGAATATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGAAGGCCAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTAGTGGCAAGCGCCATTTTGCTggggg      >  2:68293/1‑250 (MQ=255)
                                                                                                                                                                                                                                    ggCCGTTGGGCAGGATTATTCGGGTTGGAATGAGGTCTGGCTAGCGTTGCGTGGTTAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTg   <  2:183773/250‑1 (MQ=255)
                                                                                                                                                                                                                                    aGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGAATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGAGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTg   >  1:246817/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ccGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGc  >  1:116988/1‑250 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
AGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGC  >  NC_000913/4636801‑4637279

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