Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,997,650 T→A 5.8% intergenic (‑50/+42) 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,6500TA 5.8% 249.6 / 4.6 69intergenic (‑50/+42)insC1/ygeQIS2 repressor TnpA/pseudogene, glycosyl hydrolase family 15, part of T3SS PAI ETT2 remnant
Reads supporting (aligned to +/- strand):  ref base T (32/33);  new base A (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

CCTTCCTGGTATTGCTTACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTT  >  NC_000913/2997434‑2997850
                                                                                                                                                                                                                        |                                                                                                                                                                                                        
ccttcctGGTATTGCTTACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTTTCTGGAGATTNNNGGGGCAAGTCTAAAATTGCTTTTttcttt                                                                                                                                                >  3:148770/1‑274 (MQ=12)
                                   gctggctgCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaat                                                                                                              <  1:148770/274‑1 (MQ=255)
                                                                                      ttCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGTCTTCTCTGTAATGTTTATCTTTTTCCa                                                            >  1:118206/1‑273 (MQ=255)
                                                                                                                      ccTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTATAAACTTGct                                                                                                                                                         <  7:651971/148‑1 (MQ=18)
                                                                                                                          tgtgGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAAc                      >  4:496335/1‑275 (MQ=255)
                                                                                                                                      tCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCa          <  2:496336/275‑1 (MQ=255)
                                                                                                                                         gCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGt                    >  2:558/1‑262 (MQ=255)
                                                                                                                                             ttCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTAtt                <  4:558/262‑1 (MQ=255)
                                                                                                                                              tttcCGGTACTAATACATCATTTATCTGTTCTACAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTANGTGATATTGGTTGTCTGGAGATTCAGGGTGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCttt  <  3:118206/273‑1 (MQ=255)
                                                                                                                                              tctcCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCCGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCt              >  2:306643/1‑263 (MQ=255)
                                                                                                                                                 ccGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTgg                                                                                                                                             >  8:532617/1‑133 (MQ=255)
                                                                                                                                                 ccGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTgg                                                                                                                                             <  7:532617/133‑1 (MQ=255)
                                                                                                                                                  cGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCCGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCa          <  4:306642/263‑1 (MQ=255)
                                                                                                                                                        tAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGt                                                                                                                       >  5:1055747/1‑148 (MQ=255)
                                                                                                                                                                  atcatcTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGAtt                                                                                                                          >  6:491605/1‑135 (MQ=255)
                                                                                                                                                                  atcatcTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGAtt                                                                                                                          <  5:491605/135‑1 (MQ=255)
                                                                                                                                                                  atcatcTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaata                                                                                                             <  6:1055747/148‑1 (MQ=255)
                                                                                                                                                                    catcTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCtt                                                                                                                                   <  7:961045/124‑1 (MQ=255)
                                                                                                                                                                    catcTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCtt                                                                                                                                   >  8:961045/1‑124 (MQ=255)
                                                                                                                                                                          ttCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTACCTTTGAAAGATTTg                                                                                                                        >  2:35990/1‑129 (MQ=255)
                                                                                                                                                                              ccAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCtt                                                                                                                                   >  7:235669/1‑114 (MQ=255)
                                                                                                                                                                              ccAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCtt                                                                                                                                   <  8:235669/114‑1 (MQ=255)
                                                                                                                                                                              ccAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTACCTTTGAAAGATTTGTAgg                                                                                                                    <  4:35989/129‑1 (MQ=255)
                                                                                                                                                                                     cTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGa                                                                                                                   <  7:536019/123‑1 (MQ=255)
                                                                                                                                                                                     cTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGa                                                                                                                   >  8:536019/1‑123 (MQ=255)
                                                                                                                                                                                       aGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaa                                                                                                               >  7:915237/1‑125 (MQ=255)
                                                                                                                                                                                       aGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaa                                                                                                               <  8:915237/125‑1 (MQ=255)
                                                                                                                                                                                         tCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTttcttc                                                                                                                                                >  6:182369/1‑90 (MQ=255)
                                                                                                                                                                                         tCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTttcttc                                                                                                                                                <  5:182369/90‑1 (MQ=255)
                                                                                                                                                                                         tCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTa                                                                                                                      <  8:584935/116‑1 (MQ=255)
                                                                                                                                                                                         tCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTa                                                                                                                      >  7:584935/1‑116 (MQ=255)
                                                                                                                                                                                          cTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTttcttc                                                                                                                                                >  5:864963/1‑89 (MQ=255)
                                                                                                                                                                                          cTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTttcttc                                                                                                                                                <  6:864963/89‑1 (MQ=255)
                                                                                                                                                                                          cTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGgtg                                                                                                                                           >  8:390605/1‑94 (MQ=255)
                                                                                                                                                                                          cTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGgtg                                                                                                                                           <  7:390605/94‑1 (MQ=255)
                                                                                                                                                                                            aaaaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCttt                                                                                                                                                     >  8:14004/1‑82 (MQ=255)
                                                                                                                                                                                            aaaaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCttt                                                                                                                                                     <  7:14004/82‑1 (MQ=255)
                                                                                                                                                                                            aaaaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAAc                                                                                                           >  6:288983/1‑124 (MQ=255)
                                                                                                                                                                                            aaaaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAAc                                                                                                           <  5:288983/124‑1 (MQ=255)
                                                                                                                                                                                             aaaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTa                                                                                                       <  6:951675/127‑1 (MQ=255)
                                                                                                                                                                                             aaaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTa                                                                                                       >  5:951675/1‑127 (MQ=255)
                                                                                                                                                                                                   aGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTttcttc                                                                                                                                                <  7:416308/80‑1 (MQ=255)
                                                                                                                                                                                                   aGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTttcttc                                                                                                                                                >  6:410904/1‑80 (MQ=255)
                                                                                                                                                                                                   aGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTttcttc                                                                                                                                                >  8:416308/1‑80 (MQ=255)
                                                                                                                                                                                                   aGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTttcttc                                                                                                                                                <  5:410904/80‑1 (MQ=255)
                                                                                                                                                                                                   aGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTttctt                                                                                                                                                 <  8:94240/79‑1 (MQ=255)
                                                                                                                                                                                                   aGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTttctt                                                                                                                                                 >  7:94240/1‑79 (MQ=255)
                                                                                                                                                                                                    gTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGaa                                                                                                                               >  8:438896/1‑96 (MQ=255)
                                                                                                                                                                                                    gTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGaa                                                                                                                               <  7:438896/96‑1 (MQ=255)
                                                                                                                                                                                                      aTTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCt                                                                                                                                    >  5:644208/1‑89 (MQ=255)
                                                                                                                                                                                                      aTTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCt                                                                                                                                    <  6:644208/89‑1 (MQ=255)
                                                                                                                                                                                                        tAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGa                                                                                                                   <  8:348861/104‑1 (MQ=255)
                                                                                                                                                                                                        tAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGa                                                                                                                   >  7:348861/1‑104 (MQ=255)
                                                                                                                                                                                                          aGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATc                                                                                         >  4:301478/1‑128 (MQ=255)
                                                                                                                                                                                                              tATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGa                                                                                                                                >  8:978422/1‑85 (MQ=255)
                                                                                                                                                                                                              tATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGa                                                                                                                                <  7:978422/85‑1 (MQ=255)
                                                                                                                                                                                                              tATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCt                                                                                     <  2:301479/128‑1 (MQ=255)
                                                                                                                                                                                                              tATCACTTAAATAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTtct                                                                                  <  7:790544/131‑1 (MQ=255)
                                                                                                                                                                                                              tATCACTTAAATAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTtct                                                                                  >  8:790544/1‑131 (MQ=255)
                                                                                                                                                                                                               aTCACTTAAATAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATc                                                                                         <  7:758678/123‑1 (MQ=255)
                                                                                                                                                                                                               aTCACTTAAATAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATc                                                                                         >  8:758678/1‑123 (MQ=255)
                                                                                                                                                                                                                   cTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTa                                                                                                       >  8:383153/1‑105 (MQ=255)
                                                                                                                                                                                                                   cTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTa                                                                                                       <  7:383153/105‑1 (MQ=255)
                                                                                                                                                                                                                   cTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAg                                                                                                >  7:753299/1‑112 (MQ=255)
                                                                                                                                                                                                                   cTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAg                                                                                                <  8:753299/112‑1 (MQ=255)
                                                                                                                                                                                                                   cTTATTTAAGTGATATTAGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGttt                                                                       >  8:1096160/1‑137 (MQ=255)
                                                                                                                                                                                                                   cTTATTTAAGTGATATTAGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGttt                                                                       <  7:1096160/137‑1 (MQ=255)
                                                                                                                                                                                                                     tATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTa                                                                                                                                       <  7:1014953/71‑1 (MQ=255)
                                                                                                                                                                                                                     tATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTa                                                                                                                                       >  8:1014953/1‑71 (MQ=255)
                                                                                                                                                                                                                        ttAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATc                                                                                                                                     <  6:1068193/70‑1 (MQ=255)
                                                                                                                                                                                                                        ttAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATc                                                                                                                                     >  5:1068193/1‑70 (MQ=255)
                                                                                                                                                                                                                        ttAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataa                                                                                                                <  8:675108/91‑1 (MQ=255)
                                                                                                                                                                                                                        ttAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataa                                                                                                                >  7:675108/1‑91 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                        
CCTTCCTGGTATTGCTTACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTT  >  NC_000913/2997434‑2997850

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑