Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,173,076 C→A 100% pseudogene (287/442 nt) gatC ← pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,0760CA98.8% 279.9 / ‑3.9 88pseudogene (287/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (39/47);  total (39/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

ATGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCA  >  NC_000913/2172979‑2173172
                                                                                                 |                                                                                                
aTGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGc                                                                                                <  1:1078333/100‑1 (MQ=255)
aTGAAAATACCGGTCAGATAAATTGCACCGATAAAAATGAAACCGGGAATAGTTTGCGGGGAGCACACCTGAATCAGTAACCAGGTAATGGGAGAACAGc                                                                                                >  1:208370/1‑100 (MQ=255)
 tGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGcc                                                                                               <  1:1234202/100‑1 (MQ=255)
 tGAAAATACCGGTCAGATAAATTGCACCGAAAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGAGAGAACAGcc                                                                                               >  2:1686292/1‑100 (MQ=255)
   aaaaTACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCt                                                                                             <  1:2482885/100‑1 (MQ=255)
   aaaaTACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCt                                                                                             >  1:2308793/1‑100 (MQ=255)
    aaaTACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTg                                                                                            <  1:1137233/100‑1 (MQ=255)
    aaaTACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTg                                                                                            >  1:1527782/1‑100 (MQ=255)
       tACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATc                                                                                         <  1:1992444/100‑1 (MQ=255)
        aCCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATcc                                                                                        <  1:306737/100‑1 (MQ=255)
         ccGGTCAGATAAATTGCACCGATAATAATAAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCCGTAACCAGGTAATGGGAGAACAGCCCTGATCCa                                                                                       <  2:1537923/100‑1 (MQ=255)
          cGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCAt                                                                                      >  1:1252660/1‑100 (MQ=255)
              cAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGaa                                                                                  >  2:186946/1‑100 (MQ=255)
              cAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGaa                                                                                  >  1:1362222/1‑100 (MQ=255)
              cAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGAGAGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGaa                                                                                  >  2:1260633/1‑100 (MQ=255)
               aGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAg                                                                                 <  1:1866760/100‑1 (MQ=255)
               aGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAg                                                                                 >  1:281577/1‑100 (MQ=255)
               aGATAAATTCCACCGATAATAATAAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAg                                                                                 <  2:1906216/100‑1 (MQ=255)
                gATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGc                                                                                >  2:957954/1‑100 (MQ=255)
                gATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGc                                                                                >  1:2184720/1‑100 (MQ=255)
                gATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGc                                                                                >  1:2095066/1‑100 (MQ=255)
                    aaTTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGccacc                                                                            <  2:1899889/100‑1 (MQ=255)
                     aTTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGccacca                                                                           >  1:1969935/1‑100 (MQ=255)
                       tGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATa                                                                         <  1:2444784/100‑1 (MQ=255)
                        gCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATAc                                                                        >  2:2050332/1‑100 (MQ=255)
                        gCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATAc                                                                        >  1:1142894/1‑100 (MQ=255)
                        gCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATAc                                                                        >  1:2410083/1‑100 (MQ=255)
                        gCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTAAAGCCACCATAc                                                                        >  2:1992516/1‑100 (MQ=255)
                         cACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATAcc                                                                       >  1:2032868/1‑100 (MQ=255)
                           ccGATAATAATGAAACCGGGAATATTTTGCGGGGAGATAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACGATAcccc                                                                     >  2:1374699/1‑100 (MQ=255)
                           ccGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATAcccc                                                                     <  2:707288/100‑1 (MQ=255)
                             gATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCg                                                                   >  1:2184734/1‑100 (MQ=255)
                             gATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCACCg                                                                   >  2:1256827/1‑100 (MQ=255)
                               taataaTGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGc                                                                 <  1:357604/100‑1 (MQ=255)
                                 ataatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCtt                                                               <  1:2276073/100‑1 (MQ=255)
                                  taatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCttt                                                              >  2:1619002/1‑100 (MQ=255)
                                   aatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCtttt                                                             <  1:1176717/100‑1 (MQ=255)
                                   aatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCtttt                                                             >  2:978799/1‑100 (MQ=255)
                                   aatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCtttt                                                             <  1:2027305/100‑1 (MQ=255)
                                   aatGAAACCGGGAATATTTTGCGCGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCACTGAAGCCACCATACCCCCGGCtttt                                                             <  1:987476/100‑1 (MQ=255)
                                    atGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTtata                                                            >  2:1985662/1‑100 (MQ=255)
                                     tGAAACCGGGAATATTTTGCGGGGAGAACACCTGAATCAGTAACCAGGTAATGGGAGAACAACCCTGATCCATTGAAGCCACCATACCCCCGGc                                                                 >  2:1396292/1‑94 (MQ=255)
                                     tGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTaa                                                           <  2:1043976/100‑1 (MQ=255)
                                     tGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTa                                                            <  2:1979756/99‑1 (MQ=255)
                                       aaaCCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACg                                                         <  1:2080174/100‑1 (MQ=255)
                                         aCCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGcc                                                       <  1:2441294/100‑1 (MQ=255)
                                           cGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTGACGcccc                                                     <  1:1346101/100‑1 (MQ=255)
                                             ggAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAg                                                   <  2:2095052/100‑1 (MQ=255)
                                             ggAATATTATGCGGGGGGAAAACCGGAAACAGTAACCAGGTAATGGGCGAACAGCCCTGATCAATTGAGGCAACCATAACGCCGGCTTTTAACGCCCCa                                                    >  2:1727608/1‑99 (MQ=255)
                                               aatcgtttGCGGGGAGAAAGCCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCTGGCGTTTAACGCCCCAGCa                                                 <  2:2175814/95‑1 (MQ=255)
                                                atatTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCAt                                                >  2:980163/1‑100 (MQ=255)
                                                 tatTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCACAGCAtt                                               >  1:1109477/1‑100 (MQ=255)
                                                     ttGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAgcc                                           <  2:98977/100‑1 (MQ=255)
                                                       gCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAgccgc                                         <  2:550977/100‑1 (MQ=255)
                                                       gCGGGGAGAAAACATCCATCAGTAACCGGGTAATAGGACAAAAGCCCTCATCCATGGAAGCCACCATACCCCCGGCTTTTAACGCACCAGCATTAgccgc                                         <  2:508570/100‑1 (MQ=255)
                                                         ggggCGAACGCCTGAATCAGTCACCCGGTCATGGGAGCAATGGCCTGATCCATTGAAGCTACCATACCTTCGGCTTTTAACGCCCCatttgtag                                             <  2:1491549/94‑8 (MQ=255)
                                                          gggAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAg                                      <  2:1720757/100‑1 (MQ=255)
                                                            gcgaAAGCCTGAAATAGTAACCAGGAAATGGGAGACCAGACCTGACCCATTGACGCCACCCTGCCCCAGGCATATAACGCCCCAGCATCAGCCGCCAGtt                                    <  2:2397509/98‑1 (MQ=255)
                                                            gagaAAACCTGATTCAGTAACCAGGTAATGGTAGAACAGCCCTTATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGt                                     >  2:2138600/1‑99 (MQ=255)
                                                             agaAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTg                                   <  2:24112/100‑1 (MQ=255)
                                                              gaAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTgg                                  >  2:1191224/1‑100 (MQ=255)
                                                                aaaCCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCCTTTAACGCCCCAGCATTAGCCGCCAGTCGGgt                                >  2:28252/1‑100 (MQ=255)
                                                                 aaCCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCCGCATTAGCCGCCAGTTGGgtg                               >  2:144142/1‑100 (MQ=255)
                                                                 aaCCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGgtg                               <  2:683628/100‑1 (MQ=255)
                                                                 aaCCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGgtg                               <  2:216195/100‑1 (MQ=255)
                                                                 aaCCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGgtg                               <  1:1095685/100‑1 (MQ=255)
                                                                   ccTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGtgtg                             <  2:1047585/100‑1 (MQ=255)
                                                                   ccTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGtgtg                             <  2:1142886/100‑1 (MQ=255)
                                                                    cTGATTCGGTAACCAGGTACTGGGCGAACAGCCTTGATCCATTGAAACCACCATACCCCCGGCTTTTAACTCCCCAGCATTAGCCGCCAGTTGGgtgta                             <  2:2523264/99‑2 (MQ=255)
                                                                    cTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGc                            <  2:652296/100‑1 (MQ=255)
                                                                     tGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCa                           <  2:520824/100‑1 (MQ=255)
                                                                         tCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAAcc                       <  1:1375582/100‑1 (MQ=255)
                                                                           aGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCaa                     <  2:2538022/100‑1 (MQ=255)
                                                                           aGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCaa                     <  2:621157/100‑1 (MQ=255)
                                                                               aCCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGt                 >  1:2012411/1‑100 (MQ=255)
                                                                                ccAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGtt                >  1:246249/1‑100 (MQ=255)
                                                                                ccAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGtt                >  1:899574/1‑100 (MQ=255)
                                                                                  aGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTg              <  1:810993/100‑1 (MQ=255)
                                                                                   ggTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGc             <  1:1913060/100‑1 (MQ=255)
                                                                                   ggTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGc             <  1:1080292/100‑1 (MQ=255)
                                                                                    gTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCg            >  1:2244435/1‑100 (MQ=255)
                                                                                    gTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCg            <  2:1215079/100‑1 (MQ=255)
                                                                                    gTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCg            >  1:1032047/1‑100 (MQ=255)
                                                                                     tAATGGGAGAACAGCCCTGATCCATTGAAGCCCCCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTCTGc             <  2:2041858/98‑1 (MQ=255)
                                                                                        tGGGAGACCAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACCCCCCAGCAGTAGCCGCCAGTTGGGTGTGCAAACCAATAGTGTGCGTcgc        <  2:814614/100‑1 (MQ=255)
                                                                                         gggAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTCTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTcgcg       >  2:197299/1‑100 (MQ=255)
                                                                                             gaACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATcc   <  1:1512842/100‑1 (MQ=255)
                                                                                              aaCAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCca  >  2:1534242/1‑100 (MQ=255)
                                                                                                 |                                                                                                
ATGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCA  >  NC_000913/2172979‑2173172

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