Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 577,266 Δ1 bp 5.7% intergenic (‑441/‑132) 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,2660A. 5.7% 511.8 / 21.8 140intergenic (‑441/‑132)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 A (67/65);  new base . (4/4);  total (71/69)
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.92e-01

Only 100 of 141 total aligned reads displayed.

GTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAA  >  NC_000913/577122‑577414
                                                                                                                                                 |                                                                                                                                                    
gttgAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGATAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCgg                                                                                                                                                  >  1:552351/1‑149 (MQ=255)
      aTTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGtttaa                                                                                                                                            <  1:1286216/149‑3 (MQ=255)
            ttATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGtttaa                                                                                                                                            <  3:11963/143‑3 (MQ=255)
                       tAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCg                                                                                                                           >  3:683661/1‑149 (MQ=255)
                        aGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGa                                                                                                                          <  3:1590935/149‑1 (MQ=255)
                               ttttCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCAtttttt                                                                                                                                   >  4:485664/1‑133 (MQ=255)
                               ttttCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCAtttttt                                                                                                                                   <  3:485664/133‑1 (MQ=255)
                               ttttCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGc                                                                                                                   >  3:1415822/1‑149 (MQ=255)
                                 ttCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                                     >  3:356227/1‑145 (MQ=255)
                                 ttCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                                     <  4:356227/145‑1 (MQ=255)
                                  tCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGa                                                                                                                <  4:1415822/149‑1 (MQ=255)
                                     aTACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTcc                                                                                                                                          >  3:224249/1‑120 (MQ=255)
                                     aTACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTcc                                                                                                                                          <  4:224249/120‑1 (MQ=255)
                                        cTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAAccc                                                                                                          <  2:552351/149‑1 (MQ=255)
                                                cATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGc                                                                                                                                                    >  3:367244/1‑99 (MQ=255)
                                                cATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGc                                                                                                                                                    <  4:367244/99‑1 (MQ=255)
                                                 aTTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGttt                                                                                                                                              >  1:354930/1‑104 (MQ=255)
                                                     tctcAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGtctc                                                                                                                            <  4:1393582/118‑1 (MQ=255)
                                                     tctcAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTa                                                                                             >  3:1197070/1‑149 (MQ=255)
                                                      ctcAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTAt                                                                                            >  3:1421319/1‑149 (MQ=255)
                                                       tcAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGttt                                                                                                                                              >  3:745398/1‑98 (MQ=255)
                                                       tcAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGttt                                                                                                                                              <  4:745398/98‑1 (MQ=255)
                                                         aaaTATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCa                                                                                         >  3:1892800/1‑149 (MQ=255)
                                                         aaaTATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGGGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCa                                                                                         >  1:715545/1‑149 (MQ=255)
                                                            tatCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTg                                                                                      >  4:143432/1‑149 (MQ=255)
                                                             atCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGa                                                                                     <  3:143432/149‑1 (MQ=255)
                                                              tCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGAc                                                                                    >  1:835873/1‑149 (MQ=255)
                                                                aTTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGtttaa                                                                                                                                            <  1:1861475/91‑3 (MQ=38)
                                                                aTTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGtttaa                                                                                                                                            >  2:1861475/1‑89 (MQ=38)
                                                                aTTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGa                                                                                                 <  1:317474/134‑1 (MQ=255)
                                                                aTTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTc                                                                                  >  3:1392104/1‑149 (MQ=255)
                                                                 tttAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGa                                                                                                                                 <  2:1809341/101‑1 (MQ=255)
                                                                 tttAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGa                                                                                                                                 >  1:1809341/1‑101 (MQ=255)
                                                                   tAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAAc                                                                                                            <  1:1060900/120‑1 (MQ=255)
                                                                   tAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAAc                                                                                                            >  2:1060900/1‑120 (MQ=255)
                                                                     gTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGAt                                                                                                               >  2:1785802/1‑115 (MQ=255)
                                                                      tGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTc                                                                            >  1:307623/1‑149 (MQ=255)
                                                                        gTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTc                                                                                          >  4:107656/1‑133 (MQ=255)
                                                                        gTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTc                                                                                          <  3:107656/133‑1 (MQ=255)
                                                                         tttACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGt                                                                         <  4:1392104/149‑1 (MQ=255)
                                                                         tttACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGt                                                                         <  4:683661/149‑1 (MQ=255)
                                                                              gtacgtaAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGt                                                                    <  4:1892800/149‑1 (MQ=255)
                                                                                acgtaAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGttt                                                                  <  3:1380972/149‑1 (MQ=255)
                                                                                  gtaAAAAATTGCTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGa                                                                               >  1:992855/1‑134 (MQ=255)
                                                                                    aaaaaaTTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgta                                                              <  2:307623/149‑1 (MQ=255)
                                                                                      aaaaTTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                                     >  3:1581744/1‑92 (MQ=255)
                                                                                       aaaTTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGATTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGCTTCgtagta                                                           >  3:936967/1‑149 (MQ=255)
                                                                                         aTTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCa                                                                                                                                         <  2:281888/69‑1 (MQ=255)
                                                                                          ttGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                                     <  3:1587324/88‑1 (MQ=255)
                                                                                          ttGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                                     >  4:1587324/1‑88 (MQ=255)
                                                                                          ttGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTc                                                        >  1:331107/1‑149 (MQ=255)
                                                                                          ttGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTc                                                        >  1:1507651/1‑149 (MQ=255)
                                                                                           tGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTcc                                                       <  2:4906/149‑1 (MQ=255)
                                                                                           tGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTcc                                                       <  2:715545/149‑1 (MQ=255)
                                                                                           tGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTcc                                                       >  1:690889/1‑149 (MQ=255)
                                                                                            ggTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTccc                                                      <  4:936967/149‑1 (MQ=255)
                                                                                             gTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAAccc                                                                                                          >  3:1618398/1‑96 (MQ=255)
                                                                                             gTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAAccc                                                                                                          <  4:1618398/96‑1 (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)
                                                                                                      ttAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTg                                                                                      <  4:533287/107‑1 (MQ=255)
                                                                                                      ttAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTaca                                            >  2:478738/1‑149 (MQ=255)
                                                                                                      ttAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTaca                                            >  3:946999/1‑149 (MQ=255)
                                                                                                      ttAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTaca                                            >  1:16162/1‑149 (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                                                           <  4:1278344/127‑1 (MQ=255)
                                                                                                              ggTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCg                                                                                                                           >  2:1624529/1‑62 (MQ=255)
                                                                                                                ttACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgca                                                              <  1:1028726/121‑3 (MQ=255)
                                                                                                                ttACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgca                                                              >  2:1028726/1‑119 (MQ=255)
                                                                                                                 tACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGc                                       <  4:1525134/143‑1 (MQ=255)
                                                                                                                 tACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTc                                 <  2:835873/149‑1 (MQ=255)
                                                                                                                  aCTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGa                                                                                                 <  1:1642943/84‑1 (MQ=255)
                                                                                                                  aCTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGa                                                                                                 >  2:1642943/1‑84 (MQ=255)
                                                                                                                    ttCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTAt                                                          >  3:880453/1‑121 (MQ=255)
                                                                                                                    ttCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTAt                                                          <  4:880453/121‑1 (MQ=255)
                                                                                                                    ttCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGc                              >  1:266733/1‑149 (MQ=255)
                                                                                                                    ttCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGc                              >  2:398346/1‑149 (MQ=255)
                                                                                                                      cGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgtag                                                             <  4:1560089/116‑1 (MQ=255)
                                                                                                                      cGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgtag                                                             >  3:551058/1‑116 (MQ=255)
                                                                                                                      cGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGa                            <  3:334850/149‑1 (MQ=255)
                                                                                                                        tCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAgg                          >  3:500664/1‑149 (MQ=255)
                                                                                                                          acacaGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTaca                                            <  4:775246/129‑1 (MQ=255)
                                                                                                                           cacaGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTgg                       <  3:1560053/149‑1 (MQ=255)
                                                                                                                            acaGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTggg                      <  1:478738/149‑1 (MQ=255)
                                                                                                                             caGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTgggg                     <  2:565155/149‑1 (MQ=255)
                                                                                                                               gCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTacac                                           >  1:1722584/1‑125 (MQ=255)
                                                                                                                                cTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCg                                        <  1:1057292/127‑1 (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)
                                                                                                                                        cgccgTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGc                                       <  2:1605468/120‑1 (MQ=255)
                                                                                                                                        cgccgTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGc                                       >  1:1605468/1‑120 (MQ=255)
                                                                                                                                          ccgTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCAt        <  4:1197070/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAA  >  NC_000913/577122‑577414

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: