Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,794,015 G→A *427* (TAG→TAA gabT → 4‑aminobutyrate aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,794,0150GA100.0% 238.1 / NA 73*427* (TAG→TAAgabT4‑aminobutyrate aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/34);  total (39/34)

AGGCCTGATTCTTCTCTCCTGCGGCCCGTATTACAACGTGCTGCGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGG  >  NC_000913/2793883‑2794146
                                                                                                                                    |                                                                                                                                   
aGGCCTGATTCTTCTCTCCTGCGGCCCGTATTACAACGTGCTGCGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAAcgccgc                                                                                                                               <  2:80375/139‑1 (MQ=255)
                                    cGTGCTGCGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTg                                                                                                         >  2:1164271/1‑125 (MQ=255)
                                    cGTGCTGCGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTg                                                                                                         <  1:1164271/125‑1 (MQ=255)
                                      tgctgcGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTcc                                                                                                                            >  2:913260/1‑104 (MQ=255)
                                      tgctgcGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTcc                                                                                                                            <  1:913260/104‑1 (MQ=255)
                                      tgctgcGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTc                                                                                                                             <  2:1357263/103‑1 (MQ=255)
                                      tgctgcGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTc                                                                                                                             >  1:1357263/1‑103 (MQ=255)
                                      tgctgcGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGCCAGGGTCTGGAGATCATCAGCCAGTGTATTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGcc                                                                                         >  2:1116239/1‑139 (MQ=255)
                                       gctgcGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGCTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAAcgccgc                                                                                                                               <  2:1027433/100‑1 (MQ=255)
                                       gctgcGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAAcgccgc                                                                                                                               >  1:1027433/1‑100 (MQ=255)
                                           cgcATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGc                                                                                                           >  2:191076/1‑116 (MQ=255)
                                           cgcATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGc                                                                                                           <  1:191076/116‑1 (MQ=255)
                                                          cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGt                                                                                                    <  2:1528968/108‑1 (MQ=255)
                                                          cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGt                                                                                                    >  1:1528968/1‑108 (MQ=255)
                                                          cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCg                                                                                            >  1:65012/1‑116 (MQ=255)
                                                          cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCg                                                                                            <  2:65012/116‑1 (MQ=255)
                                                          cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGcc                                                                                         >  1:1506176/1‑119 (MQ=255)
                                                          cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGcc                                                                                         <  2:1506176/119‑1 (MQ=255)
                                                          cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGAt                                                                     >  2:68773/1‑139 (MQ=255)
                                                          cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGAt                                                                     >  2:1438163/1‑139 (MQ=255)
                                                          cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGAt                                                                     >  1:1595221/1‑139 (MQ=255)
                                                          cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGAt                                                                     >  1:1249892/1‑139 (MQ=255)
                                                          cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGAt                                                                     >  1:151199/1‑139 (MQ=255)
                                                          cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGAt                                                                     >  2:72799/1‑139 (MQ=255)
                                                          cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGAt                                                                     >  1:788355/1‑139 (MQ=255)
                                                          cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTACGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGAt                                                                     >  1:453357/1‑139 (MQ=255)
                                                          cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGACGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGAt                                                                     >  2:1287023/1‑139 (MQ=255)
                                                           tCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATa                                                                    >  2:214367/1‑139 (MQ=255)
                                                                  tGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCTGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGtt                                                             >  2:991846/1‑139 (MQ=255)
                                                                         gCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCt                                                                                                                           <  2:227288/70‑1 (MQ=255)
                                                                         gCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCt                                                                                                                           >  1:227288/1‑70 (MQ=255)
                                                                         gCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgccg                                                      >  1:1025274/1‑139 (MQ=255)
                                                                         gCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgccg                                                      <  1:214367/139‑1 (MQ=255)
                                                                          cTCAGATCCGCCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGACGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATTAGGCGTTTAcgccgc                                                     <  1:1116239/139‑1 (MQ=255)
                                                                                tCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                               <  1:72799/139‑1 (MQ=255)
                                                                                  cGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACa                                             <  2:788355/139‑1 (MQ=255)
                                                                                    tCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                               >  2:570883/1‑135 (MQ=255)
                                                                                    tCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                               <  1:570883/135‑1 (MQ=255)
                                                                                         gTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAAcgcc                                                                                                                                 <  1:903696/48‑1 (MQ=255)
                                                                                         gTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAAcgcc                                                                                                                                 >  2:903696/1‑48 (MQ=255)
                                                                                         gTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTa                                      >  1:1228575/1‑139 (MQ=255)
                                                                                              gagaTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgccgc                                                     >  2:996288/1‑119 (MQ=255)
                                                                                              gagaTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgccgc                                                     <  1:996288/119‑1 (MQ=255)
                                                                                                gaTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCTTTTACGCCGCATCCGACAATCTGTACGTGAAc                               <  1:991846/139‑1 (MQ=255)
                                                                                                gaTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAAc                               <  2:1249892/139‑1 (MQ=255)
                                                                                                   catcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGCCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAgg                            >  2:630489/1‑139 (MQ=255)
                                                                                                    atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGga                           >  1:893101/1‑139 (MQ=255)
                                                                                                    atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTACGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGga                           <  2:453357/139‑1 (MQ=255)
                                                                                                     tcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGCCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaa                          <  1:630489/139‑1 (MQ=255)
                                                                                                      caGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaag                         >  1:1660860/1‑139 (MQ=255)
                                                                                                        gCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCt                                                                                                                              >  1:952232/1‑36 (MQ=255)
                                                                                                        gCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCt                                                                                                                              <  2:952232/36‑1 (MQ=255)
                                                                                                          cAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAAt                     <  1:68773/139‑1 (MQ=255)
                                                                                                             tgtTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATc                                                                               >  1:1154054/1‑78 (MQ=255)
                                                                                                             tgtTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATc                                                                               <  2:1154054/78‑1 (MQ=255)
                                                                                                             tgtTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaa                          >  2:507097/1‑131 (MQ=255)
                                                                                                             tgtTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaa                          <  1:507097/131‑1 (MQ=255)
                                                                                                             tgtTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTa                  >  1:202035/1‑139 (MQ=255)
                                                                                                             tgtTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTa                  >  2:923548/1‑139 (MQ=255)
                                                                                                               ttttGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATg                <  2:1660860/139‑1 (MQ=255)
                                                                                                                tttGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGt               <  2:151199/139‑1 (MQ=255)
                                                                                                                  tgatgaGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTc                                                                                              >  2:1210743/1‑58 (MQ=255)
                                                                                                                  tgatgaGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTc                                                                                              <  1:1210743/58‑1 (MQ=255)
                                                                                                                   gatgaGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTgg            <  2:1025274/139‑1 (MQ=255)
                                                                                                                   gatgaGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTgg            <  2:893101/139‑1 (MQ=255)
                                                                                                                    atgaGGCGAAGCAGTAACGCCGCTCCTAAGCCGGGGGCGACGCTGCGCGACTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGc           >  1:1080172/1‑139 (MQ=255)
                                                                                                                       aGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTg                                                                                                         >  2:1543626/1‑42 (MQ=255)
                                                                                                                       aGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTg                                                                                                         <  1:1543626/42‑1 (MQ=255)
                                                                                                                       aGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTAc                                                                                      >  2:1524144/1‑61 (MQ=255)
                                                                                                                       aGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTAc                                                                                      <  1:1524144/61‑1 (MQ=255)
                                                                                                                       aGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGCATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaa                          <  2:1058062/121‑1 (MQ=255)
                                                                                                                       aGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGCATCCAGGTCGGAAAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaa                          >  1:1058062/1‑121 (MQ=255)
                                                                                                                             aGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAgg  <  2:202035/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                   
AGGCCTGATTCTTCTCTCCTGCGGCCCGTATTACAACGTGCTGCGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGG  >  NC_000913/2793883‑2794146

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: