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% 258.9 / NA 80pseudogene (333/2253 nt)ygaQuncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/40);  total (40/40)

CAGAAAAAGTTATTAATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAACA  >  NC_000913/2786593‑2786865
                                                                                                                                        |                                                                                                                                        
cAGAAAAAGTTATTAATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTACAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGt                                                                                                                                        <  1:1015068/139‑1 (MQ=255)
       aGTTATTAATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCaa                                                                                                                                 <  2:323854/139‑1 (MQ=255)
                          aaaaaCGAATTTAAAAATTTCCCTTCTCCAAATACAACTTACCGAAAGTTTAAAAAAGGAGCACAACGTAGCAATAAAGCGCCCGACAGGAAAAGAGAAAGAAAAACATAAGTTTATCAACGAGATAAAAGAATGAAag                                                                                                             <  2:962088/139‑1 (MQ=255)
                                      aaaaTTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                                 <  2:398417/139‑1 (MQ=255)
                                            tCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATc                                                                                                                                   >  2:757139/1‑100 (MQ=255)
                                            tCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATc                                                                                                                                   <  1:757139/100‑1 (MQ=255)
                                            tCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                                 <  2:577372/133‑1 (MQ=255)
                                            tCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                                 >  1:577372/1‑133 (MQ=255)
                                                            cAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGt                                                                           >  2:723846/1‑139 (MQ=255)
                                                             aaCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTg                                                                          <  1:581046/139‑1 (MQ=255)
                                                              aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATaaa                                                                                                                        >  2:847113/1‑91 (MQ=255)
                                                              aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATaaa                                                                                                                        <  1:847113/93‑3 (MQ=255)
                                                                  aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGaa                                                                     >  1:298001/1‑139 (MQ=255)
                                                                    cGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACa                                                                   >  1:50918/1‑139 (MQ=255)
                                                                       gagTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACAttt                                                                <  2:50918/139‑1 (MQ=255)
                                                                         gTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCg                                                                                               >  2:971900/1‑106 (MQ=255)
                                                                         gTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCg                                                                                               <  1:971900/106‑1 (MQ=255)
                                                                         gTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGgtg                                                                                            <  2:1044649/109‑1 (MQ=255)
                                                                         gTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGgtg                                                                                            >  1:1044649/1‑109 (MQ=255)
                                                                          ttGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                                                             <  2:723082/139‑1 (MQ=255)
                                                                           tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGATATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                            <  2:803322/139‑1 (MQ=255)
                                                                                  agcagcACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGgtg                                                                                            <  2:601927/100‑1 (MQ=255)
                                                                                  agcagcACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGgtg                                                                                            >  1:601927/1‑100 (MQ=255)
                                                                                    cagcaCAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATtata                                                   >  2:316042/1‑139 (MQ=255)
                                                                                       cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGa                                                                                                                <  2:769592/75‑1 (MQ=255)
                                                                                       cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGa                                                                                                                >  1:769592/1‑75 (MQ=255)
                                                                                       cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                            <  2:676531/127‑1 (MQ=255)
                                                                                       cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                            >  1:676531/1‑127 (MQ=255)
                                                                                            cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                         <  1:895502/93‑1 (MQ=255)
                                                                                            cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                         >  2:895502/1‑93 (MQ=255)
                                                                                            cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACaa                                           >  1:748917/1‑139 (MQ=255)
                                                                                            cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACaa                                           >  1:483237/1‑139 (MQ=255)
                                                                                               aGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAAcc                                        >  1:670689/1‑139 (MQ=255)
                                                                                                      cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGtt                                                                                          >  1:993240/1‑82 (MQ=255)
                                                                                                      cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGtt                                                                                          <  2:993240/82‑1 (MQ=255)
                                                                                                      cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCa                                 <  2:298001/139‑1 (MQ=255)
                                                                                                      cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCa                                 >  1:498829/1‑139 (MQ=255)
                                                                                                         cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                            >  1:1038660/1‑109 (MQ=255)
                                                                                                         cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                            <  2:1038660/109‑1 (MQ=255)
                                                                                                           ccAGACAGGCTAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGTCCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCTCCAAAAGa                            >  2:1063901/1‑139 (MQ=255)
                                                                                                           ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                                                             >  1:359378/1‑106 (MQ=255)
                                                                                                           ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                                                             <  2:359378/106‑1 (MQ=255)
                                                                                                              gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCg                                                                                               <  1:405202/69‑1 (MQ=255)
                                                                                                              gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCg                                                                                               >  2:405202/1‑69 (MQ=255)
                                                                                                              gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATaa                         >  1:276479/1‑139 (MQ=255)
                                                                                                              gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATaa                         >  2:256352/1‑139 (MQ=255)
                                                                                                              gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATAGGAAGAGGTGATAACAGCGGTCTTTATCAGTATGGTCGTGCTGAACATTTAACCCACATTATATCAGACACACCTTCCCCAATAGATaa                         >  2:261152/1‑139 (MQ=255)
                                                                                                                   gATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATg                    <  2:748917/139‑1 (MQ=255)
                                                                                                                   gATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATg                    <  2:604421/139‑1 (MQ=255)
                                                                                                                      aGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGgtg                                                                                            >  2:307919/1‑64 (MQ=255)
                                                                                                                      aGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGgtg                                                                                            <  1:307919/64‑1 (MQ=255)
                                                                                                                         agaaagaaAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCat              <  2:276479/139‑1 (MQ=255)
                                                                                                                          gaaagaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCata             <  1:316042/139‑1 (MQ=255)
                                                                                                                           aaagaaagACATAAGTTTATCAACTAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                            >  1:656300/1‑91 (MQ=255)
                                                                                                                           aaagaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                            <  2:656300/91‑1 (MQ=255)
                                                                                                                           aaagaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCatat            >  2:376039/1‑139 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                                 >  2:311594/1‑52 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                                 <  2:140793/52‑1 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                                 <  1:311594/52‑1 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                                 >  1:140793/1‑52 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                            >  1:50452/1‑89 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                            <  1:134678/89‑1 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                            <  2:50452/89‑1 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                            >  2:134678/1‑89 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATg                    <  2:267512/129‑1 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATg                    <  2:1005878/129‑1 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATg                    >  1:1005878/1‑129 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATg                    >  1:267512/1‑129 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCata             >  1:146334/1‑136 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCata             <  2:146334/136‑1 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc          >  2:481048/1‑139 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc          >  1:837942/1‑139 (MQ=255)
                                                                                                                             agaaaAACATAAGTTTAGCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTAGGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc          >  1:671260/1‑139 (MQ=255)
                                                                                                                                  aaCATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                                 >  2:1025467/1‑47 (MQ=38)
                                                                                                                                  aaCATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                                 <  1:1025467/47‑1 (MQ=38)
                                                                                                                                  aaCATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCcaca                                                        >  1:18286/1‑88 (MQ=255)
                                                                                                                                  aaCATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCcaca                                                        <  2:18286/88‑1 (MQ=255)
                                                                                                                                  aaCATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTa     <  2:837942/139‑1 (MQ=255)
                                                                                                                                  aaCATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTa     <  1:723846/139‑1 (MQ=255)
                                                                                                                                     aTAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAAcg  <  2:671260/139‑2 (MQ=255)
                                                                                                                                        |                                                                                                                                        
CAGAAAAAGTTATTAATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAACA  >  NC_000913/2786593‑2786865

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