Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 566,570 A→C 100% pseudogene (195/312 nt) peaD → pseudogene; DLP12 prophage; replication protein family;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913566,5700AC98.2% 7522.5 / 47.5 2595pseudogene (195/312 nt)peaDpseudogene; DLP12 prophage; replication protein family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base C (1333/1205);  total (1362/1233)
Fisher's exact test for biased strand distribution p-value = 8.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 2610 total aligned reads displayed.

TCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑TGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCT  >  NC_000913/566423‑566715
                                                                                                                                                    |                                                                                                                                                 
  cGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATccc                                                                                                                                                 <  6:52767/148‑1 (MQ=255)
      gAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                                                                                                                                            >  1:154018/1‑149 (MQ=255)
         aTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAg                                                                                                                                         >  5:33136/1‑149 (MQ=255)
         aTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAg                                                                                                                                         >  5:198009/1‑149 (MQ=255)
         aTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAg                                                                                                                                         <  4:319360/149‑1 (MQ=255)
             aaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTaa                                                                                                                                      >  8:228611/1‑148 (MQ=255)
              aCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTaaaa                                                                                                                                    >  1:210653/1‑149 (MQ=255)
                  gcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGCTCCCTGAACCAGTAAaacaac                                                                                                                                >  1:101471/1‑149 (MQ=255)
                  gcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAaacaac                                                                                                                                <  4:228882/149‑1 (MQ=255)
                    actACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCa                                                                                                                                          <  4:111607/135‑1 (MQ=255)
                    actACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACtt                                                                                                                              >  3:279727/3‑149 (MQ=255)
                     ctACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTc                                                                                                                             <  3:16851/148‑1 (MQ=255)
                      tACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCAG‑ATGACTGCGAGAATTAACCGTGGTGCGGCGATCCCTGAACCAGTAAAACAACTTcc                                                                                                                            >  5:154417/1‑149 (MQ=255)
                      tACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTcc                                                                                                                            >  2:279764/1‑149 (MQ=255)
                      tACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTcc                                                                                                                            >  6:75961/1‑149 (MQ=255)
                        ctggctggTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTg                                                                                                                          >  3:52493/1‑149 (MQ=255)
                            ctggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCAt                                                                                                                      >  5:65136/2‑149 (MQ=255)
                                    cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCgg                                                                                                                >  3:363046/1‑147 (MQ=255)
                                    cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTa                                                                                                              <  4:209335/149‑1 (MQ=255)
                                    cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAAACGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCgg                                                                                                                >  2:235875/1‑147 (MQ=255)
                                      acctGTATCAGAACATGCGGGCCCATGCGCTTACTGATGCGGAATTCCGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAAGTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGGAgc                                                                                                            >  1:71068/4‑148 (MQ=255)
                                      acctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTGCATATGACTGCGAGAATTCACCGTGGTGAGGCGATCCCGGAACCAGTAAAACAACGTCCGGTCATGGGCGGTAGa                                                                                                            >  6:17866/4‑149 (MQ=255)
                                       cctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGAc                                                                                                           <  6:260262/147‑1 (MQ=255)
                                           tatcAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACct                                                                                                         <  8:271914/144‑1 (MQ=255)
                                             tcAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTa                                                                                                      >  5:113166/2‑148 (MQ=255)
                                             tcAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTa                                                                                                      >  7:285511/2‑148 (MQ=255)
                                               aGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAAt                                                                                                   >  4:291255/1‑149 (MQ=255)
                                                gAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAGCAACTTCCTGTCATGGGCGGTAGCCCTCGAAATc                                                                                                  >  6:344349/1‑149 (MQ=255)
                                                gAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATc                                                                                                  >  7:339252/1‑149 (MQ=255)
                                                gAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCAAATGAGTGCGAGACTTAACCGTGGTGAGGCGATCCCTGAACCAGTAAACCAACTTCCTGTCATGGGCGGTAGACCTCTAACTc                                                                                                  <  6:373036/149‑1 (MQ=255)
                                                  aCATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGt                                                                                                >  8:199511/1‑149 (MQ=255)
                                                  aCATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGt                                                                                                <  8:360978/149‑1 (MQ=255)
                                                   cATGCTGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTg                                                                                               <  7:103137/149‑1 (MQ=255)
                                                   cATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTg                                                                                               <  2:337609/149‑1 (MQ=255)
                                                   cATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTg                                                                                               <  2:258048/149‑1 (MQ=255)
                                                   cATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTg                                                                                               <  2:202987/149‑1 (MQ=255)
                                                     tGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGca                                                                                             <  7:313516/149‑1 (MQ=255)
                                                     tGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGca                                                                                             <  6:116556/149‑1 (MQ=255)
                                                          gCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAg                                                                                                             <  2:309124/128‑1 (MQ=255)
                                                          gCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGc                                                                                        <  4:136941/149‑1 (MQ=255)
                                                          gCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGc                                                                                        >  4:59297/1‑149 (MQ=255)
                                                              aTGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTg                                                                                    <  4:178159/149‑1 (MQ=255)
                                                              aTGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTg                                                                                    >  3:308328/1‑149 (MQ=255)
                                                              aTGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTg                                                                                    <  5:279461/149‑1 (MQ=255)
                                                               tGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCTATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGGAGACCTCTAAATCGTGCACAGGCTCTgg                                                                                   <  1:283767/149‑1 (MQ=255)
                                                               tGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGGCATGGGCGGTAGACCTCTAAATCGTGCACAGGctct                                                                                     <  5:306458/147‑1 (MQ=255)
                                                                  gcttactgATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCg                                                                                                                 <  2:126160/111‑1 (MQ=255)
                                                                   cttactgATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGaa                                                                               <  7:63086/145‑1 (MQ=255)
                                                                    ttactgATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGTGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCg                                                                                 >  5:95239/4‑146 (MQ=255)
                                                                    ttactgATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCg                                                                                 >  3:334140/4‑146 (MQ=255)
                                                                    ttactgATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCg                                                                                 >  2:363712/4‑146 (MQ=255)
                                                                      actgATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGa                                                                             <  1:17773/147‑1 (MQ=255)
                                                                        tgATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATTGGAGTTAGACCTCTAAATCGTGCCCAGGATCTGGCGAAGATCg                                                                          >  2:98823/1‑149 (MQ=255)
                                                                          atgcgGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCg                                                                          >  3:223032/3‑147 (MQ=255)
                                                                          atgcgGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCa                                                                        <  7:258116/147‑1 (MQ=255)
                                                                          atgcgGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCa                                                                        >  1:396490/3‑149 (MQ=255)
                                                                             cgGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGaa                                                                     <  3:89150/149‑1 (MQ=255)
                                                                             cgGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGaa                                                                     <  8:111879/149‑1 (MQ=255)
                                                                              ggAATTACGGCGCAAGGCTGCCGATGAAC‑TGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGaaa                                                                    <  4:356351/149‑1 (MQ=255)
                                                                              ggAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGaaa                                                                    <  8:138462/149‑1 (MQ=255)
                                                                              ggAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGaaa                                                                    <  8:306100/149‑1 (MQ=255)
                                                                              ggAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGa                                                                      >  6:147889/1‑147 (MQ=255)
                                                                               gAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAAt                                                                   <  6:313324/149‑1 (MQ=255)
                                                                                   tACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAAtca                                                                 >  2:101678/1‑145 (MQ=255)
                                                                                                            ttgtccatATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGc                                        >  5:224/8‑148 (MQ=255)
                                                                                                              gtccatATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCaa                                      <  2:389334/143‑1 (MQ=255)
                                                                                                               tccatATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAg                                     <  7:230024/144‑1 (MQ=255)
                                                                                                                ccatttGACTGCGAGAATTAACCGTTGTGAGGCGATCCCTGAACCAGTAATACAACTTCCTTTCATGGGTGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAgtg                                   <  7:87084/144‑1 (MQ=255)
                                                                                                                 catATGACTGCGAGAATTAATCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGa                                                   >  1:261700/3‑132 (MQ=255)
                                                                                                                 catATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAgg                                           >  6:203643/3‑140 (MQ=255)
                                                                                                                 catATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAgtgt                                  <  1:339423/147‑1 (MQ=255)
                                                                                                                  atATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTa                                 >  7:175284/2‑149 (MQ=255)
                                                                                                                  atATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTa                                 >  1:367054/2‑149 (MQ=255)
                                                                                                                  atATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTa                                 >  7:359929/2‑149 (MQ=255)
                                                                                                                       actgcgaGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACg                            <  5:360028/146‑1 (MQ=255)
                                                                                                                        ctgcgaGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACgg                           >  3:130840/3‑149 (MQ=255)
                                                                                                                        ctgcgaGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACgg                           <  7:208574/147‑1 (MQ=255)
                                                                                                                         tgcgaGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAgtgt                                  >  8:165119/2‑141 (MQ=255)
                                                                                                                          gcgaGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTc                                                      <  4:260128/120‑1 (MQ=255)
                                                                                                                            gagAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAg                                         >  1:317850/3‑131 (MQ=255)
                                                                                                                            gagAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTAt                                <  1:240699/138‑1 (MQ=255)
                                                                                                                                aTTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAagag                   <  7:137719/149‑1 (MQ=255)
                                                                                                                                aTTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAagag                   >  6:237964/1‑149 (MQ=255)
                                                                                                                                 ttAACCGTGTTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCCCTAAATCGTGAACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAACTTCGTACTGAAAGAAGCAAGTCTATGACGGGCAAAGAgg                  >  2:6753/1‑149 (MQ=255)
                                                                                                                                 ttAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAgg                  >  2:174912/1‑149 (MQ=255)
                                                                                                                                 ttAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACATCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAgg                  >  3:196200/1‑149 (MQ=255)
                                                                                                                                   aaCCGTGGTGAGGCGATACCTGAACAAGTAAAACAAATTCCTCTCATGTGCGGTAGACCTCTAAATCGTTCACAGGCTCTGGCGAAGATCGCAGAAATCAAAACTAAGTTCGTAATCAAAGGAGCAAGTTTATTAAGGGCAAAGAGGc                 >  2:150424/1‑148 (MQ=255)
                                                                                                                                    aCCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCGGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCaa               <  7:398229/149‑1 (MQ=255)
                                                                                                                                     ccGTGGTGAGGCGAGCCCTGACCCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAat              <  5:249698/149‑1 (MQ=255)
                                                                                                                                       gtggtgAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAatta            <  1:125233/149‑1 (MQ=255)
                                                                                                                                          gtgAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCCTGGGCGTTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTc         <  4:90627/149‑1 (MQ=255)
                                                                                                                                          gtgAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTc         <  6:285700/149‑1 (MQ=255)
                                                                                                                                          gtgAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATCGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTc         <  3:252143/149‑1 (MQ=255)
                                                                                                                                             aggcGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGc                                        <  8:248658/112‑1 (MQ=255)
                                                                                                                                              ggcGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTa     <  8:71091/147‑1 (MQ=255)
                                                                                                                                               gcGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTAc    >  7:153202/2‑149 (MQ=255)
                                                                                                                                               aCGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTAc    <  5:173388/149‑1 (MQ=255)
                                                                                                                                                 gATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAattatt          >  8:45137/1‑141 (MQ=255)
                                                                                                                                                 gATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCt  <  8:335701/149‑1 (MQ=255)
                                                                                                                                                 gATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGAACTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCt  >  5:384832/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑TGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCT  >  NC_000913/566423‑566715

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑