Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,227,628 (A)5→4 intergenic (+117/‑103) metH → / → yjbB homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,227,6240A.100.0% 564.3 / NA 126intergenic (+113/‑107)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (69/57);  total (69/57)

Only 100 of 131 total aligned reads displayed.

ACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACATATAACGATAAGGA  >  NC_000913/4227526‑4227718
                                                                                                  |                                                                                              
aCTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaa                                                                                               >  2:2442932/1‑100 (MQ=255)
   tttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaac                                                                                            >  2:2551551/1‑99 (MQ=255)
    ttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaacg                                                                                           <  1:253402/100‑3 (MQ=255)
     tCCTTACAACAGACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaacga                                                                                          <  1:2365227/100‑4 (MQ=255)
       cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAGACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAAt                                                                                       <  2:979596/100‑1 (MQ=255)
       cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAAt                                                                                       <  1:300030/100‑1 (MQ=255)
       cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAAt                                                                                       >  1:1887354/1‑100 (MQ=255)
        ttACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATg                                                                                      <  2:173739/100‑1 (MQ=255)
         tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGt                                                                                     <  2:2198660/100‑1 (MQ=255)
         tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGt                                                                                     <  2:862552/100‑1 (MQ=255)
         tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGt                                                                                     >  2:1647628/1‑100 (MQ=255)
           caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCa                                                                                   >  1:2319921/1‑100 (MQ=255)
           caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAAGCTGGAATAATCTTAGCCGGAAAACGAATGTCa                                                                                   >  1:410045/1‑100 (MQ=255)
            aagagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAACCGAATGTCaa                                                                                  >  2:1675707/6‑100 (MQ=37)
            aagagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCaa                                                                                  >  2:2669907/6‑100 (MQ=37)
            aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCaa                                                                                  <  2:1464378/100‑1 (MQ=255)
             agagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCaaa                                                                                 >  1:1111142/5‑100 (MQ=37)
             acaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCaaa                                                                                 <  2:663928/100‑1 (MQ=255)
               aaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAg                                                                               <  2:505221/100‑1 (MQ=255)
               aaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAg                                                                               <  1:1360557/100‑1 (MQ=255)
                gacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCaaaa                                                                                >  2:2859703/2‑98 (MQ=38)
                gacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAgg                                                                              >  2:402735/2‑100 (MQ=38)
                 acaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGc                                                                             <  1:1467798/100‑1 (MQ=255)
                  caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAAACTTAGCCGGAAAACGAATGTCAAAAGGca                                                                            >  2:1740930/1‑100 (MQ=255)
                    gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGcaca                                                                          >  2:1338820/1‑100 (MQ=255)
                        gcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTaa                                                                      >  1:2617006/1‑100 (MQ=255)
                        gcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTaa                                                                      >  1:1589541/1‑100 (MQ=255)
                        gcTCAATCAGTACCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTaa                                                                      >  2:1518494/1‑100 (MQ=255)
                         cTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTaaa                                                                     >  2:845096/1‑100 (MQ=255)
                           cAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGa                                                                   <  2:2616987/100‑1 (MQ=255)
                           cAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGa                                                                   >  1:1276918/1‑100 (MQ=255)
                           cAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGa                                                                   >  1:1391243/1‑100 (MQ=255)
                           cAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGa                                                                   >  2:2333427/1‑100 (MQ=255)
                           cAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGa                                                                   >  1:1507840/1‑100 (MQ=255)
                            aaTGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAg                                                                    <  2:2760256/98‑1 (MQ=255)
                             aTGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGCAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATg                                                                 <  2:1843458/100‑1 (MQ=255)
                             aTGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATg                                                                 <  1:2157833/100‑1 (MQ=255)
                              tGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGc                                                                >  1:906751/1‑100 (MQ=255)
                               gAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAATCAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGcc                                                               <  2:1844666/100‑1 (MQ=255)
                               gAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGcc                                                               <  1:1062839/100‑1 (MQ=255)
                                    cctttcatctTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAACGCACAGTAAAGATGcc                                                               <  2:1445224/88‑1 (MQ=255)
                                    cccTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCa                                                          >  2:1613901/1‑100 (MQ=255)
                                     ccTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAg                                                         >  2:2012753/1‑100 (MQ=255)
                                          ctctTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCgtgt                                                    >  1:932928/1‑100 (MQ=255)
                                           tctTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGtt                                                   >  2:188829/1‑100 (MQ=255)
                                           tctTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGtt                                                   >  2:1402363/1‑100 (MQ=255)
                                           tctTTATTAATATGAAACACGTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGtt                                                   >  2:1771688/1‑100 (MQ=255)
                                             tttATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTc                                                 >  1:345985/1‑100 (MQ=255)
                                             tttATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTc                                                 >  2:687036/1‑100 (MQ=255)
                                               tattaATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCg                                               >  1:785186/1‑100 (MQ=255)
                                                attaATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCgg                                              <  2:355225/100‑1 (MQ=255)
                                                attaATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCgg                                              <  1:1068685/100‑1 (MQ=255)
                                                  taATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTc                                            <  1:2232063/100‑1 (MQ=255)
                                                  taATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTc                                            >  1:568249/1‑100 (MQ=255)
                                                      atGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCttt                                        >  2:1004748/1‑100 (MQ=255)
                                                      atGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCttt                                        <  1:1647641/100‑1 (MQ=255)
                                                        gAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTgg                                      <  1:2135115/100‑1 (MQ=255)
                                                        gAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTgg                                      <  1:2380506/100‑1 (MQ=255)
                                                         aaaCACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGt                                     >  2:492003/1‑100 (MQ=255)
                                                         aaaCACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGt                                     >  1:1698231/1‑100 (MQ=255)
                                                             acgtgtaCTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAtt                                 <  1:1253443/95‑1 (MQ=255)
                                                             acTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAtt                                 <  1:2724996/100‑1 (MQ=255)
                                                              cTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTc                                >  2:934723/1‑100 (MQ=255)
                                                              cTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTc                                >  1:2712782/1‑100 (MQ=255)
                                                              cTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTc                                <  2:1276913/100‑1 (MQ=255)
                                                                tataCTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCaa                              >  2:526972/1‑100 (MQ=255)
                                                                  taCTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAAtt                            <  1:2859724/100‑1 (MQ=255)
                                                                  taCTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAAtt                            <  1:894498/100‑1 (MQ=255)
                                                                   aCTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTg                           >  2:1032551/1‑100 (MQ=255)
                                                                   aCTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTg                           >  1:729285/1‑100 (MQ=255)
                                                                   aCTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTg                           >  2:346969/1‑100 (MQ=255)
                                                                    cTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGc                          <  1:1675720/100‑1 (MQ=255)
                                                                    cTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGc                          >  2:1868724/1‑100 (MQ=255)
                                                                    cTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGc                          >  2:735102/1‑100 (MQ=255)
                                                                    cTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGc                          >  2:557972/1‑100 (MQ=255)
                                                                     tGGAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCt                         <  2:924124/100‑1 (MQ=255)
                                                                       gAAACAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTc                       >  2:898429/1‑100 (MQ=255)
                                                                          aCAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAgg                    <  1:9288/100‑1 (MQ=255)
                                                                           cAGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGa                   <  2:1425846/100‑1 (MQ=255)
                                                                            aGGCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAAc                 <  2:1466769/100‑1 (MQ=255)
                                                                             ggCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAAc                 >  1:2215669/1‑100 (MQ=255)
                                                                             ggCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAAc                 <  2:2215656/100‑1 (MQ=255)
                                                                             ggCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAAc                 >  1:2338927/1‑100 (MQ=255)
                                                                              gCTGGAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACa                <  1:931936/100‑1 (MQ=255)
                                                                                 ggAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACatat             <  2:2754959/100‑1 (MQ=255)
                                                                                  gAATAATCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACAtata            >  2:2195660/1‑100 (MQ=255)
                                                                                    ataatCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACATATAAc          >  1:2258470/1‑100 (MQ=255)
                                                                                    ataatCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACATATAAc          <  1:1402368/100‑1 (MQ=255)
                                                                                    ataatCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACATATAAc          <  1:188830/100‑1 (MQ=255)
                                                                                    ataatCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACATATAAc          >  1:2259196/1‑100 (MQ=255)
                                                                                    ataatCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACATATAAc          >  1:2797909/1‑100 (MQ=255)
                                                                                    ataatCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACATATAAc          <  1:557974/100‑1 (MQ=255)
                                                                                        tCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACATATAACGATa      <  1:2012767/100‑1 (MQ=255)
                                                                                        tCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACATATAACGATa      <  1:1440701/100‑1 (MQ=255)
                                                                                        tCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACATATAACGATa      <  2:2100790/100‑1 (MQ=255)
                                                                                        tCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACATATAACGATa      <  1:1072099/100‑1 (MQ=255)
                                                                                        tCTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGCCCTTTGGTCATTCTATTGCTTCAGGAACATATAACGAt       >  1:1218399/1‑99 (MQ=255)
                                                                                         cTTAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACATATAACGATaa     >  2:2631670/1‑100 (MQ=255)
                                                                                          ttAGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACATATAACGATAAg    <  1:2551569/100‑1 (MQ=255)
                                                                                            aGCCGGAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACATATAACGATAAGga  <  2:2296754/100‑1 (MQ=255)
                                                                                                  |                                                                                              
ACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAACATATAACGATAAGGA  >  NC_000913/4227526‑4227718

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