Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,132 A→C G276G (GGT→GGG yaaJ ← putative transporter YaaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,1320AC97.3% 96.2 / ‑5.8 37G276G (GGT→GGGyaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (21/15);  total (22/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

AAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGAACCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGC  >  NC_000913/6907‑7376
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
aaaCTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAAc                                                                                                                                                                                                                              >  2:172526/1‑250 (MQ=255)
                gAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTAACCGAAAGTAATATCAGAATGGCGCTTGCCGTACAGCTGAACAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAAcc                                                                                                                                                                                                             >  1:591513/1‑251 (MQ=255)
                           cACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCg                                                                                                                                                                                                  >  1:83348/1‑251 (MQ=255)
                                     gCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTa                                                                                                                                                                                         >  2:544040/1‑250 (MQ=255)
                                     gCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCAAGTAATATCAACATGGCGCTTGCCGTACAGCTGACCCGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTa                                                                                                                                                                                         >  2:620114/1‑250 (MQ=255)
                                              aTAAGATGAATACCTTCCAGCGGCATGTAGGGTGTGCCGTTACCAGCCAGCAAAATCAGCAGGGCGCTTGCCGTACAGATGACCAGGGAGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTACCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACAAACTGCGCTGAAACCCGTTAGTCATCGCCTGGTTTAAGGTATATcc                                                                                                                                                                               >  1:77371/1‑251 (MQ=255)
                                                aaGCTAAATCCCTACCAGCGGAGTGTGGGTGGCGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCccg                                                                                                                                                                             <  1:542506/251‑1 (MQ=255)
                                                               cAGCGGCATGTAGGTTGTGCCGTTACCCGACAGTAATATCAGCATGGCGATTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCAGCCTCATTGGAAAACCTACTGCGCTGAAAACCGTTCGTAATCGCCTGGCTTAAGGGATATCCCGCCGCGCcgcccgccc                                                                                                                                                              >  2:312202/1‑250 (MQ=255)
                                                                     aaggtgggtgggcccGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCt                                                                                                                                                        <  1:219594/238‑1 (MQ=255)
                                                                      atagtgtgtcacgtgtacgccGCCAGGAATATCAGCATGACGCTTGCCGTTAAGATGACCATGGTGGCGAAAAAAATGCCAAATATCTGGACAACCCCTTGCGCGGCCGTATGCGGAGGCAAGGACGCCGCTGCCTTTGCCGCGTGTGTCTTCTTCCCCTTTCCCGCCTCATTGGAAAACATACTGCGCTGAAAGCCGTTAGTAATCTCCGGGCTTAAGGTATATCCCGCCGCTCCGCCTGCCGCTTACTg                                                                                                                                                       <  2:377631/232‑1 (MQ=255)
                                                                            gtttaccatcacccctccAGTAATATCGGCATGCCGCTGGCCTTGCAGATCACCATGGTTCCGATAAAAGGGCCAATCATCTGGACATTCCCTGGCGCTGCCGGATGCGGAGGCCAGGCCGCCCCTGCCGCTGCCGCGTTTGGCGCCGACCCCAGCCCCGCCTCATTGGAAAACATACTGCTCGGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCGgccagcc                                                                                                                                                 <  2:51876/235‑1 (MQ=255)
                                                                              tgtgCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTgccagcca                                                                                                                                                >  2:169284/1‑250 (MQ=255)
                                                                                 gCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGATAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCCGCCGCTTCCTGCCAGCCaaaac                                                                                                                                            >  2:463049/1‑250 (MQ=255)
                                                                                             cAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCCAAAGCACTCTCaaaa                                                                                                                                 >  2:232223/1‑250 (MQ=255)
                                                                                                        gCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCGCAAAAATAGACCAAATg                                                                                                                     <  1:172526/251‑1 (MQ=255)
                                                                                                           tGGCGCTTGCCGTACAGATGACCAGGGTGTCGACCAAAATGCCAATCATCCGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACg                                                                                                                  >  1:182476/1‑251 (MQ=255)
                                                                                                             gcgcTTGCCGTACAAATGACCAGGGGGTCGATAAAAATGCCAATCTCCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCTTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCGCAAAAATAGGCCAAATGACGTg                                                                                                                <  1:190573/251‑1 (MQ=255)
                                                                                                               gcTTGCCGTACAGCTGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTggg                                                                                                              >  1:195630/1‑251 (MQ=255)
                                                                                                                  tGCCGTCCAGATTACCATGGTGTCGATAAAAATGCCAATCATCGGGACAATCCCTTGCGCTGCGGGTGGCGGAGGCCAGGGCGCCTCTGCCGCTGCCGCGTTTGGCGTGCACCCCATTCCCGCCTCATGGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCGGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTCCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGAGTggggga                                                                                                           <  2:272819/250‑3 (MQ=255)
                                                                                                                   gCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGCCCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGaa                                                                                                           <  1:1393/250‑1 (MQ=255)
                                                                                                                      gTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGGATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCACATGACGAGGGGAGGTTg                                                                                                       >  2:144251/1‑251 (MQ=255)
                                                                                                                      aatcAGATGCCCAGGGTGGCGATAAAAATCCAAAGTATCGGGACAACACCTGGCGCTGCCGGTTGCGGAGGCCAGGACGCCGCGGCCGCTGCCGCGTTTGGCGCGGCCCCCTTTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTg                                                                                                       <  2:67441/248‑1 (MQ=255)
                                                                                                                        aCAGATGACCAGGGTGTCGATAACAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAACACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAGAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGcc                                                                                                     >  2:367133/1‑251 (MQ=255)
                                                                                                                                 cAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTg                                                                                                       >  2:570902/1‑240 (MQ=255)
                                                                                                                                 cAGGGTGTCGATAAAAATGCCAATCACCGGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTg                                                                                                       <  1:570902/240‑1 (MQ=255)
                                                                                                                                           aTAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAAt                                                                                  >  1:200336/1‑251 (MQ=255)
                                                                                                                                              aaaaTCCCAATCACCTGGACACTCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGCCCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGTAAGTTGCCCGATATTCATTACGCAAAtcac                                                                               <  1:55429/251‑4 (MQ=255)
                                                                                                                                               aaaTGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCAcc                                                                              >  1:288932/1‑251 (MQ=255)
                                                                                                                                                agtGCCAATCATCTGGACAATCCCTGGCGCTGCCGGATGCGGAGGCCAGGCCGCCCCTGCCGCGGCCGCGTTGGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCa                                                                             <  2:597039/249‑1 (MQ=255)
                                                                                                                                                   gCCAATCATCGGGACAACCCCTTGCGCGGCCGGTGGCGGAGGCCAGGCCGCCGCTGCCGCTGCCGTGTTTGGCGTCGCCCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGc                                                                          <  2:85980/251‑1 (MQ=255)
                                                                                                                                                     cAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGCACCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAACCCGTTAGTAATAGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTg                                                                        >  1:318265/1‑251 (MQ=255)
                                                                                                                                                            cggcaccagcccgtgcgcgcccGGATGCGGCGGCCAGGACGCCGCTGCCGCTGCGGCGTTTGGCGTCGTCCCCATTCCCTCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTa                                                                 <  1:304487/231‑1 (MQ=255)
                                                                                                                                                                    tCCCTTGCGCTGCCGGATGCGGAGGCAAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCACGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAACAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACAAGGCTGGTCAGTACCCAGAtt                                                         >  2:419416/1‑251 (MQ=255)
                                                                                                                                                                                      gCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGAc                                        >  2:334190/1‑250 (MQ=255)
                                                                                                                                                                                             ccAGGACGCCGCTGCCGCTGCCGCGTTTGGCGCCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCt                                <  1:290581/251‑1 (MQ=255)
                                                                                                                                                                                                  aCGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGCAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCAGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATg                           >  1:385637/1‑251 (MQ=255)
                                                                                                                                                                                                                           cGTCGACCCCATTACCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAACAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCACCGGGACAAACCCCTGCATGAGCCGGGCGACGACATGAAGACcgc  >  2:105825/1‑251 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
AAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGAACCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGC  >  NC_000913/6907‑7376

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