Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,997,689 C→T intergenic (‑89/+3) insC1 ← / ← ygeQ IS2 repressor TnpA/pseudogene, glycosyl hydrolase family 15, part of T3SS PAI ETT2 remnant

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,997,6890CT97.8% 272.7 / ‑2.4 89intergenic (‑89/+3)insC1/ygeQIS2 repressor TnpA/pseudogene, glycosyl hydrolase family 15, part of T3SS PAI ETT2 remnant
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (47/40);  total (48/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCA  >  NC_000913/2997558‑2997821
                                                                                                                                   |                                                                                                                                    
tgGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACtt                                                                                                                               >  1:1275625/1‑139 (MQ=14)
                             aaGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGaa                                                                                                  >  1:185130/1‑139 (MQ=255)
                                   tcaatcaTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTg                                                                                                                              <  2:575469/105‑1 (MQ=17)
                                   tcaatcaTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTg                                                                                                                              >  1:575469/1‑105 (MQ=17)
                                                 tcCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTg                                                                                                                              <  2:133780/91‑1 (MQ=18)
                                                 tcCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTg                                                                                                                              >  1:133780/1‑91 (MQ=18)
                                                 tcCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAAc                                                                              <  2:1275625/139‑1 (MQ=255)
                                                   cAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACAt                                                                            >  2:1389974/1‑139 (MQ=255)
                                                      tGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTa                                                                                                          >  2:1436245/1‑106 (MQ=255)
                                                      tGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTa                                                                                                          <  1:1436245/106‑1 (MQ=255)
                                                        aCTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGt                                                                       >  2:1353956/1‑139 (MQ=255)
                                                           aGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTaaa                                                                    <  1:355092/139‑1 (MQ=255)
                                                           aGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTaaa                                                                    >  2:355092/1‑139 (MQ=255)
                                                            gTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGACTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAg                                                                   <  2:1338761/139‑1 (MQ=255)
                                                             tCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGc                                                                  >  2:730885/1‑139 (MQ=255)
                                                                     cTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCgg                                                          <  1:730885/139‑1 (MQ=255)
                                                                        gTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCtt                                                       >  1:1249587/1‑139 (MQ=255)
                                                                        gTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCtt                                                       >  1:1106764/1‑139 (MQ=255)
                                                                                   aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGAttt                                                                                            <  2:1226376/91‑1 (MQ=255)
                                                                                   aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGAttt                                                                                            >  1:1226376/1‑91 (MQ=255)
                                                                                   aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGt                                            >  1:299151/1‑139 (MQ=255)
                                                                                   aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCgg                                                          <  2:1386685/125‑1 (MQ=255)
                                                                                   aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCgg                                                          >  1:1386685/1‑125 (MQ=255)
                                                                                    tctcTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTtct                                                     <  2:900950/126‑1 (MQ=255)
                                                                                    tCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTtct                                                     >  1:900950/1‑129 (MQ=255)
                                                                                          aTTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCtt                                     >  1:1086522/1‑139 (MQ=255)
                                                                                           tttAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCttt                                    >  1:527759/1‑139 (MQ=255)
                                                                                            ttAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTg                                                                                           <  2:709005/83‑1 (MQ=255)
                                                                                            ttAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTg                                                                                           >  1:709005/1‑83 (MQ=255)
                                                                                              aaGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGtt                                                                                                            <  1:441023/64‑1 (MQ=255)
                                                                                              aaGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGtt                                                                                                            >  2:441023/1‑64 (MQ=255)
                                                                                              aaGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaa                                                                                  <  2:186211/90‑1 (MQ=255)
                                                                                              aaGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaa                                                                                  >  1:186211/1‑90 (MQ=255)
                                                                                              aaGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTc                                 >  2:849737/1‑139 (MQ=255)
                                                                                                gTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTCAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaa                                                                                  <  2:372308/88‑1 (MQ=255)
                                                                                                gTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaat                                                                                 <  2:228149/89‑1 (MQ=255)
                                                                                                gTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaat                                                                                 >  1:228149/1‑89 (MQ=255)
                                                                                                gTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaa                                                                                  >  1:372308/1‑88 (MQ=255)
                                                                                                  gATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGa                             <  2:527759/139‑1 (MQ=255)
                                                                                                  gATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGa                             <  1:1353956/139‑1 (MQ=255)
                                                                                                     atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTa                                                                                                          >  1:688806/1‑59 (MQ=255)
                                                                                                     atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTa                                                                                                          <  2:688806/59‑1 (MQ=255)
                                                                                                     atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                       >  2:522007/1‑78 (MQ=255)
                                                                                                     atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                       <  1:522007/78‑1 (MQ=255)
                                                                                                     atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaa                                                                                  <  1:700496/83‑1 (MQ=255)
                                                                                                     atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaa                                                                                  >  2:700496/1‑83 (MQ=255)
                                                                                                     atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTATCTGTAATGTTTATCTTTTTCCAGATaa                          >  1:905603/1‑139 (MQ=255)
                                                                                                         gTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTg                      >  1:1196954/1‑139 (MQ=255)
                                                                                                          ttGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGgtgt                                                                                                             >  2:542499/1‑51 (MQ=255)
                                                                                                          ttGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGgtgt                                                                                                             <  1:542499/51‑1 (MQ=255)
                                                                                                           tGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGctc                                                                                                                           >  1:285850/1‑36 (MQ=255)
                                                                                                           tGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGctc                                                                                                                           <  2:285850/36‑1 (MQ=255)
                                                                                                           tGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGt                                                 <  2:856739/110‑1 (MQ=255)
                                                                                                           tGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGt                                                 >  1:856739/1‑110 (MQ=255)
                                                                                                           tGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTggg                    <  1:849737/139‑1 (MQ=255)
                                                                                                            gTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCa                                                              <  2:411445/96‑1 (MQ=255)
                                                                                                            gTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCa                                                              >  1:411445/1‑96 (MQ=255)
                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTGAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGa                  >  1:1399485/1‑139 (MQ=255)
                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGaaa                                                                                                 >  1:853258/1‑60 (MQ=255)
                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGaaa                                                                                                 <  2:853258/60‑1 (MQ=255)
                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTg                                                                                           >  1:964131/1‑66 (MQ=255)
                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTg                                                                                           <  2:964131/66‑1 (MQ=255)
                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                       >  1:1315522/1‑70 (MQ=255)
                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                       >  1:560138/1‑70 (MQ=255)
                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                       <  2:560138/70‑1 (MQ=255)
                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                       <  2:1315522/70‑1 (MQ=255)
                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCAt                                                             <  1:1127723/96‑1 (MQ=255)
                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCAt                                                             >  2:1127723/1‑96 (MQ=255)
                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGa                  >  1:1217646/1‑139 (MQ=255)
                                                                                                               tGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGgtgt                                                                                                             <  2:449557/46‑1 (MQ=255)
                                                                                                               tGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGgtgt                                                                                                             >  1:449557/1‑46 (MQ=255)
                                                                                                               tGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGAtt                >  2:579542/1‑139 (MQ=255)
                                                                                                               tGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGAtt                <  2:1217646/139‑1 (MQ=255)
                                                                                                                  agaTTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaata                                                                                >  1:1332421/1‑72 (MQ=255)
                                                                                                                  agaTTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaata                                                                                <  2:1332421/72‑1 (MQ=255)
                                                                                                                     ttCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTcccc          <  2:299151/139‑1 (MQ=255)
                                                                                                                      tCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCt                                                        >  2:556948/1‑92 (MQ=255)
                                                                                                                      tCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCt                                                        <  1:556948/92‑1 (MQ=255)
                                                                                                                      tCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTtct                                                     >  1:505879/1‑95 (MQ=255)
                                                                                                                      tCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTtct                                                     <  2:505879/95‑1 (MQ=255)
                                                                                                                       cAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCt                                                        <  1:1066942/91‑1 (MQ=255)
                                                                                                                       cAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTa                                                <  1:476732/99‑1 (MQ=255)
                                                                                                                       cAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTa                                                >  2:476732/1‑99 (MQ=255)
                                                                                                                       cAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGAt                            <  2:193353/119‑1 (MQ=255)
                                                                                                                       cAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGAt                            >  1:193353/1‑119 (MQ=255)
                                                                                                                       cAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTATGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCt                                                        >  2:1066942/1‑91 (MQ=255)
                                                                                                                            ggCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGa                             <  2:630510/113‑1 (MQ=255)
                                                                                                                            ggCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGa                             >  1:630510/1‑113 (MQ=255)
                                                                                                                             gCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCa  <  2:1086522/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                    
TGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCA  >  NC_000913/2997558‑2997821

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