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.

GTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAA  >  NC_000913/577122‑577414
                                                                                                                                                 |                                                                                                                                                    
gttgAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGATAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCgg                                                                                                                                                  >  1:552351/1‑149 (MQ=255)
      aTTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACG‑GCGGGtttaa                                                                                                                                            <  1:1286216/149‑3 (MQ=255)
            ttATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACG‑GCGGGtttaa                                                                                                                                            <  3:11963/143‑3 (MQ=255)
                       tAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCg                                                                                                                           >  3:683661/1‑149 (MQ=255)
                        aGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGa                                                                                                                          <  3:1590935/149‑1 (MQ=255)
                               ttttCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCAtttttt                                                                                                                                   >  4:485664/1‑133 (MQ=255)
                               ttttCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCAtttttt                                                                                                                                   <  3:485664/133‑1 (MQ=255)
                               ttttCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGc                                                                                                                   >  3:1415822/1‑149 (MQ=255)
                                 ttCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                                     >  3:356227/1‑145 (MQ=255)
                                 ttCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                                     <  4:356227/145‑1 (MQ=255)
                                  tCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGa                                                                                                                <  4:1415822/149‑1 (MQ=255)
                                     aTACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTcc                                                                                                                                          >  3:224249/1‑120 (MQ=255)
                                     aTACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTcc                                                                                                                                          <  4:224249/120‑1 (MQ=255)
                                        cTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAAccc                                                                                                          <  2:552351/149‑1 (MQ=255)
                                                cATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGc                                                                                                                                                    >  3:367244/1‑99 (MQ=255)
                                                cATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGc                                                                                                                                                    <  4:367244/99‑1 (MQ=255)
                                                 aTTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACG‑GCGGGttt                                                                                                                                              >  1:354930/1‑104 (MQ=255)
                                                     tctcAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGtctc                                                                                                                            <  4:1393582/118‑1 (MQ=255)
                                                     tctcAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTa                                                                                             >  3:1197070/1‑149 (MQ=255)
                                                      ctcAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTAt                                                                                            >  3:1421319/1‑149 (MQ=255)
                                                       tcAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACG‑GCGGGttt                                                                                                                                              >  3:745398/1‑98 (MQ=255)
                                                       tcAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACG‑GCGGGttt                                                                                                                                              <  4:745398/98‑1 (MQ=255)
                                                         aaaTATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCa                                                                                         >  3:1892800/1‑149 (MQ=255)
                                                         aaaTATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGGGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCa                                                                                         >  1:715545/1‑149 (MQ=255)
                                                            tatCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTg                                                                                      >  4:143432/1‑149 (MQ=255)
                                                             atCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGa                                                                                     <  3:143432/149‑1 (MQ=255)
                                                              tCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGAc                                                                                    >  1:835873/1‑149 (MQ=255)
                                                                aTTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACG‑GCGGGtttaa                                                                                                                                            <  1:1861475/91‑3 (MQ=38)
                                                                aTTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACG‑GCGGGtttaa                                                                                                                                            >  2:1861475/1‑89 (MQ=38)
                                                                aTTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGa                                                                                                 <  1:317474/134‑1 (MQ=255)
                                                                aTTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTc                                                                                  >  3:1392104/1‑149 (MQ=255)
                                                                 tttAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGa                                                                                                                                 <  2:1809341/101‑1 (MQ=255)
                                                                 tttAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGa                                                                                                                                 >  1:1809341/1‑101 (MQ=255)
                                                                   tAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAAc                                                                                                            <  1:1060900/120‑1 (MQ=255)
                                                                   tAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAAc                                                                                                            >  2:1060900/1‑120 (MQ=255)
                                                                     gTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGAt                                                                                                               >  2:1785802/1‑115 (MQ=255)
                                                                      tGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTc                                                                            >  1:307623/1‑149 (MQ=255)
                                                                        gTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTc                                                                                          >  4:107656/1‑133 (MQ=255)
                                                                        gTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTc                                                                                          <  3:107656/133‑1 (MQ=255)
                                                                         tttACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGt                                                                         <  4:1392104/149‑1 (MQ=255)
                                                                         tttACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGt                                                                         <  4:683661/149‑1 (MQ=255)
                                                                              gtacgtaAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGt                                                                    <  4:1892800/149‑1 (MQ=255)
                                                                                acgtaAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGttt                                                                  <  3:1380972/149‑1 (MQ=255)
                                                                                  gtaAAAAATTGCTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGa                                                                               >  1:992855/1‑134 (MQ=255)
                                                                                    aaaaaaTTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgta                                                              <  2:307623/149‑1 (MQ=255)
                                                                                      aaaaTTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                                     >  3:1581744/1‑92 (MQ=255)
                                                                                       aaaTTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGATTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGCTTCgtagta                                                           >  3:936967/1‑149 (MQ=255)
                                                                                         aTTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCa                                                                                                                                         <  2:281888/69‑1 (MQ=255)
                                                                                          ttGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                                     <  3:1587324/88‑1 (MQ=255)
                                                                                          ttGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                                     >  4:1587324/1‑88 (MQ=255)
                                                                                          ttGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTc                                                        >  1:331107/1‑149 (MQ=255)
                                                                                          ttGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTc                                                        >  1:1507651/1‑149 (MQ=255)
                                                                                           tGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTcc                                                       <  2:4906/149‑1 (MQ=255)
                                                                                           tGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTcc                                                       <  2:715545/149‑1 (MQ=255)
                                                                                           tGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTcc                                                       >  1:690889/1‑149 (MQ=255)
                                                                                            ggTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTccc                                                      <  4:936967/149‑1 (MQ=255)
                                                                                             gTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAAccc                                                                                                          >  3:1618398/1‑96 (MQ=255)
                                                                                             gTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAAccc                                                                                                          <  4:1618398/96‑1 (MQ=255)
                                                                                             gTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTg                                                                                      <  1:1623223/116‑1 (MQ=255)
                                                                                                aTGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCAtttt                                                                                                                                     <  3:815738/66‑1 (MQ=255)
                                                                                                aTGCTGTTAAGAGTGGTTACTCCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCAtttt                                                                                                                                     >  4:815738/1‑66 (MQ=255)
                                                                                                 tGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAAt                                                 >  4:334850/1‑149 (MQ=255)
                                                                                                     gTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTc                                                                                  >  3:1547978/1‑112 (MQ=255)
                                                                                                      ttAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTg                                                                                      <  4:533287/107‑1 (MQ=255)
                                                                                                      ttAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTaca                                            >  2:478738/1‑149 (MQ=255)
                                                                                                      ttAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTaca                                            >  3:946999/1‑149 (MQ=255)
                                                                                                      ttAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTaca                                            >  1:16162/1‑149 (MQ=255)
                                                                                                         agagTAGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGc                                                                                                                   >  1:538079/1‑75 (MQ=255)
                                                                                                            gTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGa                                                                               >  4:162593/1‑108 (MQ=255)
                                                                                                            gTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGa                                                                               <  3:162593/108‑1 (MQ=255)
                                                                                                            gTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCt                                      >  4:1560053/1‑149 (MQ=255)
                                                                                                             tGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgtagta                                                           <  4:1278344/127‑1 (MQ=255)
                                                                                                              ggTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCg                                                                                                                           >  2:1624529/1‑62 (MQ=255)
                                                                                                                ttACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgca                                                              <  1:1028726/121‑3 (MQ=255)
                                                                                                                ttACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgca                                                              >  2:1028726/1‑119 (MQ=255)
                                                                                                                 tACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGc                                       <  4:1525134/143‑1 (MQ=255)
                                                                                                                 tACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTc                                 <  2:835873/149‑1 (MQ=255)
                                                                                                                  aCTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGa                                                                                                 <  1:1642943/84‑1 (MQ=255)
                                                                                                                  aCTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGa                                                                                                 >  2:1642943/1‑84 (MQ=255)
                                                                                                                    ttCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTAt                                                          >  3:880453/1‑121 (MQ=255)
                                                                                                                    ttCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTAt                                                          <  4:880453/121‑1 (MQ=255)
                                                                                                                    ttCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGc                              >  1:266733/1‑149 (MQ=255)
                                                                                                                    ttCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGc                              >  2:398346/1‑149 (MQ=255)
                                                                                                                      cGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgtag                                                             <  4:1560089/116‑1 (MQ=255)
                                                                                                                      cGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgtag                                                             >  3:551058/1‑116 (MQ=255)
                                                                                                                      cGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGa                            <  3:334850/149‑1 (MQ=255)
                                                                                                                        tCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAgg                          >  3:500664/1‑149 (MQ=255)
                                                                                                                          acacaGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTaca                                            <  4:775246/129‑1 (MQ=255)
                                                                                                                           cacaGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTgg                       <  3:1560053/149‑1 (MQ=255)
                                                                                                                            acaGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTggg                      <  1:478738/149‑1 (MQ=255)
                                                                                                                             caGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTgggg                     <  2:565155/149‑1 (MQ=255)
                                                                                                                               gCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTacac                                           >  1:1722584/1‑125 (MQ=255)
                                                                                                                                cTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCg                                        <  1:1057292/127‑1 (MQ=255)
                                                                                                                                 ttAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCa                                                                                         <  2:1241751/77‑1 (MQ=255)
                                                                                                                                 ttAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTAt                 <  2:16162/149‑1 (MQ=255)
                                                                                                                                 ttAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTAt                 <  2:1550230/149‑1 (MQ=255)
                                                                                                                                       ccgc‑cgTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATc           >  4:939976/1‑149 (MQ=255)
                                                                                                                                        cgc‑cgTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGc                                       <  2:1605468/120‑1 (MQ=255)
                                                                                                                                        cgc‑cgTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGc                                       >  1:1605468/1‑120 (MQ=255)
                                                                                                                                          c‑cgTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCAt        <  4:1197070/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGC‑CGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAA  >  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: ‑