Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,297 T→C V221V (GTA→GTG yaaJ ← putative transporter YaaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,2970TC97.6% 131.7 / NA 41V221V (GTA→GTGyaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (21/19);  minor base A (1/0);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGAACCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGA  >  NC_000913/7055‑7523
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  
ccAATCATCTGGACAACCCCTGGCGCGGCCGGAGGCGGAGGCCAGGCCGCCGCTGCCGCGGCCGCGTTGGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCt                                                                                                                                                                                                                            <  1:116894/251‑1 (MQ=255)
  aaTCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTgg                                                                                                                                                                                                                          >  1:73526/1‑251 (MQ=255)
        cTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGt                                                                                                                                                                                                                     >  1:308985/1‑250 (MQ=255)
                        gcTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCCGCCAAAAGAACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCCACAGGCGGGTCAGGACCCAGATTACCGCCAt                                                                                                                                                                                                    >  2:591218/1‑251 (MQ=255)
                             cGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAACTGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCCGGCTGGTCAGTACCCAGATTATCGCCATCAACg                                                                                                                                                                                               >  2:526743/1‑251 (MQ=255)
                              ggATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACgg                                                                                                                                                                                              >  1:58691/1‑251 (MQ=255)
                                                                    tggcggcgcccccATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGa                                                                                                                                                         <  1:284433/240‑1 (MQ=255)
                                                                       cGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACcgc                                                                                                                                                     >  1:442586/1‑251 (MQ=255)
                                                                       cGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCAAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTGGCCCGATATTCATTACGCAAATAACCAGGCTGGTCAGTACCCAGATTATCGCCATCACCGGGCCAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACcgc                                                                                                                                                     >  2:563648/1‑251 (MQ=255)
                                                                          cgaccccATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGATTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACcgcg                                                                                                                                                    >  2:372169/5‑249 (MQ=255)
                                                                           caccccATCCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGATTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGt                                                                                                                                                 <  1:372169/248‑1 (MQ=255)
                                                                              ccactccccGCCTCATTGAAAAACATACCGCGCTGAAACCCGTTAGTAATCGCCTGGCTTAAGGTAAACCCCGCCGCGCCTCCTGCCGCGTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCAGCAACGGGACAAACCCCTGCATGAGCCGGGCGCCGCCATGAAGACCGCGAGTGAt                                                                                                                                              <  2:73526/245‑1 (MQ=255)
                                                                                     ccGCCTCATTGGAAAAAATACTGCGCTGAAAACCGCTAGCAATCGCCTGGCTTAAGGTATATACCGCAGCGCCGCCAGCCGCTTCCTGCCAGCCAAAAGCACTATCAAAAATAGACCACATGACGTCGGGCAGTTGCCCGATATTCCTTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAACCCCCGGAATGAGCCAGACGACGCCATGGAGACcgcaaggagacacaaac                                                                                                                                       >  2:263093/1‑237 (MQ=255)
                                                                                                 aaaccttACTCCGCTGAAACCCGTTAGTAACCGCCGGGCTTAGGGTATACCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGGCCGCGAGTGATCGCCAGCAGAGCAAAGACa                                                                                                                           <  2:308985/245‑1 (MQ=255)
                                                                                                 aaaaCATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGCCCGCGAGTGATCGCCAGCAGAGCAAAGAc                                                                                                                            >  2:153053/1‑250 (MQ=255)
                                                                                                           ggccggaaacccgtaataaaTCGCCTGGCTTAAGGTAAATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGaataa                                                                                                                 <  2:38219/233‑1 (MQ=255)
                                                                                                              cTGAAAACCGTTAGTAATCGCCTGGCTTATGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGAACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCAGGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCCTGAAGCCCGCGAGTGATCGCCAGCAGCGCCACGACAGCGAGaacaatac                                                                                                              >  2:148519/1‑249 (MQ=255)
                                                                                                                gAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGa                                                                                                                                                   >  2:60764/1‑212 (MQ=255)
                                                                                                                gAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGa                                                                                                                                                   <  1:60764/212‑1 (MQ=255)
                                                                                                                gAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGaataatacc                                                                                                             <  1:526743/250‑4 (MQ=255)
                                                                                                                 aaaaCCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTg                                                                                                           <  1:109751/251‑1 (MQ=255)
                                                                                                                 aaaaCAGTTAGTAACCGCCTGGCTTAAGGAATAACCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTg                                                                                                           <  2:108881/251‑1 (MQ=255)
                                                                                                                      cGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCaaa                                                                                                                               <  1:214604/226‑1 (MQ=255)
                                                                                                                      cGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCAAGCAGAGCaaa                                                                                                                               >  2:214604/1‑226 (MQ=255)
                                                                                                                          agaatacgcttggctaatgttttcccccGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGTAAGTTGCACGATATTCATTGCGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGAGCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGc                                                                                                   <  2:578512/225‑1 (MQ=255)
                                                                                                                              aTCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCggggg                                                                                              >  1:107214/1‑251 (MQ=255)
                                                                                                                                            ggTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGaa                                                                                 >  1:501607/1‑250 (MQ=255)
                                                                                                                                                              gccTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGCGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGgctcgcgca                                                               >  2:417433/1‑243 (MQ=255)
                                                                                                                                                                      ctacataccagccaAAAGCACTCTCAAAAATAGACCAAATGTCGTGGTGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAACGCAAAGACAGCGAGAATAATACCTGTCACAAGCGGGGGAAAATCAAAAGAAACACTCAGGGCGCGCGACACGGCGTTc                                                      >  2:20626/8‑251 (MQ=255)
                                                                                                                                                                                 ccaAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTcc                                           <  1:417433/251‑1 (MQ=255)
                                                                                                                                                                                 ccaAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTcc                                           <  1:57250/251‑1 (MQ=255)
                                                                                                                                                                                    aaaGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTGCGCCTGAATTCCGCt                                        <  1:167891/251‑1 (MQ=255)
                                                                                                                                                                                                        ccAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGccgta                         <  1:496671/246‑4 (MQ=255)
                                                                                                                                                                                                        ccAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTAAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGCGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTAACCAGCGGGGGAAAATCAAAAGAAAAACTAAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGccgtaagaac                    >  2:496671/1‑243 (MQ=255)
                                                                                                                                                                                                              gACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCa              <  1:265600/251‑1 (MQ=255)
                                                                                                                                                                                                               aCGTGGGGAAGTTGACCGATATTCATTACGCAACTCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGAGACGCCATGAAGACCGCGAGTGATCGCAAGCAGAGCAAAGAAAGCGCGAATAATACCTGTCACAAGCGGGGGAAAATCAAACAAAACACTCAGGGCTCGCGCAACGGCGTTAGCCTGAATTCAGCAGCAAATTAttccgaaggcactgagca              >  2:334784/1‑233 (MQ=255)
                                                                                                                                                                                                                    gggAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGa        <  1:153053/251‑1 (MQ=255)
                                                                                                                                                                                                                    gggAAGTTGCCCGAGAGTCAATACGCAAATCCCCAGGCTAGTCAGTACCCAGAATAGCGCCACCAACCGGACAAACCCCTGCATGAGCCGGGCGACGCCAAGAAGACCGCGAGTGATCGCCAGCAAAGCAAAGTCAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACGCACTGATCGCGAAACGGCGTCCGCCTGAAGTACGCTGAAAAGTATGCCGTAGGGGTTGAGCAAAAAGa        <  2:55406/251‑1 (MQ=255)
                                                                                                                                                                                                                         gTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGgct                                                                     >  1:142526/1‑184 (MQ=255)
                                                                                                                                                                                                                          ttGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCg   >  1:450525/1‑250 (MQ=255)
                                                                                                                                                                                                                           tGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTAGCCTGAATTCCGCTGCAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGa  >  1:628505/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  
CCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGAACCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGA  >  NC_000913/7055‑7523

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