Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,342 G→C G206G (GGC→GGG yaaJ ← putative transporter YaaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,3420GC92.3% 110.3 / 1.5 39G206G (GGC→GGGyaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base C (18/18);  total (20/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 = 8.11e-01

GCCGCTGCCGCTGCCGCGTTTGGCGTCGAACCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCC  >  NC_000913/7103‑7581
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
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)
                                         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)
                                                                                                                                                                                                      cAGGATGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa                                 <  2:54308/251‑1 (MQ=255)
                                                                                                                                                                                                          cTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCAGCTGAAAATTATGCCGTAGGCGATGAGCAAAACGACGGCGAACAGAACGCCCATCCAGCGCATCCCCCGccc                             >  2:486054/1‑251 (MQ=255)
                                                                                                                                                                                                                           tATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGTCGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACTGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATg            >  2:347145/1‑251 (MQ=255)
                                                                                                                                                                                                                              cGCCAGCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGTGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCg         <  1:247740/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTc      >  1:119920/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aaCGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGcc  <  1:310811/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
GCCGCTGCCGCTGCCGCGTTTGGCGTCGAACCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCC  >  NC_000913/7103‑7581

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