Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 577,263 T→A 5.6% intergenic (‑438/‑135) 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,2630TA 5.6% 494.5 / 11.2 139intergenic (‑438/‑135)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 T (67/64);  new base A (4/4);  total (71/68)
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 141 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                                                                                                                                      >  4:11963/1‑141 (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)
                                     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)
                                          aaTACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATa                                                                                                        >  4:1380972/1‑149 (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)
                                                        ttCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGa                                                                                           <  3:908966/148‑1 (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)
                                                               aaaTATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGGGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCa                                                                                   >  1:715545/1‑149 (MQ=255)
                                                                   atCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGa                                                                               <  3:143432/149‑1 (MQ=255)
                                                                    tCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGAc                                                                              >  1:835873/1‑149 (MQ=255)
                                                                      aTTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGtttaa                                                                                                                                      >  2:1861475/1‑89 (MQ=38)
                                                                      aTTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCTCGACGGCGGGtttaa                                                                                                                                      <  1:1861475/91‑3 (MQ=38)
                                                                      aTTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGa                                                                                           <  1:317474/134‑1 (MQ=255)
                                                                      aTTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTc                                                                            >  3:1392104/1‑149 (MQ=255)
                                                                       tttAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGa                                                                                                                           >  1:1809341/1‑101 (MQ=255)
                                                                       tttAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGa                                                                                                                           <  2:1809341/101‑1 (MQ=255)
                                                                         tAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAAc                                                                                                      >  2:1060900/1‑120 (MQ=255)
                                                                         tAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAAc                                                                                                      <  1:1060900/120‑1 (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: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)
                                                                                         taAAAAATTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgt                                                         >  3:448104/1‑149 (MQ=255)
                                                                                          aaaaaaTTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgta                                                        <  2:307623/149‑1 (MQ=255)
                                                                                            aaaaTTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                               >  3:1581744/1‑92 (MQ=255)
                                                                                            aaaaTTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                               <  4:1581744/92‑1 (MQ=255)
                                                                                               aTTGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCa                                                                                                                                   <  2:281888/69‑1 (MQ=255)
                                                                                                ttGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTa                                                                                                               <  3:1587324/88‑1 (MQ=255)
                                                                                                ttGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTc                                                  >  1:331107/1‑149 (MQ=255)
                                                                                                ttGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTc                                                  >  1:1507651/1‑149 (MQ=255)
                                                                                                 tGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTcc                                                 >  1:690889/1‑149 (MQ=255)
                                                                                                 tGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTcc                                                 <  2:715545/149‑1 (MQ=255)
                                                                                                 tGGTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTcc                                                 <  2:4906/149‑1 (MQ=255)
                                                                                                  ggTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTccc                                                <  4:936967/149‑1 (MQ=255)
                                                                                                   gTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAAccc                                                                                                    <  4:1618398/96‑1 (MQ=255)
                                                                                                   gTTATGCTGTTAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTg                                                                                <  1:1623223/116‑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)
                                                                                                            ttAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTg                                                                                <  4:533287/107‑1 (MQ=255)
                                                                                                            ttAAGAGTGGTTACTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTg                                                                                >  3:533287/1‑107 (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)
                                                                                                                   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)
                                                                                                                        aCTTCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCt                          <  4:1421319/149‑1 (MQ=255)
                                                                                                                          ttCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTAt                                                    >  3:880453/1‑121 (MQ=255)
                                                                                                                          ttCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTAt                                                    <  4:880453/121‑1 (MQ=255)
                                                                                                                          ttCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGc                        >  2:398346/1‑149 (MQ=255)
                                                                                                                          ttCGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGc                        >  1:266733/1‑149 (MQ=255)
                                                                                                                            cGTCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCgtag                                                       <  4:1560089/116‑1 (MQ=255)
                                                                                                                              tCACACAGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAgg                    >  3:500664/1‑149 (MQ=255)
                                                                                                                                acacaGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAAt                                           <  3:227644/124‑1 (MQ=255)
                                                                                                                                acacaGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTaca                                      <  4:775246/129‑1 (MQ=255)
                                                                                                                                acacaGCTTAAACCCGCCGTCGAGCGGGTTTTTCCATTTTTTGAGTCTCGATATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTaca                                      >  3:775246/1‑129 (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)
                                                                                                                                                    tCGAGCGGGTTTTTCCATTTTTTGAGTCTCGCTATTAGCTGATAACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTacac                                     <  3:1786372/111‑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: