Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,788,238 G→A pseudogene (1842/2253 nt) ygaQ → uncharacterized protein, putative enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,788,2380GA100.0% 258.0 / NA 79pseudogene (1842/2253 nt)ygaQuncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/42);  total (37/42)

ACAATTAGCAATGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAGCTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCTGGCCTTTAG  >  NC_000913/2788106‑2788363
                                                                                                                                    |                                                                                                                             
aCAATTAGCAATGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTATAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTtaat                                                                                                                         <  1:1099067/139‑1 (MQ=255)
      aGCAATGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGt                                                                                                                   >  1:803095/1‑139 (MQ=255)
        cAATGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTg                                                                                                                 <  2:1321750/139‑1 (MQ=255)
         aaTGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGc                                                                                                                >  2:1536839/1‑139 (MQ=255)
            ggCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGctct                                                                                                             >  2:523178/1‑139 (MQ=255)
              cATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTAc                                                                                                           >  2:947675/1‑139 (MQ=255)
                tatataGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCt                                                                                                         >  2:791318/1‑139 (MQ=255)
                     aGCCTCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                    <  1:1526416/139‑1 (MQ=255)
                     aGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGctc                                                                                                              >  2:1307245/1‑129 (MQ=255)
                     aGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGctc                                                                                                              <  1:1307245/129‑1 (MQ=255)
                     aGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                    <  1:523178/139‑1 (MQ=255)
                       ccGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTg                                                                                                                 >  2:1553390/1‑124 (MQ=255)
                       ccGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTg                                                                                                                 <  1:1553390/124‑1 (MQ=255)
                       ccGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGctc                                                                                                              >  1:706622/1‑127 (MQ=255)
                       ccGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGctc                                                                                                              <  2:706622/127‑1 (MQ=255)
                       ccGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGct                                                                                                               <  2:1316042/126‑1 (MQ=255)
                       ccGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGct                                                                                                               >  1:1316042/1‑126 (MQ=255)
                       ccGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcagca                                                                                                  >  1:907128/1‑139 (MQ=255)
                        cGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcagcag                                                                                                 <  2:1410050/139‑1 (MQ=255)
                          ctctGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGct                                                                                               >  2:671862/1‑139 (MQ=255)
                          ctctGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGct                                                                                               >  2:1292147/1‑139 (MQ=255)
                          ctctGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGct                                                                                               >  2:852600/1‑139 (MQ=255)
                            ctGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGctct                                                                                             >  2:936410/1‑139 (MQ=255)
                                cGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACaa                                                                                         <  1:947675/139‑1 (MQ=255)
                                        agGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTcc                                                                                 <  2:803095/139‑1 (MQ=255)
                                        agGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTcc                                                                                 <  1:1176566/139‑1 (MQ=255)
                                         ggCCTTTAGTTTATTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCa                                                                                >  1:1077450/1‑139 (MQ=255)
                                         ggCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCa                                                                                <  1:1292147/139‑1 (MQ=255)
                                          gCCTTTAGTTTATTGGGATGGCCAAGTATTAAAATCAACGTCTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAGTAAAGATTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                               <  2:1077450/139‑1 (MQ=255)
                                          gCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAGAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                               <  1:1287089/139‑1 (MQ=255)
                                          gCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCGGCTCTACAATGAATTCCAc                                                                               <  2:1073752/139‑1 (MQ=255)
                                          gCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                               <  1:791318/139‑1 (MQ=255)
                                          gCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                               <  1:1536839/139‑1 (MQ=255)
                                           ccTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACg                                                                              <  2:907128/139‑1 (MQ=255)
                                           ccTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACg                                                                              >  2:569981/1‑139 (MQ=255)
                                            cTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCTACGc                                                                             <  1:936410/139‑1 (MQ=255)
                                            cTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGc                                                                             <  1:569981/139‑1 (MQ=255)
                                              ttAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCAt                                                                           <  1:868190/139‑1 (MQ=255)
                                              ttAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCAt                                                                           <  2:1585718/139‑1 (MQ=255)
                                              ttAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCAt                                                                           <  2:48722/139‑1 (MQ=255)
                                                   ttACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATtatat                                                                      <  2:677861/139‑1 (MQ=255)
                                                    tACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                            <  1:1143202/132‑1 (MQ=255)
                                                    tACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                            >  2:1143202/1‑132 (MQ=255)
                                                     aCTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATtatatat                                                                    >  1:122295/1‑139 (MQ=255)
                                                        gggATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGa                                                                 <  2:122295/139‑1 (MQ=255)
                                                         ggATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGctc                                                                                                              <  1:645126/93‑1 (MQ=255)
                                                         ggATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGctc                                                                                                              >  2:645126/1‑93 (MQ=255)
                                                           aTGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                    <  2:1162060/101‑1 (MQ=255)
                                                           aTGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                    >  1:1162060/1‑101 (MQ=255)
                                                           aTGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                               >  1:1334933/1‑122 (MQ=255)
                                                           aTGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                               <  2:1334933/122‑1 (MQ=255)
                                                             ggCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAg                                                            <  1:852600/139‑1 (MQ=255)
                                                              gCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcagcag                                                                                                 <  2:1103100/101‑1 (MQ=255)
                                                              gCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcagcag                                                                                                 >  1:1103100/1‑101 (MQ=255)
                                                                        aaaaTCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATgg                                                 >  2:1608539/1‑139 (MQ=255)
                                                                              aaCGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATaaa                                                                                                                      <  1:756657/64‑1 (MQ=255)
                                                                              aaCGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATaaa                                                                                                                      >  2:756657/1‑64 (MQ=255)
                                                                              aaCGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                               <  1:1559250/103‑1 (MQ=255)
                                                                              aaCGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                               >  2:1559250/1‑103 (MQ=255)
                                                                              aaCGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                            >  1:878364/1‑106 (MQ=255)
                                                                              aaCGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                            <  2:878364/106‑1 (MQ=255)
                                                                                                         gTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAgc                                              <  1:684289/109‑1 (MQ=255)
                                                                                                         gTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAgc                                              >  2:684289/1‑109 (MQ=255)
                                                                                                          tCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                    >  2:881690/1‑54 (MQ=255)
                                                                                                          tCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                    <  1:881690/54‑1 (MQ=255)
                                                                                                          tCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATc                                                                                                      >  2:560082/1‑52 (MQ=255)
                                                                                                          tCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATc                                                                                                      <  1:560082/52‑1 (MQ=255)
                                                                                                              cgcgTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCt           >  2:602751/1‑139 (MQ=255)
                                                                                                               gcgTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCTg          <  1:602751/139‑1 (MQ=255)
                                                                                                                   gtCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCTGGCCt      >  1:1591038/1‑139 (MQ=255)
                                                                                                                    tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCa                                                                                >  1:1174840/1‑64 (MQ=255)
                                                                                                                    tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCa                                                                                <  2:1174840/64‑1 (MQ=255)
                                                                                                                    tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCAt                                                                           <  2:848177/69‑1 (MQ=255)
                                                                                                                    tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCAt                                                                           >  1:848177/1‑69 (MQ=255)
                                                                                                                    tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCTGGCCtt     >  2:1198251/1‑139 (MQ=255)
                                                                                                                    tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCTGGCCtt     >  2:903179/1‑139 (MQ=255)
                                                                                                                    tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCTGGCCtt     >  2:47235/1‑139 (MQ=255)
                                                                                                                      ggcggTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCTGGCCTTTa   <  1:1608539/139‑1 (MQ=255)
                                                                                                                       gcggTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTAGTAGCAACTAAAAACGTCCTGGACTTTAg  >  1:1333871/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                             
ACAATTAGCAATGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAGCTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCTGGCCTTTAG  >  NC_000913/2788106‑2788363

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