Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,786,729 G→A pseudogene (333/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,786,7290GA100.0% 252.9 / NA 77pseudogene (333/2253 nt)ygaQuncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/40);  total (37/40)

AGAAAAAGTTATTAATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTA  >  NC_000913/2786594‑2786862
                                                                                                                                       |                                                                                                                                     
aGAAAAAGTTATTAATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGtt                                                                                                                                    >  1:529045/1‑139 (MQ=255)
                gATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATaaat                                                                                                                    >  2:932577/1‑136 (MQ=255)
                 aTCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATaaata                                                                                                                   >  2:239115/1‑135 (MQ=255)
                        aaaaaaCGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGaaa                                                                                                           >  1:1018081/1‑139 (MQ=255)
                          aaaaCGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAaga                                                                                                         <  1:932577/139‑1 (MQ=255)
                              cGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTg                                                                                                     >  1:592631/1‑139 (MQ=255)
                                aTTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATc                                                                                                                                >  2:547504/1‑111 (MQ=255)
                                aTTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATc                                                                                                                                <  1:547504/111‑1 (MQ=255)
                                aTTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGAt                                                                                                   >  2:307064/1‑139 (MQ=255)
                                     aaaaTTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                              <  1:307064/139‑1 (MQ=255)
                                          ttCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGgtg                                                                                         <  1:1177/139‑1 (MQ=255)
                                          ttCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGgtg                                                                                         <  1:1215595/139‑1 (MQ=255)
                                          ttCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGgtg                                                                                         <  1:106889/139‑1 (MQ=255)
                                               tttCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTAt                                                                                    <  2:592631/139‑1 (MQ=255)
                                               tttCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTAt                                                                                    <  1:35817/139‑1 (MQ=255)
                                                 tctcCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATc                                                                                                                                <  2:1162612/94‑1 (MQ=255)
                                                 tctcCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATc                                                                                                                                >  1:1162612/1‑94 (MQ=255)
                                                           cAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGt                                                                        >  2:172725/1‑139 (MQ=255)
                                                             aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCgg                                                                                           >  2:1164845/1‑118 (MQ=255)
                                                             aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCgg                                                                                           <  1:1164845/118‑1 (MQ=255)
                                                              cTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCt                                                                     <  2:691091/139‑1 (MQ=255)
                                                                 aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGgtg                                                                                         >  2:1297150/1‑116 (MQ=255)
                                                                 aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGgtg                                                                                         <  1:1297150/116‑1 (MQ=255)
                                                                  ccGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGAGGTGTTTATCAGTATGGTCGTGCTGAAc                                                                 >  2:1221853/1‑139 (MQ=255)
                                                                   cGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCg                                                                                            <  1:323134/111‑1 (MQ=255)
                                                                   cGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCg                                                                                            >  2:323134/1‑111 (MQ=255)
                                                                       agTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTc                                                            >  2:693216/1‑139 (MQ=255)
                                                                        gTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCa                                                           >  1:1258699/1‑139 (MQ=255)
                                                                         ttGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                                                          >  2:1131253/1‑139 (MQ=255)
                                                                          tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                         <  1:693216/139‑1 (MQ=255)
                                                                          tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                         <  1:1017354/139‑1 (MQ=255)
                                                                           gAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAccc                                                        <  2:873776/139‑1 (MQ=255)
                                                                            aaTAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTAt                                                                                    >  2:666662/1‑110 (MQ=255)
                                                                            aaTAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTAt                                                                                    <  1:666662/110‑1 (MQ=255)
                                                                                 agcagcACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGGAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATta                                                  >  1:680634/1‑139 (MQ=255)
                                                                                 agcagcACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATta                                                  >  2:1073214/1‑139 (MQ=255)
                                                                                   cagcaCAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATtata                                                <  2:1097021/139‑1 (MQ=255)
                                                                                      cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCcaca                                                     <  2:406279/131‑1 (MQ=255)
                                                                                      cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCcaca                                                     >  1:406279/1‑131 (MQ=255)
                                                                                           cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                      <  2:624393/93‑1 (MQ=255)
                                                                                           cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                      >  1:624393/1‑93 (MQ=255)
                                                                                           cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACaa                                        >  2:1144919/1‑139 (MQ=255)
                                                                                            gTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACaaa                                       <  2:991649/139‑1 (MQ=255)
                                                                                                     cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCcaca                                                     <  2:767072/116‑1 (MQ=255)
                                                                                                     cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCcaca                                                     >  1:767072/1‑116 (MQ=255)
                                                                                                      aGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaa                             <  2:529045/139‑1 (MQ=255)
                                                                                                          ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaa                             <  2:430337/135‑1 (MQ=255)
                                                                                                          ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaa                             >  1:430337/1‑135 (MQ=255)
                                                                                                             gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATaa                      >  1:148355/1‑139 (MQ=255)
                                                                                                             gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATaa                      <  1:1026973/139‑1 (MQ=255)
                                                                                                             gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATaa                      <  2:1258699/139‑1 (MQ=255)
                                                                                                             gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATaa                      >  2:1026973/1‑139 (MQ=255)
                                                                                                                  gATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATg                 <  2:935618/139‑1 (MQ=255)
                                                                                                                   aTAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACAttt                                                             <  1:373468/94‑1 (MQ=255)
                                                                                                                   aTAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACAttt                                                             >  2:373468/1‑94 (MQ=255)
                                                                                                                      ggAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGc             <  1:1073214/139‑1 (MQ=255)
                                                                                                                         gaaagaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCata          <  2:148355/139‑1 (MQ=255)
                                                                                                                         gaaagaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCata          <  1:1144919/139‑1 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGtt                                                                                       <  2:1087825/59‑1 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGtt                                                                                       >  1:1087825/1‑59 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                         <  2:1217873/89‑1 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                         >  2:925292/1‑89 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                         >  1:1217873/1‑89 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                         <  1:925292/89‑1 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATta                                                  >  1:410162/1‑96 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATta                                                  <  2:410162/96‑1 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAatat                  >  2:1153978/1‑128 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAatat                  <  1:1153978/128‑1 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGc             <  1:1185500/133‑1 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGc             >  2:542597/1‑133 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGc             <  1:542597/133‑1 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGc             >  2:1185500/1‑133 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc       >  1:686203/1‑139 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc       >  2:226154/1‑139 (MQ=255)
                                                                                                                                 aaCATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGgtg                                                                                         <  2:824036/52‑1 (MQ=255)
                                                                                                                                 aaCATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGgtg                                                                                         >  1:824036/1‑52 (MQ=255)
                                                                                                                                 aaCATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTa  <  2:92814/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
AGAAAAAGTTATTAATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTA  >  NC_000913/2786594‑2786862

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