Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,086,634 G→A A143A (GCC→GCT plaP ← putrescine:H(+) symporter PlaP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,086,6340GA95.2% 200.6 / NA 65A143A (GCC→GCTplaPputrescine:H(+) symporter PlaP
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (28/32);  minor base C (3/0);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGCTCAGAGCCGTAA  >  NC_000913/2086390‑2086872
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
ccGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTc                                                                                                                                                                                                                                          >  1:154522/1‑251 (MQ=255)
 cGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTc                                                                                                                                                                                                                                          >  1:30631/1‑250 (MQ=255)
             gATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGccac                                                                                                                                                                                                                               <  2:419277/249‑1 (MQ=255)
                gCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAgcg                                                                                                                                                                                                                          >  1:188219/1‑251 (MQ=255)
                         gCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACcaca                                                                                                                                                                                                                  >  1:172104/1‑250 (MQ=255)
                               aaaGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGCCGAGTGCTCGCCAACGTACCAGCGCCTTCACCTTCAAATACGCCATAAAAAACCATGCCCAGAATCACCACGATCAGCACTACCTGCAACACTGCGATTACCGTGATGACGTTCGCTACGAATTTCAGCCGACGCAGGTTAAAAGTGGACATAAAGGCCACCAGCGCCACCACAAAcctac                                                                                                                                                                                                           >  1:236036/1‑247 (MQ=255)
                               aaaGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATcc                                                                                                                                                                                                           >  1:454770/1‑251 (MQ=255)
                                     gaAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACAGGCGCATCGCCAGACCAGAATGGAAGAGTGCACGCCAGCGTACCAGCGCCTTAACCTTCAAATACGCCATAAACAACCATGCCCAGAAGCACCGCGAGCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGGCATAAAGGCCACCAGCGCCACCACAAACATCCACgatg                                                                                                                                                                                                     <  2:395088/251‑1 (MQ=255)
                                          acaGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACgatgggatg                                                                                                                                                                                                <  1:356438/251‑1 (MQ=255)
                                           caGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCAGCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACgatgggatgg                                                                                                                                                                                               <  2:188219/251‑1 (MQ=255)
                                              gATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGa                                                                                                                                                                                              >  2:277727/1‑249 (MQ=255)
                                                 tGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGc                                                                                                                                                                                          <  1:47452/250‑1 (MQ=255)
                                                         cGGTAAGCATCTGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCATTCACCTGCCAATACTCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTTCCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGAAGGTGAAAAGCGGTCACAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCCCCAGAGcg                                                                                                                                                                                 <  2:206797/251‑2 (MQ=255)
                                                          gttAATCATCGGGATAACATGCGCAGCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTTTGGAAGTACGCTACGGATTTCATACTACGCATGTTAAAAGCGGTCAGAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCtt                                                                                                                                                                                <  2:567845/249‑1 (MQ=255)
                                                           gTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCAATAAACAACAATGCCCAGAATCACCGCGATCAGCAATACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTCAATGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCCCCAGAGCTTc                                                                                                                                                                               >  2:329501/1‑251 (MQ=255)
                                                             aaTCATAGGGAGAACAGGCGAATCGGCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCTCCGCGATCAGCCCTCCCTGCTACACGACGATTACAGTTTTGGAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCaa                                                                                                                                                                             <  1:267537/251‑1 (MQ=255)
                                                                   cGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAAt                                                                                                                                                                          >  2:412639/1‑248 (MQ=255)
                                                                   cGGGATAACATGCGCATCGCCAGACCAGAACGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGAGGGGATGGAAGGCACCAGATCTTCAAAATAGa                                                                                                                                                                       <  2:50050/251‑1 (MQ=255)
                                                                    gggctcaCATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCCGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAAt                                                                                                                                                                          >  2:390356/7‑247 (MQ=255)
                                                                          aCATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGc                                                                                                                                                                <  1:269745/251‑1 (MQ=255)
                                                                           gaTGCGCATCGCCATACCAGAAATGACGAATGCTCTCCAGCGTACCAGCGCCTTCCCCTTCAAATACTCCATAAAAAACCAGGCCCAGAATCACCTCGACCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTCCGGATTTCAGACTACGCAGGTTGAGGTCGGTCATACAGGCCACCAGCGCCACCACAAACATCCCCGCTGGGATGGGAGGCACCCGAGCTACAAAATAGACTTTCGcc                                                                                                                                                               >  2:517753/2‑251 (MQ=255)
                                                                            aTGCGCCTCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTGCCCCTTCAAATACGCCAGAAACAACCAGGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTTCGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCTCCAGATCTTCAAAATAGATTTTCGCCa                                                                                                                                                              <  2:183755/251‑1 (MQ=255)
                                                                              gcgcATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATATATTTGCGCCAGt                                                                                                                                                            <  2:284449/251‑1 (MQ=255)
                                                                              gcgcATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGt                                                                                                                                                            >  1:194081/1‑251 (MQ=255)
                                                                               cgcATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTa                                                                                                                                                           <  1:446282/251‑1 (MQ=255)
                                                                                gcATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAg                                                                                                                                                          <  1:9681/251‑1 (MQ=255)
                                                                                    cGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATg                                                                                                                                                      <  1:136852/251‑1 (MQ=255)
                                                                                      ccagaccagaATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGtt                                                                                                                                                    <  2:299490/251‑1 (MQ=255)
                                                                                       cagaccagaATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCAAAAACAACAAAGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAACGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTg                                                                                                                                                   >  1:75698/1‑251 (MQ=255)
                                                                                        agaccagaATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGa                                                                                                                                                  <  1:412639/251‑1 (MQ=255)
                                                                                           ccagaATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCGTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGTAAGGCACCAGAGCTTCAAAATAGAGTTGCGCCAGTAGAATTTTGatca                                                                                                                                               <  2:373668/251‑1 (MQ=255)
                                                                                                      gAGTGCTCGCCAGCGTCCCAGCCCGATCACCTTCAAATACGCCATAACCAACCATGCCCACGATCACCGCGATCAGCAACACCCGCAACACGACGATTACAGTGTGGAAGCTCGCTACGGATTACAGACTACGCAGGTTAAAAGCGGTCAGAAAGTCCACCAGCGCCACCCAAAACATCCAGGATGGGATGGAAGGCACCAGAGCTTCAAAAGAGATTTTCGCAAGGAGAATGTTGATCATCGGCGCGAAc                                                                                                                                    <  2:289300/251‑1 (MQ=255)
                                                                                                              gCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTc                                                                                                                            >  1:132133/1‑251 (MQ=255)
                                                                                                                cagagaaCCAGCGCCTTCACCCTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCATCACTCCCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACTTCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTCGCGCCAGTAGAATGTTGATCATCGTCGCGAACAGATAGTCGa                                                                                                                          <  2:153454/245‑1 (MQ=255)
                                                                                                                cAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGCGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGCTGGGATGGAAGGCACAAGAGCTTCAAACTAGGTTTTCGCCAGTAGAATGTTGATCATCGTCGCGAACAGAAAGTCg                                                                                                                           >  2:45574/1‑250 (MQ=255)
                                                                                                                              cttcaccttcaAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCa                                                                                                             >  2:217854/1‑250 (MQ=255)
                                                                                                                                   ccttcaAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCa                                                                                                                       <  1:511236/235‑1 (MQ=255)
                                                                                                                                   ccttcaAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCa                                                                                                                       >  2:511236/1‑235 (MQ=255)
                                                                                                                                      tcaAATACGCCATAACCCACCATGCCCAGAATCACCGCGATCAGCACTAACTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGCTTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCAt                                                                                                    >  1:201831/1‑251 (MQ=255)
                                                                                                                                                cATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGATACGGATTTCAGACTACGCAGGTTCAAAGCGGTCATAAAGGCCACCAGCGGCACCACAAAAATCCACGCTGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGCAGACCAACCCACCATAAAGCCGAc                                                                                           >  1:585067/1‑250 (MQ=255)
                                                                                                                                                   aaaCAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTc                                                                                       <  1:507221/251‑1 (MQ=255)
                                                                                                                                                        acCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCt                                                                                  <  1:37995/251‑1 (MQ=255)
                                                                                                                                                                 aGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTGCGCTACGGATGTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGAttt                                                                         <  1:79704/251‑1 (MQ=255)
                                                                                                                                                                               cAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAgtgt                                                           <  1:217854/251‑1 (MQ=255)
                                                                                                                                                                                   ccTACCTGCAACACGTAGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCCCCACAAACATCCACGATGGGATGGAAGGCTCCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCTCCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGc                                                       <  1:498327/250‑1 (MQ=255)
                                                                                                                                                                                   aCTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGc                                                       >  1:239194/1‑251 (MQ=255)
                                                                                                                                                                                      aCCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCaga                                                    <  1:277727/251‑1 (MQ=255)
                                                                                                                                                                                       ccTGCAACAAGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGGTAAAAGCGGTCATAAAGGCCACAAGCGCCACCACAAACATCCACGATGGGATGGACGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAAAAGATAGTCTAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGGAAGTGTATGCag                                                     >  2:209768/1‑249 (MQ=255)
                                                                                                                                                                                        cTGCAACACGACGATTACAGTGTTGAAGTTCAATACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACAAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGc                                                  >  1:515426/1‑251 (MQ=255)
                                                                                                                                                                                                cgacgaTTACAGTGTTTAAGTGCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCTCCAGCGCCACCACAAACATCCACGAGGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATTTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAATACCAACCTCCCATAAAGCAGACAGTCGGTTTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGaa                                          <  1:192491/251‑1 (MQ=255)
                                                                                                                                                                                                cgacgaTTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGAc                                                                                           <  1:474959/202‑1 (MQ=255)
                                                                                                                                                                                                cgacgaTTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGAc                                                                                           >  2:474959/1‑202 (MQ=255)
                                                                                                                                                                                                cgacgaTTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGCGTATGCAGAGCCAGCAGa                                           >  2:350363/1‑250 (MQ=255)
                                                                                                                                                                                                  acgaTTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAgg                                        >  2:398384/1‑251 (MQ=255)
                                                                                                                                                                                                       tACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGATCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATATTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTa                                                                                 <  1:554249/205‑1 (MQ=255)
                                                                                                                                                                                                       tACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAgc                                   >  1:107033/1‑251 (MQ=255)
                                                                                                                                                                                                       tACAGAGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGACACCACAAACATCCACGATGGGATGGACGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTa                                                                                 >  2:554249/1‑205 (MQ=255)
                                                                                                                                                                                                         cAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAACGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAgcgg                                 >  1:224078/1‑251 (MQ=255)
                                                                                                                                                                                                                            ggATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGc              >  1:224735/1‑251 (MQ=255)
                                                                                                                                                                                                                            ggATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGc              >  1:182302/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaaga                                                                                                                >  1:87680/1‑143 (MQ=255)
                                                                                                                                                                                                                                      cTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaaga                                                                                                                <  2:87680/143‑1 (MQ=255)
                                                                                                                                                                                                                                       tACGCAGGTTAAACGCGGTCATAAAGGCCACCAGCGCCACCACACACATCCAAGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGACTGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGCTCAGAGCCGt    >  2:269835/1‑250 (MQ=255)
                                                                                                                                                                                                                                       tACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGCTCAGAGCCGt    >  1:280851/1‑250 (MQ=255)
                                                                                                                                                                                                                                        aCGCAGGTTAAAAGCGGTCATAAAGGCCCCCAGCGCCACCACAAACATCCACGATGGGATGGAAGTCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGGCGAGCAGAGAAGACCAACCCACCATAAAGCCTACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGCTCAGAGCCGTaa  <  1:371542/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
CCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGCTCAGAGCCGTAA  >  NC_000913/2086390‑2086872

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