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,2970TC92.7% 115.2 / ‑0.0 41V221V (GTA→GTGyaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (20/18);  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.08e-01

ACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGAACCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCC  >  NC_000913/7067‑7538
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
ccagccccgtGCGCGGCCGGAGGCGGAGCCCGGGCCCCCGCGCCCGCTGCGGCGTTGGGGGCGGACCCCATTCCCGCCTCATGGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCa                                                                                                                                                                                                                                <  2:237657/242‑1 (MQ=255)
 cAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGTCGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGAAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGCCCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAACTCCACAGGCAGGTCAGTAccc                                                                                                                                                                                                                                 >  2:247740/1‑249 (MQ=255)
  aaTCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGCCAGTAccc                                                                                                                                                                                                                                 >  2:360331/1‑248 (MQ=255)
             cgcccGGAGGCGGAGGCCGGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATc                                                                                                                                                                                                                   <  1:360331/248‑1 (MQ=255)
                   gATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCCGTACCCAGATTATCGCCATCAACggg                                                                                                                                                                                                             >  1:93977/1‑251 (MQ=255)
                    aTGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACggg                                                                                                                                                                                                             >  2:93021/1‑250 (MQ=255)
                     gggggagcccggtcgcccgcgcccgcggccgcTTTTGGCGTCGTCCCCTCCCCCGCCTCATTGAAAAACATACTGCGCTGAAAACCGTTAGTAAGCGCCTGGCTAAAGGTATACCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGAc                                                                                                                                                                                                           <  2:51707/229‑1 (MQ=255)
                                    cccccgcccgcggccggGTTGGGCGTCCACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAg                                                                                                                                                                                           <  2:146066/244‑1 (MQ=255)
                                      cgctgcccgctgccGGTTTTGGCGCCCACCCCTTCCCCGCCTCATTGGAAAACATACTGCGCGGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGa                                                                                                                                                                                            <  2:61678/249‑1 (MQ=255)
                                          gccgctgccgcGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGc                                                                                                                                                                                     >  1:331976/1‑251 (MQ=255)
                                          gccgctgccgcGTTGGGCGTCGCCCCCATCCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCggg                                                                                                                                                                                      <  1:245919/250‑1 (MQ=255)
                                              gatgccgcGTTTGGCGACGGCCCCATTCCCGACTCATTGTAAAACATCCTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGGATATACCGCCGCGCCGACTGCCGCTTCCGGCCAGCCAAAAGCACGCTCAAAAATCGCCCCAATGACGTGAGGAAGTTGCCCGATATTAATTACGCAACTCACCCGGCAGGACCGTACCCAGATTACCGCCAACCACGGACAAAACCCCTGCATGAGACGGGCGAc                                                                                                                                                                                  >  2:139540/3‑251 (MQ=255)
                                                      gTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCa                                                                                                                                                                                                                                                >  2:75047/1‑182 (MQ=255)
                                                                             ccttagtGGAAAACATACGGCGGTTATAACCTTTAGATATCGCGTTGCTTAAGGTATCGCCCTCCGCGCCGCCTGCCTCTCCCGGCCAGCCAAAAGCACGCTCAAAAATAGGCCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCTGGCTTGTCAGTACCCAGATTATCGCCATCAACGGGACAATCCCCTGCATGAGCCGGGCGTTGCCATGAAGGCCGCGAGTGTTCGCCAGCaga                                                                                                                                                   <  2:481707/245‑1 (MQ=255)
                                                                                     gAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACAGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCagagccaagcc                                                                                                                                           >  1:30380/1‑244 (MQ=255)
                                                                                           aTACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCg                                                                                                                                       >  1:393739/1‑249 (MQ=255)
                                                                                             aCTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCgaga                                                                                                                                    >  2:66423/1‑250 (MQ=255)
                                                                                             aCTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGaa                                                                                                                                   <  1:394854/251‑1 (MQ=255)
                                                                                                gcgcTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGaat                                                                                                                                  <  1:66423/249‑1 (MQ=255)
                                                                                                     gAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGaataatacc                                                                                                                            <  1:373523/250‑4 (MQ=255)
                                                                                                     gAAAACCGTTAGTAACCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGaataatacc                                                                                                                            <  1:84996/250‑4 (MQ=255)
                                                                                                              tAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACAAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGTTCGCACGCAGAGCAAAGACAGCGAGAATAATACCTGTCCCCAgg                                                                                                                  >  2:142317/1‑250 (MQ=255)
                                                                                                                gTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCgg                                                                                                                <  1:93021/251‑1 (MQ=255)
                                                                                                                      gCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCaga                                                                                                                                                   >  2:140498/1‑210 (MQ=255)
                                                                                                                      gCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCaga                                                                                                                                                   <  1:140498/210‑1 (MQ=255)
                                                                                                                        cTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAAt                                                                                                        >  1:411608/1‑251 (MQ=255)
                                                                                                                                    atatCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAATACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAAc                                                                                            <  1:410861/251‑1 (MQ=255)
                                                                                                                                                  cgccTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGgctcgcgca                                                                              <  1:113950/251‑8 (MQ=255)
                                                                                                                                                                 tccagccaAAAGCACCCTCAAAAATAGACCAAATGACGTGGGGAAGTTACCCGATAATCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATGATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGAGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAAGAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCgact                                                                 <  2:30380/248‑4 (MQ=255)
                                                                                                                                                                 gccagccaAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGcct                                                                 <  1:346452/249‑3 (MQ=255)
                                                                                                                                                                                      aaaTAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGCAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGc                                            >  2:245681/1‑249 (MQ=255)
                                                                                                                                                                                       aaTAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCaaaa                                                                                                   <  1:423176/193‑1 (MQ=255)
                                                                                                                                                                                       aaTAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCAAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCaaaa                                                                                                   >  2:423176/1‑193 (MQ=255)
                                                                                                                                                                                           gACCAAATGACGTGGGGAGGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGCCGCAATGAAGACCGCGAGTGATTGCCAGCAGAGAAAAGACAGCGAGACTAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCACGGGACACGGCGTACGCTTAACCTCAGCTGAAAATTCAGCCATAGGc                                     >  2:15719/1‑251 (MQ=255)
                                                                                                                                                                                           gACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGc                                     >  1:54308/1‑251 (MQ=255)
                                                                                                                                                                                              cAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCCACAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGCGTGATCGCCAGCAGAGCAATGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCACAAGAAAAACTCAGGGCTCGCGCAACGGCGTCCGCCTGCCTTCCGCTGTAAATTATGCCGTAGGCGAt                                  >  2:123262/1‑251 (MQ=255)
                                                                                                                                                                                                       aggggAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCaaaa                          >  2:95815/2‑250 (MQ=255)
                                                                                                                                                                                                                         ttCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCCGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGACAAGAACGcc       >  2:362665/1‑251 (MQ=255)
                                                                                                                                                                                                                                cGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGCACGCCCATcc  >  1:340463/1‑249 (MQ=255)
                                                                                                                                                                                                                                    aaTTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAATAATGACTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAg                                                                                        >  2:186683/1‑159 (MQ=255)
                                                                                                                                                                                                                                    aaTTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAAGAATGACTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAg                                                                                        <  1:186683/159‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
ACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGAACCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCC  >  NC_000913/7067‑7538

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: