Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,346,191 A→G V379A (GTA→GCA)  fumB ← anaerobic class I fumarate hydratase (fumarase B)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,346,1910AG98.3% 218.3 / ‑3.7 59V379A (GTA→GCA) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (34/24);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

GGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTC  >  NC_000913/4346046‑4346336
                                                                                                                                                 |                                                                                                                                                 
ggTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCggg                                                                                                                                              >  2:771731/1‑151 (MQ=255)
 gTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGt                                                                                                                                             <  1:764464/151‑1 (MQ=255)
 gTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGt                                                                                                                                             <  1:819464/151‑1 (MQ=255)
   tttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTAt                                                                                                                                           >  1:278373/1‑151 (MQ=255)
     tCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTg                                                                                                                                         >  1:705543/1‑151 (MQ=255)
      cGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGc                                                                                                                                        >  1:527152/1‑151 (MQ=255)
        ccGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCg                                                                                                                                       >  2:668699/1‑150 (MQ=255)
          ggACTCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGACAAACAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGAAAATAGTTCCAGAGAGCTTAAAACGAGTGGATGAAGGATATTGCGaa                                                                                                                                     >  2:211540/1‑150 (MQ=255)
             cccGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCt                                                                                                                                 <  1:312279/151‑1 (MQ=255)
             cccGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCt                                                                                                                                 <  1:999167/151‑1 (MQ=255)
              ccGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTg                                                                                                                                <  2:517667/151‑1 (MQ=255)
                  gtagtaGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGc                                                                                                                             <  2:791375/150‑1 (MQ=255)
                        gATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGAtc                                                                                                                      >  1:647778/1‑151 (MQ=255)
                          tCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGAtctc                                                                                                                    >  2:607365/1‑151 (MQ=255)
                            gggTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCtt                                                                                                                  >  1:119016/1‑151 (MQ=255)
                               tGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTc                                                                                                                >  2:322516/1‑150 (MQ=255)
                                  tCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCg                                                                                                            >  1:627032/1‑151 (MQ=255)
                                   cTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCgg                                                                                                           >  1:835279/1‑151 (MQ=255)
                                   cTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCgg                                                                                                           >  1:391543/1‑151 (MQ=255)
                                     ttGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCggg                                                                                                          <  1:614614/150‑1 (MQ=255)
                                     ttGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCgg                                                                                                           <  2:648676/149‑1 (MQ=255)
                                     ttGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCg                                                                                                            <  2:446677/148‑1 (MQ=255)
                                      tGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCg                                                                                                        >  2:887095/1‑151 (MQ=255)
                                       gATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCgg                                                                                                       >  1:820213/1‑151 (MQ=255)
                                          gTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGtt                                                                                                     <  1:795925/150‑1 (MQ=255)
                                              tGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGCGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAgg                                                                                                 >  2:473182/1‑150 (MQ=255)
                                               gCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTc                                                                                               <  1:912864/151‑1 (MQ=255)
                                                 ggAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCa                                                                                              >  2:101278/1‑150 (MQ=255)
                                                    aGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACtt                                                                                          <  1:629135/151‑1 (MQ=255)
                                                      ttCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTc                                                                                        >  1:8285/1‑151 (MQ=255)
                                                       tCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTc                                                                                        <  2:180742/150‑1 (MQ=255)
                                                           tACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCg                                                                                    <  1:145952/150‑1 (MQ=255)
                                                            aCCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCt                                                                                  >  2:540210/1‑151 (MQ=255)
                                                             ccGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCtt                                                                                 <  2:295182/151‑1 (MQ=255)
                                                                  gTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCt                                                                            <  1:988283/151‑1 (MQ=255)
                                                                    cAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCtt                                                                           <  2:822732/150‑1 (MQ=255)
                                                                    cAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTc                                                                          >  1:830441/1‑151 (MQ=255)
                                                                     aaTCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTc                                                                          >  1:146718/1‑150 (MQ=255)
                                                                        cAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCg                                                                      <  1:101294/151‑1 (MQ=255)
                                                                               ttCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGgc                                                               >  2:941256/1‑151 (MQ=255)
                                                                                tCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGgc                                                               <  2:67680/150‑1 (MQ=255)
                                                                                      tGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGtt                                                         >  2:326267/1‑150 (MQ=255)
                                                                                      tGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTc                                                        >  1:256608/1‑151 (MQ=255)
                                                                                                aaTATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATg                                              <  2:828932/151‑1 (MQ=255)
                                                                                                aaTATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATg                                              <  2:695939/151‑1 (MQ=255)
                                                                                                 atatCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGt                                             >  1:382350/1‑151 (MQ=255)
                                                                                                     ctcGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTg                                         >  2:12207/1‑151 (MQ=255)
                                                                                                         gCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCt                                     <  2:57061/151‑1 (MQ=255)
                                                                                                               gATAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTtg                               >  1:842451/1‑151 (MQ=255)
                                                                                                                 tAATGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTtgt                              >  2:397691/1‑150 (MQ=255)
                                                                                                                    tGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTc                          >  2:378509/1‑151 (MQ=255)
                                                                                                                    tGGTGCCGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTc                          >  2:923595/1‑151 (MQ=255)
                                                                                                                         ccGGTGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGtt                     <  1:194408/151‑1 (MQ=255)
                                                                                                                             tGAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTc                 >  2:796286/1‑151 (MQ=255)
                                                                                                                              gAGCGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCg                <  1:10799/151‑1 (MQ=255)
                                                                                                                                 cGACAAACGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATc             >  1:950669/1‑151 (MQ=255)
                                                                                                                                     aaaCGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATcc            >  2:311937/1‑148 (MQ=255)
                                                                                                                                        cGAGTGGATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGACAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATCCAGATAc      >  1:229087/1‑151 (MQ=255)
                                                                                                                                             ggATGCCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTc  <  2:705463/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTC  >  NC_000913/4346046‑4346336

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑