Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,725,370 A→C L126* (TTA→TGA)  kgtP ← alpha‑ketoglutarate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,725,3700AC98.4% 255.3 / ‑3.9 63L126* (TTA→TGA) kgtPalpha‑ketoglutarate transporter
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (33/29);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

TAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTG  >  NC_000913/2725222‑2725517
                                                                                                                                                    |                                                                                                                                                   
tAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCaatcat                                                                                                                                                   >  2:822923/1‑151 (MQ=255)
     gTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCgg                                                                                                                                              >  1:415722/1‑151 (MQ=255)
      ttGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCgg                                                                                                                                              >  2:43405/1‑150 (MQ=255)
        gTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAg                                                                                                                                            <  2:840368/150‑1 (MQ=255)
          ccgccgATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGccc                                                                                                                                         >  2:589474/1‑151 (MQ=255)
            gccgATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCAc                                                                                                                                       <  2:533830/151‑1 (MQ=255)
            gccgATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCAc                                                                                                                                       <  1:498040/151‑1 (MQ=255)
             ccgATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACg                                                                                                                                      >  2:435704/1‑151 (MQ=255)
               gATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCAGAGCCCACGt                                                                                                                                     >  2:796664/1‑150 (MQ=255)
                   aaCGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCta                                                                                                                                <  2:323825/151‑1 (MQ=255)
                   aaCGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCt                                                                                                                                 >  2:747624/1‑150 (MQ=255)
                        cacaTACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGtt                                                                                                                           <  1:278498/151‑1 (MQ=255)
                          caTACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTc                                                                                                                         <  1:421192/151‑1 (MQ=255)
                            tACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCAt                                                                                                                       >  1:131855/1‑151 (MQ=255)
                             aCTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCAt                                                                                                                       >  2:243191/1‑150 (MQ=255)
                               tGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAAc                                                                                                                    >  1:306053/1‑151 (MQ=255)
                                      atgCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTgg                                                                                                                >  2:976167/1‑148 (MQ=255)
                                         cGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGaggc                                                                                                           >  1:827247/1‑150 (MQ=255)
                                             taaaCCTGTGCGACCTTCAACCGAAATTTTACTCATACAGGTCGCGCTAGGGCCATATTCTCCGCCAACAGATTATCACTGAAATAAACGAGCGAGCTTATATCACGCCAGTGCCAACGGACGTATAGTTTTATAACCTGGGagt                                                                                                            >  1:968324/2‑144 (MQ=255)
                                              aaaCCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTGTAGTTTCATAACCTGGGAGGCAGGCGa                                                                                                     <  1:613778/151‑1 (MQ=255)
                                                aCCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGAt                                                                                                    >  1:215280/1‑150 (MQ=255)
                                                    ttGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAAc                                                                                                 <  1:163953/149‑1 (MQ=255)
                                                     tGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCa                                                                                               <  2:411805/150‑1 (MQ=255)
                                                       cgcCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCg                                                                                            <  1:698991/151‑1 (MQ=255)
                                                       cgcCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCg                                                                                            >  2:280665/1‑151 (MQ=255)
                                                           cTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATcc                                                                                        >  2:591628/1‑151 (MQ=255)
                                                             tCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATcc                                                                                        <  2:617959/149‑1 (MQ=255)
                                                              cAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGa                                                                                      <  1:879064/150‑1 (MQ=255)
                                                               aaCGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGaa                                                                                     >  2:965326/1‑150 (MQ=255)
                                                                aCGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAac                                                                                   >  2:609272/1‑151 (MQ=255)
                                                                   gCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAAcaca                                                                                <  2:361200/151‑1 (MQ=255)
                                                                   gCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAAcaca                                                                                <  2:306132/151‑1 (MQ=255)
                                                                        ttCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACAcatcat                                                                            >  1:723689/1‑150 (MQ=255)
                                                                           aCTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATacac                                                                        >  1:540730/1‑151 (MQ=255)
                                                                            cTCATATAGGTGGCGCTGGTGCCATATTCTCCGTCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGCACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATac                                                                          >  1:256281/1‑148 (MQ=255)
                                                                             tCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATacacac                                                                      >  2:666548/1‑151 (MQ=255)
                                                                                tataGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCg                                                                    >  1:364765/1‑150 (MQ=255)
                                                                                   aGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATa                                                                 >  1:979820/1‑150 (MQ=255)
                                                                                   aGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATa                                                                 <  2:805153/150‑1 (MQ=255)
                                                                                   aGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATa                                                                 <  1:741512/150‑1 (MQ=255)
                                                                                   aGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATa                                                                 <  2:483377/150‑1 (MQ=255)
                                                                                   aGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAAGCATGCCGGAGCCCACGTCCCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATaa                                                                <  1:700546/151‑1 (MQ=255)
                                                                                            tGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCg                                                        >  2:561747/1‑150 (MQ=255)
                                                                                             ggTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGAt                                                      <  1:867027/151‑1 (MQ=255)
                                                                                               tGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGAttt                                                    >  2:176382/1‑151 (MQ=255)
                                                                                               tGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGAttt                                                    >  1:679047/1‑151 (MQ=255)
                                                                                                 ccATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGAttttt                                                  <  2:220314/151‑1 (MQ=255)
                                                                                                      ttCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCg                                              >  2:15468/1‑150 (MQ=255)
                                                                                                          ccGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGAc                                            <  2:750642/148‑1 (MQ=255)
                                                                                                             ccAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCTGGAGCCCACGTGCCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGtt                                      <  1:956866/151‑1 (MQ=255)
                                                                                                              cAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGttt                                     <  1:405669/151‑1 (MQ=255)
                                                                                                               aaCAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTa                                    <  1:506676/151‑1 (MQ=255)
                                                                                                                aCAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTAt                                   >  2:47788/1‑151 (MQ=255)
                                                                                                                aCAGATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTAt                                   <  1:588439/151‑1 (MQ=255)
                                                                                                                   gATAATCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCgg                                <  1:964914/151‑1 (MQ=255)
                                                                                                                       aTCCCTGAAATAAACGAGCGAGAAGCAATCATGCCAGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTAt                            <  1:796515/151‑1 (MQ=255)
                                                                                                                        tCCCTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATg                           >  1:47411/1‑151 (MQ=255)
                                                                                                                         cccTGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATg                           <  2:201074/150‑1 (MQ=255)
                                                                                                                            tGAAATAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGc                       >  2:621493/1‑151 (MQ=255)
                                                                                                                                 tAAACGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAACCACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAAt                  >  1:793577/1‑151 (MQ=255)
                                                                                                                                  aaaCGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAAt                  <  1:739962/150‑1 (MQ=255)
                                                                                                                                   aaCGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAg                <  1:739797/151‑1 (MQ=255)
                                                                                                                                    aCGAGCGAGAAGCAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGc               >  2:294322/1‑151 (MQ=255)
                                                                                                                                                cAATCATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTAtt   >  1:190628/1‑151 (MQ=255)
                                                                                                                                                 aatcatGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTg  <  1:747481/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TAGCAGTTGTCCGCCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTG  >  NC_000913/2725222‑2725517

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