Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,555,324 A→G 100% S116P (TCT→CCT)  maeA ← malate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,555,3240AG97.3% 263.7 / ‑2.4 85S116P (TCT→CCT) maeAmalate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (47/34);  total (48/37)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

CTTTAATATTATGGTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAGGGTTTCGTT  >  NC_000913/1555225‑1555420
                                                                                                   |                                                                                                
cTTTAATATTATGGTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACAAGCCGCGTGAACGGCGGTAGATCTCAGa                                                                                                  >  2:2595041/1‑100 (MQ=255)
    aaTATTATGGTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGgaaaa                                                                                              <  1:545824/100‑1 (MQ=255)
     atatTATGGTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAAc                                                                                             >  2:2137488/1‑100 (MQ=255)
     atatTATGGTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCATGAACAGCGGTAGATCTCAGgaaaa                                                                                              >  1:601877/1‑99 (MQ=255)
     atatTATGGTTCGGCACGTTTTGCAAAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAAc                                                                                             <  1:129061/100‑1 (MQ=255)
      tattatGGTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTAAACGGCGGTAGATCTCAGGAAAACg                                                                                            >  1:1750904/1‑100 (MQ=255)
          atGGTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTca                                                                                        >  1:1285842/1‑100 (MQ=255)
             gTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAg                                                                                     >  1:1389010/1‑100 (MQ=255)
             gTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGAACTCAGGAAAACGCTCACAg                                                                                     >  1:2892629/1‑100 (MQ=255)
               tCGGCACGGTTTGCAGTATATCGTCCATATTGTGTGGGTTCTGGTAAGCGATAAACACGCCGCGTGAACGGCGGTAGATTTCAGGAAAACGCTCACAGGc                                                                                   <  2:1775353/100‑1 (MQ=255)
                cGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGACCGGCGGTAGATCTCAGGAAAACGCTCACAGGc                                                                                   >  2:2731585/1‑99 (MQ=255)
                cGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCt                                                                                  >  1:1972506/1‑100 (MQ=255)
                 ggCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTg                                                                                 >  1:816160/1‑100 (MQ=255)
                 ggCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTg                                                                                 >  2:2307067/1‑100 (MQ=255)
                  gCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTgc                                                                                >  1:1271666/1‑100 (MQ=255)
                  gCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTgc                                                                                <  2:2254949/100‑1 (MQ=255)
                   cACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTgcg                                                                               <  2:1972483/100‑1 (MQ=255)
                   cACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTgcg                                                                               <  1:1657191/100‑1 (MQ=255)
                        tttGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATATCAGGAAAACGCTCACAGGCTGCGCCGAc                                                                          >  1:360593/1‑100 (MQ=255)
                             aGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGACGGTAGATCTCAGAAAAACGCTCACAAACTGCGCCGACGGTTg                                                                     >  1:1040341/1‑100 (MQ=255)
                                atatCGTCCATATTGTGCCGTTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTgggg                                                                  <  1:680414/100‑1 (MQ=255)
                                atatCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTTCGCCGACGGTTgggg                                                                  <  2:881566/100‑1 (MQ=255)
                                atatCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTgggg                                                                  <  2:211548/100‑1 (MQ=255)
                                atatCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTgggg                                                                  <  1:680468/100‑1 (MQ=255)
                                atatCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTgggg                                                                  <  1:2572492/100‑1 (MQ=255)
                                  atCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGGTGGGGta                                                                >  2:289290/1‑100 (MQ=255)
                                    cGTCCATATTGTGCCGTTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTGGGGtata                                                              <  2:1767283/100‑1 (MQ=255)
                                    cGTCCATATTGTGCCGGTTCTGGTAAGAGCTAAACCCGCCGCGTGAACGACGGAAGATCTAAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGtata                                                              >  1:1624237/1‑100 (MQ=255)
                                    cGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGtata                                                              <  1:1756680/100‑1 (MQ=255)
                                    cGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGtata                                                              >  1:1211424/1‑100 (MQ=255)
                                     gTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTataa                                                             >  2:1957803/1‑100 (MQ=255)
                                     gTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTataa                                                             >  2:1698397/1‑100 (MQ=255)
                                      tCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTataaa                                                            >  1:2456817/1‑100 (MQ=255)
                                        cATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTataaata                                                          <  2:1152116/100‑1 (MQ=255)
                                         atatTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTataaataa                                                         >  1:1171869/1‑100 (MQ=255)
                                            ttGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAg                                                      <  1:1516237/100‑1 (MQ=255)
                                              gtgCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGc                                                    >  2:2528604/1‑100 (MQ=255)
                                              gtgCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGc                                                    >  1:1051930/1‑100 (MQ=255)
                                              gtgCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGc                                                    >  2:2698276/1‑100 (MQ=255)
                                              gtgCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGc                                                    >  1:335455/1‑100 (MQ=255)
                                               tgCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTGGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGca                                                   <  1:1574340/100‑1 (MQ=255)
                                               tgCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGc                                                    >  2:828218/1‑99 (MQ=255)
                                                gCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTAAGGAAAACGCTCACAGGCTGCGCCAACGGTTGGGGAATAAATAACAgg                                                     >  2:630627/1‑97 (MQ=255)
                                                   ggTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGcatcat                                               >  2:1641605/1‑100 (MQ=255)
                                                    gTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCatcatc                                              <  1:1273487/100‑1 (MQ=255)
                                                         ggTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCAtcatctcatc                                         >  2:2097525/1‑100 (MQ=255)
                                                               agaTAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAgatg                                   <  2:1904474/100‑1 (MQ=255)
                                                                gaTAAACACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAgat                                    <  1:2189046/98‑1 (MQ=255)
                                                                  tAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGAATAAATAACAGGCATCATCTCATCAAgatgat                                 >  1:1545085/1‑99 (MQ=255)
                                                                  tAAACACGCCGCGTGAACGACGTTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCGTCTCATCAAGATGAtt                                >  2:1766296/1‑100 (MQ=255)
                                                                   aaaCACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTg                               <  2:2187420/100‑1 (MQ=255)
                                                                   aaaCACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTg                               <  1:101208/100‑1 (MQ=255)
                                                                   aaaCACGCCGCGTGAACGGCGGTAGAGCTCAGGAAAAAGCTCACAGGTTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTg                               <  2:3035435/100‑1 (MQ=255)
                                                                    aaCACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCAGACGGTTGGGGTATAAATAACAGGCATCATCTCATAAAGATGATTGt                              >  1:455233/1‑100 (MQ=255)
                                                                     acacGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGtt                             <  1:2553875/100‑1 (MQ=255)
                                                                       acGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGTTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTa                           <  2:1847320/100‑1 (MQ=255)
                                                                        cgccgcGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTAc                          <  1:1539609/100‑1 (MQ=255)
                                                                        cgccgcGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGAGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTAAc                          >  1:328175/1‑100 (MQ=255)
                                                                         gccgcGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTAcc                         <  1:818869/100‑1 (MQ=255)
                                                                         gccgcGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTAcc                         <  2:851488/100‑1 (MQ=255)
                                                                         gccgcGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTAcc                         >  1:2107805/1‑100 (MQ=255)
                                                                         gccgcGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTAcc                         >  2:1346691/1‑100 (MQ=255)
                                                                         gccgcGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTAcc                         <  1:274800/100‑1 (MQ=255)
                                                                          ccgcGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCa                        <  1:1252357/100‑1 (MQ=255)
                                                                          ccgcGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCa                        >  2:2986705/1‑100 (MQ=255)
                                                                          ccgcGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCa                        >  2:2731466/1‑100 (MQ=255)
                                                                           cgcgTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCa                        <  2:1749441/99‑1 (MQ=255)
                                                                           cgcgTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAg                       >  2:1228515/1‑100 (MQ=255)
                                                                            gcgTGAACGGCGGTAGATCTCAGTAATACTTTCACAGGCTGCGCCGACGGTTTGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCa                        <  2:399788/98‑1 (MQ=255)
                                                                                gAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGt                  >  2:1922838/1‑100 (MQ=255)
                                                                                  aCGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTag                <  1:2505211/100‑1 (MQ=255)
                                                                                  aCGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTag                >  2:2336282/1‑100 (MQ=255)
                                                                                   cggcggTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTaga               <  2:895951/100‑1 (MQ=255)
                                                                                    ggcggTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTagaa              >  2:521491/1‑100 (MQ=255)
                                                                                      cggTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTagaaga            >  1:356642/1‑100 (MQ=255)
                                                                                         tAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAggg         <  2:610520/100‑1 (MQ=255)
                                                                                         tAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAggg         <  1:1976867/100‑1 (MQ=255)
                                                                                         tAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAggg         <  1:1670330/100‑1 (MQ=255)
                                                                                         naGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAggg         >  1:1805470/2‑100 (MQ=255)
                                                                                           gATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGCTTGTTTACCAGACGGTAGAAGAGGGtt       >  1:482673/1‑100 (MQ=255)
                                                                                             tctcAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAGGGTTTc     >  2:1762108/1‑100 (MQ=255)
                                                                                              ctcAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAGGGTTTCg    >  1:2737867/1‑100 (MQ=255)
                                                                                              ctcAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAGGGTTTCg    >  1:2588919/1‑100 (MQ=255)
                                                                                              ctcAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACAAGACGGTAGAAGAGGATTTCg    >  1:466599/1‑100 (MQ=255)
                                                                                               tcAGGAAAACGCTCACAGTCTGTGTTGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAGGGTTTCGt   <  1:2677555/100‑1 (MQ=255)
                                                                                               tcAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTATTTACCAGACGGTAGAAGAAGGATTCg    >  1:1835138/1‑99 (MQ=255)
                                                                                                cAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAGGGTTTCGtt  >  1:1306448/1‑100 (MQ=255)
                                                                                                cAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAGGGTTTCGtt  <  1:505466/100‑1 (MQ=255)
                                                                                                   |                                                                                                
CTTTAATATTATGGTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAGGGTTTCGTT  >  NC_000913/1555225‑1555420

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