Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 577,260:1 +T 6.0% intergenic (‑435/‑138) nmpC ← / → essD DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage/DLP12 prophage; putative phage lysis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913577,2601.T 6.0% 582.7 / 30.0 134intergenic (‑435/‑138)nmpC/essDDLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage/DLP12 prophage; putative phage lysis protein
Reads supporting (aligned to +/- strand):  ref base . (65/61);  new base T (4/4);  total (69/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 139 total aligned reads displayed.

TTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCAT  >  NC_000913/577116‑577408
                                                                                                                                                 |                                                                                                                                                    
ttATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTc                                                                                                                                                  >  4:1590935/1‑149 (MQ=255)
      gttgAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGATAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCgg                                                                                                                                            >  1:552351/1‑149 (MQ=255)
            aTTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGtttaa                                                                                                                                      <  1:1286216/149‑3 (MQ=255)
                  ttATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGtttaa                                                                                                                                      <  3:11963/143‑3 (MQ=255)
                  ttATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGtttaa                                                                                                                                      >  4:11963/1‑141 (MQ=255)
                                     ttttCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCAtttttt                                                                                                                             >  4:485664/1‑133 (MQ=255)
                                     ttttCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCAtttttt                                                                                                                             <  3:485664/133‑1 (MQ=255)
                                     ttttCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGc                                                                                                             >  3:1415822/1‑149 (MQ=255)
                                      tttCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCt                                                                                                            <  4:788308/149‑1 (MQ=255)
                                       ttCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                               <  4:356227/145‑1 (MQ=255)
                                        tCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGa                                                                                                          <  4:1415822/149‑1 (MQ=255)
                                          aaTACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATa                                                                                                        >  4:1380972/1‑149 (MQ=255)
                                           aTACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTcc                                                                                                                                    <  4:224249/120‑1 (MQ=255)
                                           aTACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTcc                                                                                                                                    >  3:224249/1‑120 (MQ=255)
                                              cTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAAccc                                                                                                    <  2:552351/149‑1 (MQ=255)
                                                   tAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATAc                                                                                               <  1:730002/149‑1 (MQ=255)
                                                      cATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGc                                                                                                                                              >  3:367244/1‑99 (MQ=255)
                                                       aTTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGttt                                                                                                                                        >  1:354930/1‑104 (MQ=255)
                                                        ttCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGa                                                                                           <  3:908966/148‑1 (MQ=255)
                                                           tctcAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGtctc                                                                                                                      >  3:1393582/1‑118 (MQ=255)
                                                            ctcAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTAt                                                                                      >  3:1421319/1‑149 (MQ=255)
                                                             tcAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGttt                                                                                                                                        <  4:745398/98‑1 (MQ=255)
                                                             tcAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGttt                                                                                                                                        >  3:745398/1‑98 (MQ=255)
                                                             tcAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTAtt                                                                                     >  1:322260/1‑149 (MQ=255)
                                                               aaaTATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGGGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCa                                                                                   >  1:715545/1‑149 (MQ=255)
                                                                  tatCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTg                                                                                >  4:143432/1‑149 (MQ=255)
                                                                    tCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGAc                                                                              >  1:835873/1‑149 (MQ=255)
                                                                      aTTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGtttaa                                                                                                                                      <  1:1861475/91‑3 (MQ=38)
                                                                      aTTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGtttaa                                                                                                                                      >  2:1861475/1‑89 (MQ=38)
                                                                      aTTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGa                                                                                           >  2:317474/1‑134 (MQ=255)
                                                                      aTTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGa                                                                                           <  1:317474/134‑1 (MQ=255)
                                                                         tAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAAc                                                                                                      <  1:1060900/120‑1 (MQ=255)
                                                                           gTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGAt                                                                                                         <  1:1785802/115‑1 (MQ=255)
                                                                           gTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGAt                                                                                                         >  2:1785802/1‑115 (MQ=255)
                                                                            tGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTc                                                                      >  1:307623/1‑149 (MQ=255)
                                                                              gTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTc                                                                                    <  3:107656/133‑1 (MQ=255)
                                                                              gTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTc                                                                                    >  4:107656/1‑133 (MQ=255)
                                                                               tttACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGt                                                                   <  4:683661/149‑1 (MQ=255)
                                                                               tttACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGt                                                                   <  4:1392104/149‑1 (MQ=255)
                                                                                    gtacgtaAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGt                                                              <  4:1892800/149‑1 (MQ=255)
                                                                                      acgtaAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGttt                                                            <  3:1380972/149‑1 (MQ=255)
                                                                                        gtaAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGa                                                                         <  2:992855/134‑1 (MQ=255)
                                                                                         taAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTTGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgt                                                         >  1:565155/1‑149 (MQ=255)
                                                                                         taAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgt                                                         >  3:448104/1‑149 (MQ=255)
                                                                                          aaaaaaTTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgta                                                        <  2:307623/149‑1 (MQ=255)
                                                                                            aaaaTTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                               >  3:1581744/1‑92 (MQ=255)
                                                                                             aaaTTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGATTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGCTTCgtagta                                                     >  3:936967/1‑149 (MQ=255)
                                                                                                ttGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                               >  4:1587324/1‑88 (MQ=255)
                                                                                                ttGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                               <  3:1587324/88‑1 (MQ=255)
                                                                                                ttGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTc                                                  >  1:1507651/1‑149 (MQ=255)
                                                                                                 tGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTcc                                                 <  2:4906/149‑1 (MQ=255)
                                                                                                 tGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTcc                                                 <  2:715545/149‑1 (MQ=255)
                                                                                                  ggTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTccc                                                <  4:936967/149‑1 (MQ=255)
                                                                                                   gTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAAccc                                                                                                    >  3:1618398/1‑96 (MQ=255)
                                                                                                   gTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTg                                                                                >  2:1623223/1‑116 (MQ=255)
                                                                                                   gTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTg                                                                                <  1:1623223/116‑1 (MQ=255)
                                                                                                      aTGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCAtttt                                                                                                                               <  3:815738/66‑1 (MQ=255)
                                                                                                      aTGCTGTTAAGAGTGGTTACTCCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCAtttt                                                                                                                               >  4:815738/1‑66 (MQ=255)
                                                                                                       tGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAAt                                           >  4:334850/1‑149 (MQ=255)
                                                                                                           gTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTc                                                                            >  3:1547978/1‑112 (MQ=255)
                                                                                                           gTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTc                                                                            <  4:1547978/112‑1 (MQ=255)
                                                                                                           gTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTac                                       <  2:1507651/149‑1 (MQ=255)
                                                                                                            ttAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTg                                                                                <  4:533287/107‑1 (MQ=255)
                                                                                                            ttAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTg                                                                                >  3:533287/1‑107 (MQ=255)
                                                                                                            ttAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTaca                                      >  1:16162/1‑149 (MQ=255)
                                                                                                               agagTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGc                                                                                                             <  2:538079/75‑1 (MQ=255)
                                                                                                               agagTAGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGc                                                                                                             >  1:538079/1‑75 (MQ=255)
                                                                                                                  gTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGa                                                                         >  4:162593/1‑108 (MQ=255)
                                                                                                                  gTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGa                                                                         <  3:162593/108‑1 (MQ=255)
                                                                                                                  gTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCt                                >  4:1560053/1‑149 (MQ=255)
                                                                                                                   tGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgtagta                                                     >  3:1278344/1‑127 (MQ=255)
                                                                                                                      ttACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgca                                                        <  1:1028726/121‑3 (MQ=255)
                                                                                                                      ttACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgca                                                        >  2:1028726/1‑119 (MQ=255)
                                                                                                                       tACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGc                                 >  3:1525134/1‑143 (MQ=255)
                                                                                                                        aCTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGa                                                                                           <  1:1642943/84‑1 (MQ=255)
                                                                                                                        aCTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGa                                                                                           >  2:1642943/1‑84 (MQ=255)
                                                                                                                        aCTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCt                          <  4:1421319/149‑1 (MQ=255)
                                                                                                                          ttCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTAt                                                    <  4:880453/121‑1 (MQ=255)
                                                                                                                          ttCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGc                        >  1:266733/1‑149 (MQ=255)
                                                                                                                          ttCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGc                        >  2:398346/1‑149 (MQ=255)
                                                                                                                            cGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgtag                                                       >  3:1560089/1‑116 (MQ=255)
                                                                                                                            cGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgtag                                                       >  3:551058/1‑116 (MQ=255)
                                                                                                                            cGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgtag                                                       <  4:1560089/116‑1 (MQ=255)
                                                                                                                            cGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGa                      <  3:334850/149‑1 (MQ=255)
                                                                                                                                acacaGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAAt                                           <  3:227644/124‑1 (MQ=255)
                                                                                                                                acacaGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTaca                                      >  3:775246/1‑129 (MQ=255)
                                                                                                                                  acaGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTggg                <  1:478738/149‑1 (MQ=255)
                                                                                                                                  acaGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTggg                <  2:266733/149‑1 (MQ=255)
                                                                                                                                  acaGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTggg                >  1:1550230/1‑149 (MQ=255)
                                                                                                                                     gCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTacac                                     >  1:1722584/1‑125 (MQ=255)
                                                                                                                                     gCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTacac                                     <  2:1722584/125‑1 (MQ=255)
                                                                                                                                     gCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGtt             <  4:448104/149‑1 (MQ=255)
                                                                                                                                      cTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCg                                  <  1:1057292/127‑1 (MQ=255)
                                                                                                                                      cTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCg                                  >  2:1057292/1‑127 (MQ=255)
                                                                                                                                       ttAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCa                                                                                   <  2:1241751/77‑1 (MQ=255)
                                                                                                                                       ttAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTAt           <  2:16162/149‑1 (MQ=255)
                                                                                                                                       ttAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTAt           <  2:1550230/149‑1 (MQ=255)
                                                                                                                                             ccgccgTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATc     >  4:939976/1‑149 (MQ=255)
                                                                                                                                                ccgTCGAGCGGGTTTTTCCATTTTTTTAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCAt  <  2:690889/149‑1 (MQ=255)
                                                                                                                                                ccgTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCAt  <  4:1197070/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCAT  >  NC_000913/577116‑577408

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