Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 11,695 G→A A31V (GCC→GTC)  yaaI ← DUF2541 domain‑containing protein YaaI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091311,6950GA96.4% 152.7 / ‑3.3 57A31V (GCC→GTC) yaaIDUF2541 domain‑containing protein YaaI
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (30/24);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

TTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAG‑AAGATTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATATGAATAAAAT  >  NC_000913/11452‑11935
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
tttGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCg                                                                                                                                                                                                                                              >  1:171420/1‑249 (MQ=255)
         cAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCGGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGGGGGGTTTATg                                                                                                                                                                                                                                   >  2:89729/1‑251 (MQ=255)
             gTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCg                                                                                                                                                                                                                               <  1:436019/251‑1 (MQ=255)
               ttATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGACGtt                                                                                                                                                                                                                             >  2:672777/1‑251 (MQ=255)
               ttATTGACGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCGTTCTTTTATTTCTGATGGAATATTCAGACTCTGACTGGCGCTACGGGAGGCTTTGATATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCACCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAAGGAAGCCGAGGAGTTTATGGTCGtt                                                                                                                                                                                                                             >  2:553861/1‑251 (MQ=255)
                       gTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTgcg                                                                                                                                                                                                                     >  2:578647/1‑251 (MQ=255)
                                  tAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGcagcac                                                                                                                                                                                                          >  1:227507/1‑251 (MQ=255)
                                        tGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCAt                                                                                                                                                                                                     <  1:578647/250‑1 (MQ=255)
                                         gATCCAGGCAGGTGTTTGCCCTTCTTTTATTTCTGAAGGTATATTCAGACGCTGACTGGCGCTACGGGCGGCTGTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATc                                                                                                                                                                                                    <  2:187195/250‑1 (MQ=255)
                                            ccAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGt                                                                                                                                                                                                                 >  1:246056/1‑234 (MQ=255)
                                            ccAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGc                                                                                                                                                                                                <  1:94184/251‑1 (MQ=255)
                                            ccAGTCAGTGGTTTGCCCTTCTTGTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGt                                                                                                                                                                                                                 <  2:246056/234‑1 (MQ=255)
                                                 cagtGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGCGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATc                                                                                                                                                                                           >  1:224807/1‑251 (MQ=255)
                                                 cagtGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATc                                                                                                                                                                                           <  1:672777/251‑1 (MQ=255)
                                                   gtgttttGCCCTTCTTTGATTTCTGAAGGAATATGCAGACTCTGACTGGCGCTCCGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTTT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGTGTTTCACAATGCGACAAACAGGATGTTTCAGCGCCAGATTGTTGGTTTCGTTACGGGGCATTGCAATGACGTCGAGGAGGTTATGGTCGTGTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGc                                                                                                                                                                                         <  2:224807/247‑1 (MQ=255)
                                                        ttGCCCTGTTTTTATGTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCCTTGAAATAAACCGATGCACCGCTTAACCGT‑AAATCGCCATGGTCGGCAGAGAGGTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAgg                                                                                                                                                                                     <  2:374904/250‑1 (MQ=255)
                                                                         cTGAAGGAATATTCAGACTCTGACTGGCGCTACGAGCGGCTTTGAAATAAACCGATACACCGCTTAACTGT‑AAATCGCCCTGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGTTCGTTTGCCTGCGCCGTGCAGCAAAGCATCAGGCAAATCGCCAGGCTGGCGGAAATctgaac                                                                                                                                                                    >  2:418587/1‑245 (MQ=255)
                                                                          tGAAGGAATATTCAGACACTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGTTGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACACTGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGACTGGCGGAAATCGTacaa                                                                                                                                                                   >  2:479614/1‑247 (MQ=255)
                                                                            aaGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGGCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACg                                                                                                                                                                <  1:182749/251‑1 (MQ=255)
                                                                             aGGAATATTCAGACTCTGAATGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTCCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGAGGAGATCGTAAAAACgg                                                                                                                                                               >  2:139421/1‑251 (MQ=255)
                                                                                aaTATTCAGACAATGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGATCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACCGCATCAGGCTAATCGCCAGTCTAGCGGGAATctgcaacaaggatt                                                                                                                                                            >  2:256095/1‑238 (MQ=255)
                                                                                  tatTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCACTGACACCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGAGGAAATCGTAAAAACGGATTTc                                                                                                                                                          >  1:472454/1‑251 (MQ=255)
                                                                                   atTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCa                                                                                                                                                         <  2:48544/251‑1 (MQ=255)
                                                                                                 tGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAAGCGCCAGGCGGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCtt                                                                                                                                            >  2:196838/1‑250 (MQ=255)
                                                                                                       tACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCCTGGACGGCAGAGAGTTGTATGCGTTTCACAATGCGCCCAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGAGCCGTGCAGCACAGCATCAGGCTAATCGCCAGTCTGGCGGAAATCGTAAAAACGGATTTAATAAGGATTCTATTAGTggg                                                                                                                                      >  1:379295/1‑250 (MQ=255)
                                                                                                         cGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGAt                                                                                                                                                             <  1:58936/225‑1 (MQ=255)
                                                                                                         cGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGAt                                                                                                                                                             >  2:58936/1‑225 (MQ=255)
                                                                                                           ggcggcTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGttt                                                                                                                                                                                                                                      <  1:593480/150‑1 (MQ=255)
                                                                                                           ggcggcATTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGTGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGttt                                                                                                                                                                                                                                      >  2:593480/1‑150 (MQ=255)
                                                                                                                     aaaTAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAg                                                                                                                             >  1:318193/1‑245 (MQ=255)
                                                                                                                     aaaTAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAg                                                                                                                             <  2:318193/245‑1 (MQ=255)
                                                                                                                       aTAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAGATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATg                                                                                                                      >  2:466302/1‑250 (MQ=255)
                                                                                                                                tGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATc                                                                                                                                                                         <  2:562552/190‑1 (MQ=255)
                                                                                                                                tGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATc                                                                                                                                                                         >  1:562552/1‑190 (MQ=255)
                                                                                                                                 gCACCGCTTAACTGT‑ACATCGCCATGTTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGTAAGAGGTAGAGGGATGAATACCCACt                                                                                                            <  1:553861/250‑1 (MQ=255)
                                                                                                                                 gCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACt                                                                                                            <  2:317993/250‑1 (MQ=255)
                                                                                                                                    ccGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATAGCCAGGCTGGCGGAAATCGAAAAAACGGATTTCATAAGGATTCTCTTAGTAGGAAGAGGTAGAGGGATGAATACCCACTAGtt                                                                                                        >  1:369488/1‑251 (MQ=255)
                                                                                                                                    ccGCTTAAC‑GTGAAATCTCCATGGGCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTGTCAGCGCCAGATCGTTGGGTGCGTTACGCGGCATTGCAATGTCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTTCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAg                                                                                                          <  2:600325/249‑1 (MQ=255)
                                                                                                                                          aaCTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTActgg                                                                                                  <  1:466302/251‑2 (MQ=255)
                                                                                                                                          aaCTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTActgc                                                                                                  >  1:216068/1‑251 (MQ=255)
                                                                                                                                                   aTCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATaa                                                                                            >  2:544046/1‑249 (MQ=255)
                                                                                                                                                                                  ttCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCATTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCCCAGCATCAGGCTAATCGCCAGGCTGGCGGGAATCGTAACAACGAAGTTCATAAGGCTTCTCTTAGTGGGAAGAGGAAGAGGGTTGAATTCCCACTAGTTTCTTGATGATATAGAA‑AAGGTTCAGGTACGGAATCTTTTCttttt                                                          >  2:495844/1‑251 (MQ=255)
                                                                                                                                                                                             aCAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAG‑AAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAAtttt                                               >  2:676361/1‑251 (MQ=255)
                                                                                                                                                                                                 ccAGGAAGGTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGTCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAG‑AAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTtga                                           <  1:43514/250‑1 (MQ=255)
                                                                                                                                                                                                      aaGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAG‑AAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAAt                                      <  1:264661/251‑1 (MQ=255)
                                                                                                                                                                                                          ttCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGGTTTCATAAGGATTCTCTTAGTGGGAAGAGGGCGAGGGATGAATACCCACTAGTTTACTGCTGGTAA‑GAGGAAGGTTCAGGTACGTAATCTTTTTTTTTTATTACACTTTTTTGATGAAGGCCt                                  >  2:473955/1‑251 (MQ=255)
                                                                                                                                                                                                             aGCGCCAGATCGATGGTTTCGTTACGCGGCAGTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTACTCGCCAGGCTGGCGGCAGTCGTAAAAACGGTTTTCATCAAGATTCTCTTAGTGGATAGAGGTAGAGGGTTGATTCCCCACTCGTTTACTGCTGATAAAGAG‑ACGGTTAAGGAACTTACTCTTGTTTTTTTTTTACAATTTTTTGATGAATTCCtttg                               >  2:533352/1‑249 (MQ=255)
                                                                                                                                                                                                              gcgcCAGATCGGTGGTTTCGTTACGCGGCATTGCAATGACGACGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAG‑AAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGc                              <  2:171420/251‑1 (MQ=255)
                                                                                                                                                                                                               cgcCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAG‑AAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCt                             <  1:196838/251‑1 (MQ=255)
                                                                                                                                                                                                                        gTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATaaa                                                                                           >  2:648660/1‑181 (MQ=255)
                                                                                                                                                                                                                        gTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATaaa                                                                                           <  1:648660/181‑1 (MQ=255)
                                                                                                                                                                                                                          tGGTTTCGTTACGCGTCATTGCAATAACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGACAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCCCTAGTTTACTGCTGATAAGGAG‑AAGGATAAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGattcattc                  >  2:424231/1‑251 (MQ=255)
                                                                                                                                                                                                                            gttgcGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATTGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAG‑AAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCtt                <  1:473471/247‑1 (MQ=255)
                                                                                                                                                                                                                                cGTTACGCGGCATTGAAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAG‑AAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCTATTCATTCTTtat             <  1:537259/250‑1 (MQ=255)
                                                                                                                                                                                                                                  ttACGCAGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAG‑AAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTtatttg          >  2:652294/1‑248 (MQ=255)
                                                                                                                                                                                                                                    aCGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGCCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAG‑AAGATTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATATGaa        <  1:89729/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cATTGCAAAGACGCCGAGGAGTTTATGGACGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTACAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAG‑AAGGTTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAAt  >  1:680082/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
TTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGT‑AAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAG‑AAGATTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATATGAATAAAAT  >  NC_000913/11452‑11935

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