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
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        
tACGCTGTAACGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGCCCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAAcccc                                                                                                                                                                                                                                            >  2:344576/1‑249 (MQ=255)
        aaCGCTTTCAGCTCGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGCGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGCAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCcac                                                                                                                                                                                                                                    >  2:24065/1‑250 (MQ=255)
                cAGCTCGTTCTGCATTTGCTTACGCAGTTTGCAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGCCGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGa                                                                                                                                                                                                                            <  2:34554/250‑1 (MQ=255)
                     cGTTCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATGTTCGAACACgg                                                                                                                                                                                                                      <  1:374678/251‑1 (MQ=255)
                        tCTGCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGGAATTTCAGCAGCGGGGTTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTc                                                                                                                                                                                                                    <  1:48894/250‑1 (MQ=255)
                           gCATTTGCTTACGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCAt                                                                                                                                                                                                                  <  2:94161/249‑1 (MQ=255)
                                tGCTTCCGCAGTTTGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTgggg                                                                                                                                                                                                            <  1:116193/250‑1 (MQ=255)
                                    tACGCAGTTTGCAATCAAGCGATGAGAGAGACTCATCCAGCAACAACAGACGAGGCTTGTTAAACACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGAGGTTTGCGTTGAGCGAAGGTTTCCAACTTCACCATCCGCAGGGCTTCCATCACGCAGGTCGTCATTTCAGCAGCGGGGGGTTTTGCATGCGCAACACAAAGGACACAGTTTCGATCACGGTCAtggtgtggaca                                                                                                                                                                                                        >  2:393019/1‑241 (MQ=255)
                                            ttGTAATCCAGCGCTTAGAGCGACTCATCCAGCAACACCAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGGTGCTGTTGACCACCAGAGAGCTGATGCGGTTTTCGTTGAGCGAAGGTTTCCAACTGCGCCATCCGCAGGGCTTCCAGCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGGCACATTTTCGAACTCGGTCATGTGGGGGAAAAGTGCGtt                                                                                                                                                                                                <  2:241820/250‑2 (MQ=255)
                                             tGTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGAGTTGCTGTTGACCACCCGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCAGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCCGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGc                                                                                                                                                                                              >  2:319162/1‑251 (MQ=255)
                                              gTAATCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCt                                                                                                                                                                                             <  1:69536/251‑1 (MQ=255)
                                                aaTCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTgg                                                                                                                                                                                                              <  2:47280/232‑1 (MQ=255)
                                                aaTCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTgg                                                                                                                                                                                                              >  1:47280/1‑232 (MQ=255)
                                                 aTCCAGCGCTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTg                                                                                                                                                                                          >  1:122074/1‑251 (MQ=255)
                                                       cgcTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGaaaa                                                                                                                                                                                     >  2:282596/1‑250 (MQ=255)
                                                        gcTGAGAGCGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGAACACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCCCAGGCCACATTTTCGCACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAAc                                                                                                                                                                                    >  1:173557/1‑250 (MQ=255)
                                                                cGACTCATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTcac                                                                                                                                                                           >  1:28334/1‑251 (MQ=255)
                                                                     cATCCAGCAACAACAGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAg                                                                                                                                                                       <  2:89206/250‑1 (MQ=255)
                                                                                 acaGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGCAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTCATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCggag                                                                                                                                                                                               >  2:326501/1‑211 (MQ=255)
                                                                                 acaGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAg                                                                                                                                                                                               <  1:326501/214‑1 (MQ=255)
                                                                                  caGACGAGGCTTGTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCgg                                                                                                                                                         >  1:202426/1‑251 (MQ=255)
                                                                                              gTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGAt                                                                                                                                             <  2:346672/251‑1 (MQ=255)
                                                                                              gTTAACCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGAt                                                                                                                                             <  1:255240/251‑1 (MQ=255)
                                                                                                 aaCCACCGCGCGAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGt                                                                                                                                           <  2:339583/250‑1 (MQ=255)
                                                                                                       cgcgcgAGCAATGGCGCCGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGtt                                                                                                                                    <  2:273120/251‑1 (MQ=255)
                                                                                                          gcgAGCAATGGCGACGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATc                                                                                                                                 <  1:359503/251‑1 (MQ=255)
                                                                                                                       aCGCGTTGCTGTTGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGCTGCGt                                                                                                                    >  2:372445/1‑251 (MQ=255)
                                                                                                                                   tGACCACCAGAGAGCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAAc                                                                                                         >  1:357856/1‑250 (MQ=255)
                                                                                                                                    cgtttcaaaggagCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAg                                                                                                       <  1:24336/241‑1 (MQ=255)
                                                                                                                                              gagCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGAcc                                                                                             <  1:344576/251‑1 (MQ=255)
                                                                                                                                                gCTGATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTg                                                                                           >  1:350789/1‑251 (MQ=255)
                                                                                                                                                gCTGATGCGGTTTGCGTTGAGCGAAGGATTCCACCTGCACCATCCGCAGGGCTTCCCTCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCCACCCAAAGTCCACATATTCGAACACGGTCATGTGGGGGCAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCACCGGACCGGGGGGGATGTCCTCGTTATCCAGCCTGATTCGTCCGGAATCAACAGATTCCAGACCTg                                                                                           >  2:43084/1‑251 (MQ=255)
                                                                                                                                                   gATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCaa                                                                                        <  1:372445/251‑1 (MQ=255)
                                                                                                                                                   gATGCGGTTTGCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCa                                                                                         <  2:350789/250‑1 (MQ=255)
                                                                                                                                                             gCGTTGAGCGAAGGTTTCCAACTGCACCATCCGCAGGTCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGTTTTTTGCATGCGCAACCCAAAGGCCACATTGTCGAACACGGTCATGTGGGGGAAAAGTGCTTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGTGTGATGTCCTCGTTATCCAGCATGATGGGTTCGGAATCAACAGTTTCCAGACCTGCAATCAGGCgaa                                                                               <  1:286176/250‑1 (MQ=255)
                                                                                                                                                                      gAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGgttgt                                                                     >  2:129891/1‑251 (MQ=255)
                                                                                                                                                                       aaGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGCAGCTCTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGACGCGTCCGGAATCAACAGTTTCAAGACCCGCAATCAGGCGAAGAACGgttgtt                                                                    >  1:7555/1‑251 (MQ=255)
                                                                                                                                                                        aGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGGTTTTGCATGCGCAACCCAAAGGCCACATTTGCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGttt                                                                   <  2:122074/251‑1 (MQ=255)
                                                                                                                                                                                    gCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGCACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAg                                                       <  2:118886/251‑1 (MQ=255)
                                                                                                                                                                                    gCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAg                                                       <  1:129891/251‑1 (MQ=255)
                                                                                                                                                                                     cACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGcc                                                         >  1:86896/1‑248 (MQ=255)
                                                                                                                                                                                         aTCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCGTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAGAAggg                                                  <  2:237324/251‑1 (MQ=255)
                                                                                                                                                                                               aGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAGAAGGGCCacgc                                            >  2:295003/1‑248 (MQ=255)
                                                                                                                                                                                                      ccATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAGAAGGGCCAAGCAGCGTg                                      <  1:24065/250‑1 (MQ=255)
                                                                                                                                                                                                       cATCACGCGGGTCGTAATTTCAGCAGCGGGGGTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGGTGATGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACAGTTTCCAGACCTGCAATCAGGCGAAGAACGGTTGTTTTACCGCAGCCAGAAGGGCCAAGCAGCGTGAg                                    <  2:375387/251‑1 (MQ=255)
                                                                                                                                                                                                                            agcagcGGGGGTTTGTGCCTGCTCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAATCAGTGTTCACATAGCGGTGTTCCGCCGGAACGTGGGTGTTGTCCTCGTTATCCAGCATGATGCGTCCGGAATCAACACTTTCCAGGCCTGCATTCAGGCGAAGAACGGTTGTTTTACCGCAGCCAGAAGGGCCAAGCAGCGTGATGATCTCGCCATTGTTGATAGt               <  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: