Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,185,245 (T)6→5 coding (350/1137 nt) potA ← spermidine preferential ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,185,2450T.91.7% 178.2 / 13.2 48coding (350/1137 nt)potAspermidine preferential ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base . (19/25);  total (20/28)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

TACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAGAAGGGCCAAGCAGCGTGAGGAACTCGCCATTGTTGATAGTCAGATCCAGCTGG  >  NC_000913/1185009‑1185493
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        
tACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGCCCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAAcccc                                                                                                                                                                                                                                            >  2:344576/1‑249 (MQ=255)
        aaCGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGCGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGCAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCcac                                                                                                                                                                                                                                    >  2:24065/1‑250 (MQ=255)
                cAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGCCGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGa                                                                                                                                                                                                                            <  2:34554/250‑1 (MQ=255)
                     cGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATGTTCGAACACgg                                                                                                                                                                                                                      <  1:374678/251‑1 (MQ=255)
                        tCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGGAATTTCAGCAGCGGGG‑TTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTc                                                                                                                                                                                                                    <  1:48894/250‑1 (MQ=255)
                           gCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCAt                                                                                                                                                                                                                  <  2:94161/249‑1 (MQ=255)
                                tGCTTCCGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTgggg                                                                                                                                                                                                            <  1:116193/250‑1 (MQ=255)
                                    tACGCAGTTTGCAATCAAGCGATGAGAGAGACTCATCCAGCAACAACAGACGAGGCTTGTTAAACACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGAGGTTTGCGTTGAGCGAAGGTTTCCAACTTCACCATCCGCAGGGCTTCCATCACGCAGGTCGTCATTTCAGCAGCGGGGGGTTTT‑GCATGCGCAACACAAAGGACACAGTTTCGATCACGGTCAtggtgtggaca                                                                                                                                                                                                        >  2:393019/1‑241 (MQ=255)
                                            ttGTAATCCAGCGCTTAGAGCGACTCATCCAGCAACACCAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGGTGCTGTTGACCACCAGAGAGCTGATGCGGTTTTCGTTGAGCGAAGGTTTCCAACTGCGCCATCCGCAGGGCTTCCAGCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGGCACATTTTCGAACTCGGTCATGTGGGGGAAAAGTGCGtt                                                                                                                                                                                                <  2:241820/250‑2 (MQ=255)
                                             tGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGAGTTGCTGTTGACCACCCGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCAGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCCGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGc                                                                                                                                                                                              >  2:319162/1‑251 (MQ=255)
                                              gTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCt                                                                                                                                                                                             <  1:69536/251‑1 (MQ=255)
                                                aaTCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTgg                                                                                                                                                                                                              <  2:47280/232‑1 (MQ=255)
                                                aaTCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTgg                                                                                                                                                                                                              >  1:47280/1‑232 (MQ=255)
                                                 aTCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTg                                                                                                                                                                                          >  1:122074/1‑251 (MQ=255)
                                                       cgcTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGaaaa                                                                                                                                                                                     >  2:282596/1‑250 (MQ=255)
                                                        gcTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGAACACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCCCAGGCCACATTTTCGCACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAAc                                                                                                                                                                                    >  1:173557/1‑250 (MQ=255)
                                                                cGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTcac                                                                                                                                                                           >  1:28334/1‑251 (MQ=255)
                                                                     cATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAg                                                                                                                                                                       <  2:89206/250‑1 (MQ=255)
                                                                                 acaGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGCAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTCATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCggag                                                                                                                                                                                               >  2:326501/1‑211 (MQ=255)
                                                                                 acaGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAg                                                                                                                                                                                               <  1:326501/214‑1 (MQ=255)
                                                                                  caGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCgg                                                                                                                                                         >  1:202426/1‑251 (MQ=255)
                                                                                              gTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGAt                                                                                                                                             <  2:346672/251‑1 (MQ=255)
                                                                                              gTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGAt                                                                                                                                             <  1:255240/251‑1 (MQ=255)
                                                                                                 aaCCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGt                                                                                                                                           <  2:339583/250‑1 (MQ=255)
                                                                                                       cgcgcgAGCAATGGCGCCGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGtt                                                                                                                                    <  2:273120/251‑1 (MQ=255)
                                                                                                          gcgAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATc                                                                                                                                 <  1:359503/251‑1 (MQ=255)
                                                                                                                       aCGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGCTGCGt                                                                                                                    >  2:372445/1‑251 (MQ=255)
                                                                                                                                   tGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAAc                                                                                                         >  1:357856/1‑250 (MQ=255)
                                                                                                                                    cgtttcaaaggagCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAg                                                                                                       <  1:24336/241‑1 (MQ=255)
                                                                                                                                              gagCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGAcc                                                                                             <  1:344576/251‑1 (MQ=255)
                                                                                                                                                gCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTg                                                                                           >  1:350789/1‑251 (MQ=255)
                                                                                                                                                gCTGATGCGGTTTGCGTTGAGCGAAGGATTCCACCTGCACCATCCGCAGGGCTTCCCTCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCCACCCAAAGTCCACATATTCGAACACGGTCATGTGGGGGCAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCACCGGACCGGGGGGGATGTCCTCGTTATCCAGCCTGATTCGTCCGGAATCAACAGATTCCAGACCTg                                                                                           >  2:43084/1‑251 (MQ=255)
                                                                                                                                                   gATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCaa                                                                                        <  1:372445/251‑1 (MQ=255)
                                                                                                                                                   gATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCa                                                                                         <  2:350789/250‑1 (MQ=255)
                                                                                                                                                             gCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGTCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGG‑TTTTTTGCATGCGCAACCCAAAGGCCACATTGTCGAACACGGTCATGTGGGGGAAAAGTGCTTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGTGTGATGTCCTCGTTATCCAGCATGATGGGTTCGGAATCAACAGTTTCCAGACCTGCAATCAGGCgaa                                                                               <  1:286176/250‑1 (MQ=255)
                                                                                                                                                                      gAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGgttgt                                                                     >  2:129891/1‑251 (MQ=255)
                                                                                                                                                                       aaGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGCAGCTCTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGACGCGTCCGGAATCAACAGTTTCAAGACCCGCAATCAGGCGAAGAACGgttgtt                                                                    >  1:7555/1‑251 (MQ=255)
                                                                                                                                                                        aGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGGTTTT‑GCATGCGCAACCCAAAGGCCACATTTGCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGttt                                                                   <  2:122074/251‑1 (MQ=255)
                                                                                                                                                                                    gCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGCACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAg                                                       <  2:118886/251‑1 (MQ=255)
                                                                                                                                                                                    gCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAg                                                       <  1:129891/251‑1 (MQ=255)
                                                                                                                                                                                     cACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGcc                                                         >  1:86896/1‑248 (MQ=255)
                                                                                                                                                                                         aTCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCGTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAGAAggg                                                  <  2:237324/251‑1 (MQ=255)
                                                                                                                                                                                               aGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAGAAGGGCCacgc                                            >  2:295003/1‑248 (MQ=255)
                                                                                                                                                                                                      ccATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAGAAGGGCCAAGCAGCGTg                                      <  1:24065/250‑1 (MQ=255)
                                                                                                                                                                                                       cATCACGCGGGTCGTAATTTCAGCAGCGGGGGTTTTT‑GCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAGAAGGGCCAAGCAGCGTGAg                                    <  2:375387/251‑1 (MQ=255)
                                                                                                                                                                                                                            agcagcGGGGGTTTGT‑GCCTGCTCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAATCAGTGTTCACATAGCGGTGTTCCGCCGGAACGTGGGTGTTGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACACTTTCCAGGCCTGCATTCAGGCGAAGAACGGTTGTTTTACCGCAGCCAGAAGGGCCAAGCAGCGTGATGATCTCGCCATTGTTGATAGt               <  1:294322/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cggggGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAGAAGGGCCAAGCAGCGTGAGGAACTCGCCATTGTTGATAGTc              >  1:60697/2‑247 (MQ=255)
                                                                                                                                                                                                                                   ggggGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAGAAGGGCCAAGCAGCGTGAGGAACTCGCCATTGTTGGTAGTCAGATc         <  1:192201/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ggtttttGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAGAAGGGCCAAGCAGCGTGAGGAACTCGCCATTGTTGATAGTCAGATcc        >  2:222242/3‑249 (MQ=255)
                                                                                                                                                                                                                                           ttGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAGAAGTGCCAAGCAGCGTTAGGAACTCGCCATTGTTTATATTCAGATCCAGCTgg  <  2:225238/250‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        
TACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAGAAGGGCCAAGCAGCGTGAGGAACTCGCCATTGTTGATAGTCAGATCCAGCTGG  >  NC_000913/1185009‑1185493

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑