Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,151 A→G A152A (GCA→GCG creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,637,1510AG96.0% 141.4 / ‑6.4 50A152A (GCA→GCGcreCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (25/23);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCG  >  NC_000913/4636908‑4637389
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
tAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGAtt                                                                                                                                                                                                                                          >  1:273576/1‑250 (MQ=255)
 aaTCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTAt                                                                                                                                                                                                                                        <  1:212171/251‑1 (MQ=255)
                 gCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAGAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCtt                                                                                                                                                                                                                        <  1:474038/251‑1 (MQ=255)
                     gTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTg                                                                                                                                                                                                                    <  2:208109/251‑1 (MQ=255)
                      tttCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTgg                                                                                                                                                                                                                   >  2:19923/1‑251 (MQ=255)
                          gcgcCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGTCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCACAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGtt                                                                                                                                                                                                               >  2:181083/1‑251 (MQ=255)
                                              aaaaaaGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTGGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCg                                                                                                                                                                                                                                             <  2:205565/198‑1 (MQ=255)
                                              aaCAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCg                                                                                                                                                                                                                                             >  1:205565/1‑201 (MQ=255)
                                                      gcgcAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACgcggcg                                                                                                                                                                                    >  1:177890/1‑250 (MQ=255)
                                                                              taccgatttgcaggttaaaGTATTGTGCGATTCGGCAAATAAAGCCGTTGGGCAGGAGTAGTCGCGCTGGAATGACGTCTGTCTAACGTTGCGTGGTCAGTTTGGTGCGCGCAGCACGTTGCAAAAGCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGTCAAACCGACCGCGGCGATGGCTCCGGTCATTAAGCGTagc                                                                                                                                                            <  1:209253/234‑1 (MQ=255)
                                                                                  gattcgcAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGTCGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCAt                                                                                                                                                                      <  1:107272/232‑1 (MQ=255)
                                                                                  gattcgcAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCCCCGGTCAtagatcggagagcgac                                                                                                                                                       >  2:107272/5‑236 (MQ=255)
                                                                                         ggggtaaaggagtggtagatgcgtcAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGGGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAAtt                                                                                                                                                <  2:139283/227‑1 (MQ=255)
                                                                                             taaaGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGGGCGCTCAGAACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGACATTAAGCGTAGCGGGCGGCGAATTTTAt                                                                                                                                            >  2:191535/2‑251 (MQ=255)
                                                                                             taaaGAATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGATGGACTGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTAGGTGATGTATGTTGCCGCGCCGCTTATGGACGGCTCGCGGCTTATTGGAGTTTTGAGCGTAGGCAAACCGAACGAGGCGATGGCTCCGGACATTAAGCGTAGCGAGCGGCGAATTTTAt                                                                                                                                            >  2:341816/2‑251 (MQ=255)
                                                                                                      gTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAgcg                                                                                                                                   >  1:333505/1‑251 (MQ=255)
                                                                                                             tCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGGGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCAttt                                                                                                                             <  1:116036/250‑1 (MQ=255)
                                                                                                             tCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGAAGTCTGTCTAACGTTGCGTGGTCAGTAGGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCAttt                                                                                                                             <  1:8288/250‑1 (MQ=255)
                                                                                                                gCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTttgct                                                                                                                         <  1:209857/251‑3 (MQ=255)
                                                                                                                  aaagaaaGCCGTTTGGCAGGATTATTCGCGCTGTAATGACGTATGGCTAACGTTGCGGGGACAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGCTCCAGAAAGTTCGGTGATGTCTCTTGCCGCGCCGTTTATGGACGGCTCGCGACTTATTAGCGTTTTGATCGTAGTCAAACCGAACGCGGCGATGGCTGCGTTCATTCAGCGTATCGAACGGCGAATTTTATGGGCCAACGCCCtttttaggg                                                                                                                       >  2:264198/5‑246 (MQ=255)
                                                                                                                  aaaTAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGAGATGTATGTTGCCGCGACGATTATGGAAGGCTCGCGGCTTATAGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGAGGGCTCCGGTAATTAAGCGTAGCGAGCGGTAGATTTTATGGGCCAGCGCCTTTTTGCTgg                                                                                                                       >  2:264172/1‑251 (MQ=255)
                                                                                                                    aTAAAGCCGTTGGGCAGGATTATTCGCGCGGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCAGATCCCGAAAGTTAGGTGCTGCATGTTGACGCGCCGATAATGGACGGCCCGCGGCTTATTGGCGTTTTGAGCGTAGACAAACCGAACGCGGCGATGGTACCGGTCATTAAGCGTGGCGCGCGGGGGTTTTTATGGGCCAGCACCATTttgcggcgg                                                                                                                     >  2:346808/1‑245 (MQ=255)
                                                                                                                     tAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCtt                                                                                                                                                                                                                        <  1:359581/151‑1 (MQ=255)
                                                                                                                     tAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCtt                                                                                                                                                                                                                        >  2:359581/1‑151 (MQ=255)
                                                                                                                      aaaGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGGCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCCTTTTGCTGGGGAt                                                                                                                   >  2:289172/1‑251 (MQ=255)
                                                                                                                       aaGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCAtttt                                                                                                                            >  2:277429/1‑241 (MQ=255)
                                                                                                                       aaGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCAtttt                                                                                                                            <  1:277429/241‑1 (MQ=255)
                                                                                                                               gggcAGGATTAGCCGGGCTGGAATAACGTTTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAACCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTTCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTgg                                                                                                           <  1:325083/247‑1 (MQ=255)
                                                                                                                               ggccAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGGGCAGCACTTTGCAAATTCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGt                                                                                                          <  1:447657/248‑1 (MQ=255)
                                                                                                                                ggcAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTAAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGGGATGTATGTTGCCGCGCCGCTTATGGACGTCTCGCGGCTTATTGGCGTTTTGAGCGGTGGCAAACCGGACGCAGCGATGGCTCCGGTCATTAAGCGTAGCGGGGAGCGAATTTTATGGGCAAGCGCACTATTGCTGGGGATTGCACTgggg                                                                                                         >  2:344438/3‑249 (MQ=255)
                                                                                                                                        tattCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGgcg                                                                                                  <  1:66851/250‑1 (MQ=255)
                                                                                                                                                         cGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGCTTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCCGGCATGGTTtggtggg                                                                                >  1:66083/1‑250 (MQ=255)
                                                                                                                                                           tCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATc                                                                              >  2:60250/1‑251 (MQ=255)
                                                                                                                                                                     gtggagtggtgtggtgggggcgcgcAGCACGTTGCAAAACCCGGCCGATCCCGAAAGGTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATAGGCGTAGGCATGGTTTGGTGGATCAACCGctcca                                                                    <  2:79054/232‑3 (MQ=255)
                                                                                                                                                                         cGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGc                                                                <  1:140497/251‑1 (MQ=255)
                                                                                                                                                                            ggTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCg                                                             >  1:362393/1‑251 (MQ=255)
                                                                                                                                                                            ggTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCg                                                             <  1:268461/251‑1 (MQ=255)
                                                                                                                                                                                aGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACAGCTCCATTGCCCGGCTc                                                         >  2:273551/1‑251 (MQ=255)
                                                                                                                                                                                        gcggtcagcaAGTTGAAAAACCCTGCCGATCCCGAAAGTTCGGTGATGTATGGTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCATTCGCTa                                                 <  2:177890/246‑1 (MQ=255)
                                                                                                                                                                                         cgcgcAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTAt                                                >  2:224559/1‑251 (MQ=255)
                                                                                                                                                                                                      gCAAAACCCTGCCGTTCCCGAATGTTCGGTGATGTATGTTGCCGCGCCGAGTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCa                                   <  2:130788/251‑1 (MQ=255)
                                                                                                                                                                                                           aTCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCCTGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGGTTCCGTCACTGAc                              >  1:411405/1‑251 (MQ=255)
                                                                                                                                                                                                             ccTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACa                             >  1:404134/1‑250 (MQ=255)
                                                                                                                                                                                                             ccTGCCGATCCCGAAAGTGCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACaa                            <  2:280453/251‑1 (MQ=255)
                                                                                                                                                                                                                    aTCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGCTTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGccc                     >  2:75796/1‑251 (MQ=255)
                                                                                                                                                                                                                    aTCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTagcg                                                                                                                                                           >  1:454757/1‑117 (MQ=255)
                                                                                                                                                                                                                    aTCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTagcg                                                                                                                                                           <  2:454757/117‑1 (MQ=255)
                                                                                                                                                                                                                                cggtgATGTTTGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTAGTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCctcttcct         <  2:66083/249‑5 (MQ=255)
                                                                                                                                                                                                                                  gtgATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGACTCCGGTCATTAAGCGTAGCGCGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGCTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCctcttcctg        >  1:239357/1‑245 (MQ=255)
                                                                                                                                                                                                                                       gtatgtTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCg  >  2:1657/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
TAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCG  >  NC_000913/4636908‑4637389

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