Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,767,179 G→A *281* (TAG→TAA yibA → putative immunity protein for polymorphic toxin RhsA, HEAT‑domain protein, lethality reduction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,767,1790GA100.0% 246.8 / NA 76*281* (TAG→TAAyibAputative immunity protein for polymorphic toxin RhsA, HEAT‑domain protein, lethality reduction protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (44/32);  total (44/32)

GTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCG  >  NC_000913/3767045‑3767311
                                                                                                                                      |                                                                                                                                    
gTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTg                                                                                                                                  >  1:1079988/1‑139 (MQ=255)
      gatgatATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                            <  2:1079988/139‑1 (MQ=255)
          atatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTAt                                                                                                                        <  2:416427/139‑1 (MQ=255)
           tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATg                                                                                                                       >  1:565255/1‑139 (MQ=255)
           tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGAAAATGAAATTATAACTTCCGCTATTGATAACCTGAAGCGTTCATAACGTGAGTTGCCTATg                                                                                                                       >  2:764946/1‑139 (MQ=255)
                ttGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                               <  2:1159967/126‑1 (MQ=255)
                ttGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                               >  1:1159967/1‑126 (MQ=255)
                     gCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                               >  1:363551/1‑121 (MQ=255)
                     gCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                               <  2:363551/121‑1 (MQ=255)
                     gCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTACCTATGCACAGTgggg                                                                                                             >  2:25690/1‑139 (MQ=255)
                         cGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTCGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGa                                                                                                                                 <  1:177293/115‑1 (MQ=255)
                         cGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTCGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGa                                                                                                                                 >  2:177293/1‑115 (MQ=255)
                          gggTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                               >  2:732582/1‑116 (MQ=255)
                          gggTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                               <  1:732582/116‑1 (MQ=255)
                          gggTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTc                                                                                                        >  2:176735/1‑139 (MQ=255)
                                   aGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                               <  2:565255/139‑1 (MQ=255)
                                         tAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGc                                                                                         <  1:764946/139‑1 (MQ=255)
                                         tAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAAATTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGc                                                                                         >  1:799994/1‑139 (MQ=255)
                                              cGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                               <  1:384985/96‑1 (MQ=255)
                                              cGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                               >  2:384985/1‑96 (MQ=255)
                                                   cTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCgg                                                                                                 <  2:845853/121‑1 (MQ=255)
                                                   cTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCgg                                                                                                 >  1:845853/1‑121 (MQ=255)
                                                   cTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGa                                                                               >  2:414790/1‑139 (MQ=255)
                                                               gATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                     <  1:1027940/89‑1 (MQ=255)
                                                               gATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                     >  2:1027940/1‑89 (MQ=255)
                                                                         tgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGa                                                                                                                                 <  1:609862/67‑1 (MQ=255)
                                                                         tgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGa                                                                                                                                 >  2:609862/1‑67 (MQ=255)
                                                                         tgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTg                                                         >  1:893293/1‑139 (MQ=255)
                                                                         tgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTg                                                         >  2:353200/1‑139 (MQ=255)
                                                                          gtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTgg                                                        <  2:893293/139‑1 (MQ=255)
                                                                           tACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                               >  2:474336/1‑67 (MQ=255)
                                                                           tACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                               <  1:474336/67‑1 (MQ=255)
                                                                           tACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGc                                                       >  2:354328/1‑139 (MQ=255)
                                                                           tACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGc                                                       >  1:317688/1‑139 (MQ=255)
                                                                            aCAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGCTGCTCGACCGGCTGGAAAGTGAAATCCTGGCt                                                      <  2:799994/139‑1 (MQ=255)
                                                                                gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGa                                                                                                                                 <  2:269348/60‑1 (MQ=255)
                                                                                gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGa                                                                                                                                 >  1:269348/1‑60 (MQ=255)
                                                                                gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTggggg                                                                                                            >  1:166856/1‑81 (MQ=255)
                                                                                gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTggggg                                                                                                            <  2:166856/81‑1 (MQ=255)
                                                                                gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGAcc                                                  >  1:109788/1‑139 (MQ=255)
                                                                                   tgatgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCggg                                               <  2:634540/139‑1 (MQ=255)
                                                                                    gatgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgt                                              <  2:109788/139‑1 (MQ=255)
                                                                                     atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                               >  2:1097777/1‑57 (MQ=255)
                                                                                     atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                               <  1:1097777/57‑1 (MQ=255)
                                                                                     atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCgg                                                >  1:786746/1‑136 (MQ=255)
                                                                                     atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCgg                                                <  2:786746/136‑1 (MQ=255)
                                                                                     atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtg                                             >  1:1094935/1‑139 (MQ=255)
                                                                                          aaTGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                     >  2:1220087/1‑62 (MQ=255)
                                                                                          aaTGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                     <  1:1220087/62‑1 (MQ=255)
                                                                                               aaTTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGaaa                                   >  1:1323/1‑139 (MQ=255)
                                                                                                  tataACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAgcc                                >  2:1135983/1‑139 (MQ=255)
                                                                                                    taACTTCCGCTATTGATAATCTGAAGCGTGCATAACGTGAGTTGCCTATGCACAGTGGGGGCTTCACGCCGGCACTGAGCAGATGCTCGACGGGGTGGAAAGTGAAATTCTGGCTGACCGGGGGAGTGAGGAAAgccgc                              >  2:360116/1‑139 (MQ=25)
                                                                                                      aCTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                            >  2:1224691/1‑43 (MQ=255)
                                                                                                      aCTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                            <  1:1224691/43‑1 (MQ=255)
                                                                                                      aCTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                            <  2:348744/139‑1 (MQ=255)
                                                                                                          ccGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctgg                        >  2:1266190/1‑139 (MQ=255)
                                                                                                           cGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                               >  2:311364/1‑67 (MQ=255)
                                                                                                           cGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                               <  1:311364/67‑1 (MQ=255)
                                                                                                           cGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                               <  1:704399/67‑1 (MQ=255)
                                                                                                           cGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                               >  2:704399/1‑67 (MQ=255)
                                                                                                           cGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctggc                       >  2:664084/1‑139 (MQ=255)
                                                                                                           cGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctggc                       >  1:1273561/1‑139 (MQ=255)
                                                                                                                   aTAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGAGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCg               <  1:354328/139‑1 (MQ=255)
                                                                                                                        cTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                               >  2:400930/1‑54 (MQ=255)
                                                                                                                        cTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                               <  1:400930/54‑1 (MQ=255)
                                                                                                                          gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                            <  1:1143020/119‑1 (MQ=255)
                                                                                                                          gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                            >  2:1143020/1‑119 (MQ=255)
                                                                                                                          gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGcc        >  2:1190303/1‑139 (MQ=255)
                                                                                                                          gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGcc        >  1:410568/1‑139 (MQ=255)
                                                                                                                           aaGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCt       >  2:669037/1‑139 (MQ=255)
                                                                                                                             gCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCggg                                               <  2:621689/97‑1 (MQ=255)
                                                                                                                             gCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCggg                                               >  1:621689/1‑97 (MQ=255)
                                                                                                                                ttCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTg             <  1:1136067/128‑1 (MQ=255)
                                                                                                                                ttCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTg             >  2:1136067/1‑128 (MQ=255)
                                                                                                                                ttCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCg  >  1:105365/1‑139 (MQ=255)
                                                                                                                                ttCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGACGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCg  <  1:1190303/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
GTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCG  >  NC_000913/3767045‑3767311

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: