Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 324,318 G→A 6.8% A187A (GCG→GCA)  ykgG → DUF162 domain‑containing lactate utilization protein YkgG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913324,3180GA 6.8% 159.4 / 2.4 60A187A (GCG→GCA) ykgGDUF162 domain‑containing lactate utilization protein YkgG
Reads supporting (aligned to +/- strand):  ref base G (23/32);  new base A (2/2);  total (26/34)
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.86e-01

CGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCATGGCCCGGTGAAAGCGGTGTATCTGATTATTGAGGATTGTTGAGAGAGAT  >  NC_000913/324179‑324460
                                                                                                                                           |                                                                                                                                              
cGGGAGGCGTGGGTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAg                                                                                                                                             <  1:298102/143‑1 (MQ=255)
  ggAGGCGTGGTTCTTTTTTCCGCCGCCGCGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCACCTCGCAGAAAAATTGCACCAGAAAGCGCAGGc                                                                                                                                           >  2:258635/1‑143 (MQ=255)
   gAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGACAAATTGCATCAGAAAGCGAAGGcc                                                                                                                                          >  2:304760/1‑143 (MQ=255)
      gCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGt                                                                                                                                       <  1:754577/143‑1 (MQ=255)
      gCGTGGGTCTTTTTTCCGCCGCCGGGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCTGCGTGTAGCGCAGCTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGt                                                                                                                                       <  1:528661/143‑1 (MQ=255)
       cGTGGTTCTTTGTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTg                                                                                                                                      <  1:406388/143‑1 (MQ=255)
         tGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGaa                                                                                                                                    <  2:84225/143‑1 (MQ=255)
          ggTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAAc                                                                                                                                   <  1:85594/143‑1 (MQ=255)
           gTTCTTTTTTCCGCCGCCGAGCGCGGGGGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACg                                                                                                                                  <  1:489603/143‑1 (MQ=255)
             tcctttgtCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACGATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGaa                                                                                                                                <  1:534619/140‑1 (MQ=255)
             tCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCGCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGaa                                                                                                                                <  1:239607/143‑1 (MQ=255)
             tCTTTTTCCCGCCCCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGaa                                                                                                                                <  1:618316/143‑1 (MQ=255)
              ctcttttCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAAt                                                                                                                               <  2:263960/140‑1 (MQ=255)
              cTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAAt                                                                                                                               <  1:356511/143‑1 (MQ=255)
               ttttttCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATg                                                                                                                              <  2:227635/143‑1 (MQ=255)
                  tttCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTACAAGCACTATCCTGCCGCGTGTAGCGCACCTCGCCGAAACATTGCATCAGAAGGCGCAGGCCGGTGAACGCATGCCt                                                                                                                           >  1:74568/1‑143 (MQ=255)
                        ccgccgAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCTCAGGCCGGTGAACGAATGCCTTCCTGc                                                                                                                     >  2:801241/1‑143 (MQ=255)
                        ccgccgAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGc                                                                                                                     <  1:807794/143‑1 (MQ=255)
                         cgccgAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTCGACGCATGCCTTCCTGCa                                                                                                                    >  2:20692/1‑143 (MQ=255)
                          gccgAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCAt                                                                                                                   <  1:674317/143‑1 (MQ=255)
                               gcgcGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAAca                                                                                                              <  1:516745/143‑1 (MQ=255)
                                  cgGGCGTTCACTGCGCCTGCTCCCGGAAGATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAAcatca                                                                                                           <  1:506172/143‑1 (MQ=255)
                                     gCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGAATCGGTGGGCGCAGACCGGTGCGCGCAGGCCCTCCTGCATTAACATCATTa                                                                                                        >  1:181684/1‑143 (MQ=255)
                                           aTTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCACAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGcc                                                                                                  >  1:615451/1‑143 (MQ=255)
                                           aTTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCACAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGcc                                                                                                  >  1:479664/1‑143 (MQ=255)
                                            ttGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCCTTAACATCCTTAGCGGccc                                                                                                 >  2:783093/1‑143 (MQ=255)
                                               aGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGGGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACAGCATTAGCGGCCCCAg                                                                                              >  2:317835/1‑143 (MQ=255)
                                                     cTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGCAAGCGCAGGCCGGTGAACGAATGCGTTCCTGCATTACCATCATTAGCGGCCCCAGCTCAAc                                                                                        >  2:390436/1‑143 (MQ=255)
                                                        ccGGAATATTCTCTTTTTATCCTGCGGAAAAGCACTATCCTGCCGCGTGTAGCGCACCTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCCAtggc                                                                                     <  2:762300/143‑1 (MQ=255)
                                                          ggAATATTCTCTTTTTATCCTGCGTAAAAGCACGATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAAcggcgg                                                                                   <  1:659755/143‑1 (MQ=255)
                                                                 gctctTTTTATCCTGCGTAAAAGCACGATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGa                                                                            <  1:472183/142‑1 (MQ=255)
                                                                   tctTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGc                                                                          >  1:147776/1‑143 (MQ=255)
                                                                           tCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCaaa                                                                  <  1:763402/143‑1 (MQ=255)
                                                                             cTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTCTCAAAgt                                                                >  1:415415/1‑143 (MQ=255)
                                                                              tGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAgtc                                                               >  1:311296/1‑143 (MQ=255)
                                                                              tGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAgtc                                                               <  1:98063/143‑1 (MQ=255)
                                                                                   aaaaGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCgg                                                          >  2:253261/1‑143 (MQ=255)
                                                                                     aaGGACTATCCTGGGGGGTGTAGCGCAACTCGCAGAGAAATTGCATCAGAAGGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAg                                                        <  2:731372/143‑1 (MQ=255)
                                                                                        cACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCACGTTCAACGGCGGATATTGAGCTTATCAAAGTCGACGGAGTTc                                                     >  2:131104/1‑143 (MQ=255)
                                                                                          cTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCAt                                                   <  2:683705/143‑1 (MQ=255)
                                                                                          cTACCCTGCCGCGTGGAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCAt                                                   <  1:93368/143‑1 (MQ=255)
                                                                                                tGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCCTGGCCCg                                             >  1:567199/1‑143 (MQ=255)
                                                                                                tGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAAAGGCGGATATTGAGCTTATCAAAGTCGCCGGAGTTCATGCCCCg                                             >  2:546215/1‑143 (MQ=255)
                                                                                                       tgtAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCATGGCCCGGTGAAAg                                      <  1:638804/143‑1 (MQ=255)
                                                                                                         tAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTGATCAAAGTCGTCAGAGTTCGGGGCGCGGTGGAAACg                                    >  2:690138/1‑143 (MQ=255)
                                                                                                           gcgcAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCATGGCCCGGTGAAAGCGgt                                  >  1:349195/1‑143 (MQ=255)
                                                                                                           gcgcAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAAAATCATTAGCGGACCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCATGGCACGGTGAAAGCGgt                                  >  1:324286/1‑143 (MQ=255)
                                                                                                             gcagcgcGCAGAAAAAGTGCATCAGACAGCGCAGGCCGGTGTACGTATCCCTGCCTGCATTAACATGAATAGCGGCTCCAGTGCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCATGGCCCGGTGAAAGCGgtgt                                <  2:485231/137‑1 (MQ=255)
                                                                                                                acgcGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCATGGCCCGGTGAAAGCGGTGTATc                             <  1:282101/140‑1 (MQ=255)
                                                                                                                   cGCAGAAAGAGTGCATCAGAAAGCGCAGGCCGGTGGACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCATGGCCCGGTGAAAGCGGTGTATCTGa                          <  2:238129/143‑1 (MQ=255)
                                                                                                                       gAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCAGGGCCCGGTGACAGCGGCGTATCCGattat                      >  2:639782/1‑143 (MQ=255)
                                                                                                                        aaaaaTTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCATGGCCCGGTGAAAGCGGTGTATCTGattatt                     <  1:74840/143‑1 (MQ=255)
                                                                                                                         aaaaTTGCATCAGGAAGCGGAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCATGGCCCGGTGAAAGCGGTGTATCTGATTATTg                    <  2:286839/143‑1 (MQ=255)
                                                                                                                             ttGCATCAGAAAGCACAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCATGGCCCGGTGAAAGCGGTGTATCTGATTATTGAGGa                <  2:615451/143‑1 (MQ=255)
                                                                                                                             ttGCATCAGAAAGCACAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCATGGCCCGGTGAAAGCGGTGTATCTGATTATTGAGGa                <  2:479664/143‑1 (MQ=255)
                                                                                                                              tGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCACCGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCATGGCCCGGTGAAAGCGGTGTATCTGATTATCGAGGAt               >  2:540946/1‑143 (MQ=255)
                                                                                                                                cATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTCCATGGCCCGGTGAGAGCGGTGTCTCTGACTCTTGCGGAttg             >  2:380325/1‑143 (MQ=255)
                                                                                                                                 aTCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCAGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGGTCATGGCCCGGTGAAAGCGGTGTATCTGAGTAATGAGGagggg            >  1:469113/1‑139 (MQ=255)
                                                                                                                                  tCAGAGGGCGCCGGCCGGTGAACGAATGCCTTCCTGGATTAACATCAGTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGGTCATGGGCCGGTGAAAGCGGTGTATCTGATTATTGAGGAttgtt           <  2:194902/143‑1 (MQ=255)
                                                                                                                                        cgcgcAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTGATCAAAGTCGTCGGAGTTCATGGCCCGGTGAAAGCGGTGTATCTGATTATTGAGGATTGTTgagaga     <  2:268190/142‑1 (MQ=255)
                                                                                                                                        aGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCATGGCCCGGTGAAAGCGGGGTATCTGATTATTGAGGATTGTTgagaga     >  1:562407/1‑143 (MQ=255)
                                                                                                                                          cgcAGGCCGGTGAACGAATGCCTTCCTGCATTAACGTCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCATGGCCCGGTGAAAGCGGTGTATCTGATTATTGAGGATTGTTgagagaga   <  1:783093/143‑1 (MQ=255)
                                                                                                                                           gcAGGCCGGTGGACGCATGCCTTCCTGCATTAACATCAGTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCATGGCCCGGTGAAAGCGGTGTATCTGATTATTGAGGATTGTTGAGAGAGat  <  1:363317/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
CGGGAGGCGTGGTTCTTTTTTCCGCCGCCGAGCGCGGGCGTTCATTGAGCCTGCTCCCGGAATATTCTCTTTTTATCCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTCGTCGGAGTTCATGGCCCGGTGAAAGCGGTGTATCTGATTATTGAGGATTGTTGAGAGAGAT  >  NC_000913/324179‑324460

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