Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,909,886 G→A L171L (CTC→CTT pstC ← phosphate ABC transporter membrane subunit PstC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,909,8860GA97.8% 147.9 / NA 45L171L (CTC→CTTpstCphosphate ABC transporter membrane subunit PstC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (29/15);  minor base C (1/0);  total (30/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGG  >  NC_000913/3909649‑3910129
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
ccagaagaaaaggcccccccccATGTTTTCACCGAGCGCGCGGCCCACCCCCAGCATGATGCCGCCGATAACCCCATTTTTGGTGAAGGAAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGCCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCCCGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGgataccgatac                                                                                                                                                                                                                                        <  1:305136/231‑1 (MQ=255)
 ccGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGgataccgatacc                                                                                                                                                                                                                                       >  2:76185/1‑251 (MQ=255)
           aGGTCACCGCCATGGTTTCACCGAGCGCGAGGCACAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAAAACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCgggg                                                                                                                                                                                                                             >  1:16220/1‑251 (MQ=255)
               cACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGACACTGCCGCTACCACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAACCAATACCGTAGGCCGCATCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCCCGCCTGCCGCCAGGATACCGATACCAAATGCGGGGCCac                                                                                                                                                                                                                         >  2:215293/1‑250 (MQ=255)
                cccccaatgttttCCCCGGGCGCGCGCCCCAGCCCGAGCAAGATGCCGCAGATAACACCATTTTGGGTGAACGGAAGAACGATACGCCAGAAACCTTCCCAGGGGGTGCAGCCAATCCCGTAGGCCGCCTCTTTCACCAGCCCCGGGTTTTTTTCGAACACATCCCGCATTACCGCCGCAATGTACGGAATAATCATGATGTCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCaga                                                                                                                                                                                                                        <  1:52647/241‑1 (MQ=255)
                               ccGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCCATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGTACAGCGCGCcaac                                                                                                                                                                                                         >  2:95476/1‑251 (MQ=255)
                                     gcgcgGCCCACCCCCAGCATGATGCCGCCGATAACCCCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATCCCGTAGGCCGACTCTTTCATCATCCCCGGGGTTTGTCCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCgg                                                                                                                                                                                                   <  1:76185/251‑1 (MQ=255)
                                             cacccccagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACTATACGCCATATAACTCCCCAGGTGGTGCAGCCAATCCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTGCAAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCg                                                                                                                                                                                           <  2:80091/248‑1 (MQ=255)
                                              agccccagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCg                                                                                                                                                                                           >  1:319918/1‑250 (MQ=255)
                                               gccccagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGAc                                                                                                                                                                                         >  1:12000/1‑251 (MQ=255)
                                                                    aaaccccATTTTTGGTGACGGAAACAACGTTACGCAAGAAATCTCCCCGGGTGGGGCAGCCAATACCTTAGCCCGCCCTTTTCACCATAACCGGGTTTTTTTCGAACACACCGCGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTcc                                                                                                                                                                    <  2:319774/246‑1 (MQ=255)
                                                                       cacCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGAGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATTCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATCTTGCCGACCGGCTCCTg                                                                                                                                                                  >  1:189752/1‑250 (MQ=255)
                                                                           aTTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCAAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGCTACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATAATATTGCCGACCGGCTCCTGAAAGt                                                                                                                                                             >  2:52672/1‑251 (MQ=255)
                                                                           aTTTTTGGTGAACGGAAGAACGAGACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGACGAGGATCACGCATGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATGTTGCCGACCGGCTCCTgca                                                                                                                                                                >  2:89221/1‑246 (MQ=255)
                                                                                 ggggAACTGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCTAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGAACACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGATCCTGCAAGTAAActg                                                                                                                                                       >  2:291508/4‑249 (MQ=255)
                                                                                  gTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGACCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGc                                                                                                                                                      >  1:54778/1‑251 (MQ=255)
                                                                                      aCGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAAc                                                                                                                                                  >  1:365767/1‑251 (MQ=255)
                                                                                                     cGCCAGATAACTTCCAGGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCCCCGGGTTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGAGACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGAt                                                                                                                                    <  1:250992/251‑1 (MQ=255)
                                                                                                              acttcccggggggtcaacccAATCCCGTAGGCCGACTCTTTCATCATCCCCGGGTTTTGTTCGAACACATCACTCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGcc                                                                                                                           <  1:391293/233‑1 (MQ=255)
                                                                                                                       gtggggCAGCCAATCCCGTAGGCCGACCTTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGcc                                                                                                                  <  1:289274/251‑1 (MQ=255)
                                                                                                                        tggtgCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCg                                                                                                                 >  1:268767/1‑251 (MQ=255)
                                                                                                                        gggtgCCGCAAATACCGTAGGCCGCCTCTTTCATCACCGCCGGGGTTTGTTCAAACACATCCCGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCg                                                                                                                 <  1:307451/250‑1 (MQ=255)
                                                                                                                             cgcccacaccagagagccaaCTCTTTCATCCCCACGGGGTTTTGTGCGACCAAATCACGCTTCACCGCCGCAATGTACGGAATAATCATGATGGCGAGGTTCCCGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCTACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAAc                                                                                                            <  2:87959/232‑1 (MQ=255)
                                                                                                                                   aTACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACt                                                                                                      >  1:175496/1‑251 (MQ=255)
                                                                                                                                    tACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACtt                                                                                                     >  2:314431/1‑251 (MQ=255)
                                                                                                                                     aCCGTCGGCCGACTATTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTg                                                                                                    >  1:302893/1‑251 (MQ=255)
                                                                                                                                     aCCGTAGGCCGCCTCTTTCATAATCACCGGGGGTTGTTTGACCGCATAAAGCATTACCGCCGCAATGTACGGAATAATCATGATGGAGAGGCTCACAAATGCCGCAAGGGTACCGATACCAAATGCGGGGCCAGAGACCAGCACGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGTGAAAAGTGGCGCAAAGATAAAAAAGCACAACAGGCCGTAAACGATAcgtg                                                                                                    >  2:255566/1‑248 (MQ=255)
                                                                                                                                            gCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGc                                                                                             >  1:317180/1‑251 (MQ=255)
                                                                                                                                                 ctctTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACAGCCGCAATGTAAGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGAGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGcc                                                                                         >  1:366710/1‑250 (MQ=255)
                                                                                                                                                         tcatcaCCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTc                                                                                 >  1:153084/1‑250 (MQ=255)
                                                                                                                                                         tcatcaCCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATCCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTc                                                                                 >  2:310789/1‑250 (MQ=255)
                                                                                                                                                            ccccgggggTTTTTTCAAACAAATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATg                                                                             <  1:312243/246‑1 (MQ=255)
                                                                                                                                                                          tagaaacactcctgatgtcccgccgCAATGTACGGAATAACCAGGAGGGCGAGGATCACGCCTGCGGCAAGGATACCGATACCAAAGGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATACGACATGATATTGCCGACCGGCTCCGGAAAGTAAACGGCGAGCAGCGGCGCAAAGAGAAACAGGCCCCGCATGGCGTAAACGATACTTGGAATGGCTGCCAACAGCTCAATGGCGATACCCAGCg                                                                <  2:55794/231‑1 (MQ=255)
                                                                                                                                                                                acaTCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCg                                                                          >  1:192345/1‑234 (MQ=255)
                                                                                                                                                                                    cACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTcagc                                                     >  2:313281/1‑251 (MQ=255)
                                                                                                                                                                                      cGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGATCAATGGCGATACCCAGCGGGCGTTTcagcc                                                    >  1:95597/1‑250 (MQ=255)
                                                                                                                                                                                          ttACCGCCGAAATGCACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAAAGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCagccagcc                                                >  2:169003/1‑250 (MQ=255)
                                                                                                                                                                                               gccgcAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCCGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACACAGCGGGCGTTTCAGCCAGCCAGgcgc                                          >  1:368116/1‑251 (MQ=255)
                                                                                                                                                                                                       gTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCtc                                    >  1:231499/1‑249 (MQ=255)
                                                                                                                                                                                                                     tgattgggAGGATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTACCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGAGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCg                    >  2:60049/8‑251 (MQ=255)
                                                                                                                                                                                                                           cggggATCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCg              <  1:53433/248‑1 (MQ=255)
                                                                                                                                                                                                                             gggATCCCGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCACCAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGaa            <  2:231499/250‑1 (MQ=255)
                                                                                                                                                                                                                             gggATCACGCCGGCCGCAAGGAAACCGATACCAAATGCGGGGCCAGAGAACGGCCCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAAGGGCGAGACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCGCAGTCAGGAACAGGGCGATACCGaa            <  1:215293/250‑1 (MQ=255)
                                                                                                                                                                                                                                aTCACGCCTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACt         >  1:192380/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cTGCCGCAAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCggg  >  1:46913/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
ACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGG  >  NC_000913/3909649‑3910129

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