Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 580,115 C→T 100% intergenic (+29/+136) ybcW → / ← ylcI DLP12 prophage; uncharacterized protein/DUF3950 family protein, DLP12 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,1150CT99.0% 12354.8 / 23.4 4259intergenic (+29/+136)ybcW/ylcIDLP12 prophage; uncharacterized protein/DUF3950 family protein, DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base C (25/21);  new base T (2127/2062);  total (2166/2093)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 4312 total aligned reads displayed.

GATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTC  >  NC_000913/579970‑580259
                                                                                                                                                 |                                                                                                                                                
 aTGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAATCATTCAGAGCAATTGAGGCAGCGTTGGTGTAGCACTATAATAATATGAAGGATTATTCCCTGGTGTTTGACTGATCACCATAACTTCTAATCATTCAAACTATTTa                                                                                                                                               <  4:108404/148‑1 (MQ=255)
    ttgtAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTagtct                                                                                                                                           >  7:16310/1‑145 (MQ=255)
       tAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGT                                                                                                                                         <  2:133765‑M1/148‑1 (MQ=255)
         aTTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGA                                                                                                                                       <  4:66865‑M1/148‑1 (MQ=255)
            gCATGTTTTGAACATAAGGAGTCTCTGGAAGCATTCAGAGCTATTGAGGCAGCGTTGGTGAAGCACGTTATTATTATGAAGGATTATTCCCTGGTGGTTGACTGATCACCTTAACTGCTAATCTTTCAAACTATTTAGTCTGTGACAG                                                                                                                                    <  4:314070‑M1/148‑1 (MQ=255)
             cttGTATAGCACATAAGGTGTCTCTGGATGCATTGAGAGGAATTGAGGCAGCGTTGGTGAGGCTCGATAATAATATGAAGGATTAGTGCCTGGTGGTTGACTGATCCCCATAGCTGGTAATCATTCAGCCTATGTAGTCTGTGACagag                                                                                                                                  >  8:262708/3‑149 (MQ=255)
                  ataGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCaa                                                                                                                              >  4:158552/1‑148 (MQ=255)
                   taGAACATAAGTTGTCTCTGGAACCATTAATAGCAATTGAGGCAGCGTTGTTGAAGCACGATAATAATATGAAGGATTATTCCCTGTTGGATGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAac                                                                                                                             <  3:289887/148‑1 (MQ=255)
                   taGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCCCGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAaca                                                                                                                            >  2:230933/1‑149 (MQ=255)
                   taGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCGac                                                                                                                             <  4:18762/148‑1 (MQ=255)
                   taGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAaca                                                                                                                            <  2:259171/149‑1 (MQ=255)
                     gAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACg                                                                                                                          >  1:189823/1‑149 (MQ=255)
                      aaCATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGGAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGCGCCAACACg                                                                                                                          >  6:173989/1‑148 (MQ=255)
                      aaCATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGGACGATAATAATATGAATGATTATTACCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACg                                                                                                                          >  5:306273/1‑148 (MQ=255)
                         aTAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGt                                                                                                                      >  5:230531‑M2/1‑148 (MQ=255)
                            aGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGtct                                                                                                                    >  5:150414‑M2/1‑145 (MQ=255)
                               tgtCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGtctgtca                                                                                                                >  6:353545‑M2/1‑142 (MQ=255)
                                 tctcTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCTCGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGActgcgtgttggctctgtcaca                                                                                                              <  2:103585‑M1/149‑22 (MQ=255)
                                 tctcTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTctgtcact                                                                                                              >  5:121056/1‑149 (MQ=255)
                                 tctcTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTctgtcact                                                                                                              <  8:169817/149‑1 (MQ=255)
                                   tctGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTctgtcactgt                                                                                                            >  8:286238/1‑149 (MQ=255)
                                     tGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGtct                                                                                                                    <  6:331794‑M2/139‑4 (MQ=255)
                                       gAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGCTCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAg                                                                                                         >  4:3453/1‑148 (MQ=255)
                                       gAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATTACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAg                                                                                                         >  8:398350/1‑148 (MQ=255)
                                       gAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAg                                                                                                         >  3:319543/1‑148 (MQ=255)
                                       gAAGCATTCAGAGCAAATGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTCATCACCATAACTGCTAATCATTCAAACTATTTCGTCTGTGACAGAGCCAACACGCAGTctgtcactgtca                                                                                                          <  7:185362/147‑1 (MQ=255)
                                        aaGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATATTAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTTGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGa                                                                                                       <  8:144930/149‑1 (MQ=255)
                                         aGCATTCAGCGCAATTGAGGCAGCGTTGGGGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGaa                                                                                                      <  1:338145/149‑1 (MQ=255)
                                               cAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTgg                                                                                                 <  1:74930/148‑1 (MQ=255)
                                                agagCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGt                                                                                                >  2:40903/4‑148 (MQ=255)
                                                agagCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGt                                                                                                >  6:235420/4‑148 (MQ=255)
                                                agagCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGt                                                                                                >  7:346231/4‑148 (MQ=255)
                                                  agCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTaa                                                                                              >  7:317216/2‑148 (MQ=255)
                                                    cAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAAc                                                                                           >  3:54437/1‑149 (MQ=255)
                                                    cAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAAc                                                                                           >  7:295247/1‑149 (MQ=255)
                                                      aTTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTg                                                                                         <  3:208211/149‑1 (MQ=255)
                                                      aTTGAGGCAGCGTGTGTGAAGCACGATAATAATATGAATGATTATTCCCTGTGGGGTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACAAGCAGTCTGTCACTGTCAGGAAAGGGGTAACACTg                                                                                         <  5:258656/149‑1 (MQ=255)
                                                         gAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCaa                                                                                      >  7:259002/1‑149 (MQ=255)
                                                          aGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGTTAAAACTGCAAc                                                                                     >  7:213315/1‑149 (MQ=255)
                                                            gCAGCGTTGGTGAAGCACGATAATAATATGAATGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTc                                                                                   <  6:38266/149‑1 (MQ=255)
                                                            gCAGCGATGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTTACAGAGCCAAAACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACt                                                                                    >  5:95091/1‑148 (MQ=255)
                                                              aGCGTTGGTGAAGCACTATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCaa                                                                                 <  6:269284/149‑1 (MQ=255)
                                                                        aaGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTTGTAAAACTGCAACTCAATTACTGCaa                                                                        <  1:101338/148‑1 (MQ=255)
                                                                            aCGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGccc                                                                   <  6:147546/149‑1 (MQ=255)
                                                                              gATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTc                                                                 >  6:61744/1‑149 (MQ=255)
                                                                              gATAATAATATGAAGGATTATTCCCTGGTGGTAGACTGATCACCATAACTCCAAATCATTAAAACTATTTAGTCTGTGAAAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTc                                                                 >  1:165858/1‑149 (MQ=255)
                                                                                taataataTGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAAGTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCAGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCg                                                                >  4:209029/1‑148 (MQ=255)
                                                                                 aataataTGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGt                                                               >  8:196835/1‑148 (MQ=255)
                                                                                  ataataTGAAGGATTATTCCCTGGTTGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTaa                                                             <  1:41802/149‑1 (MQ=255)
                                                                                   taataTGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGCGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAAt                                                            >  5:169689/1‑149 (MQ=255)
                                                                                   taataTGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAAt                                                            <  8:69224/149‑1 (MQ=255)
                                                                                     atatGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAAtt                                                           >  3:356887/1‑148 (MQ=255)
                                                                                      tatGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTaa                                                         <  7:396266/149‑1 (MQ=255)
                                                                                        tGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAg                                                        <  2:369936/148‑1 (MQ=255)
                                                                                        tGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGt                                                       <  5:205507/149‑1 (MQ=255)
                                                                                        tGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGt                                                       <  3:161434/149‑1 (MQ=255)
                                                                                         gAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGc                                                                          <  6:85623/129‑1 (MQ=255)
                                                                                          aaGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGa                                                     <  4:2320/149‑1 (MQ=255)
                                                                                           aGGATTATTCCCTGGTGGTTTACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGaa                                                    <  5:283910/149‑1 (MQ=255)
                                                                                           aGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTaa                                                         <  3:102803/144‑1 (MQ=255)
                                                                                           aGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGaa                                                    >  1:282662/1‑149 (MQ=255)
                                                                                           aGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGaa                                                    >  3:128190/1‑149 (MQ=255)
                                                                                               ttattCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTa                                                <  5:355043/149‑1 (MQ=255)
                                                                                                 attCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTAc                                               <  7:70974/148‑1 (MQ=255)
                                                                                                 attCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTAc                                               <  5:63283/148‑1 (MQ=255)
                                                                                                   tCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTAc                                               <  4:42149/146‑1 (MQ=255)
                                                                                                    cccTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAAt                                                            <  2:72374/132‑1 (MQ=255)
                                                                                                     ccTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAatat                                          >  7:118503/1‑149 (MQ=255)
                                                                                                        ggtggtTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTg                                                                                                  >  2:345000/1‑90 (MQ=255)
                                                                                                         gtgggTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTc                                      <  7:371450/149‑1 (MQ=255)
                                                                                                            gtTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCt                                    <  1:249883/148‑1 (MQ=255)
                                                                                                            gtTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGACACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTg                                   >  4:131402/1‑149 (MQ=255)
                                                                                                             ttGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGt                                  >  1:220381/1‑149 (MQ=255)
                                                                                                               gACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTa                                                <  8:372110/133‑1 (MQ=255)
                                                                                                                aCTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCg                               >  3:239712/1‑149 (MQ=255)
                                                                                                                aCTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCg                               <  3:77107/149‑1 (MQ=255)
                                                                                                                  tGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCg                               >  7:142100/1‑147 (MQ=255)
                                                                                                                  tGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGa                             <  2:167235/149‑1 (MQ=255)
                                                                                                                  tGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCg                               >  6:303149/1‑147 (MQ=255)
                                                                                                                     tCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAggg                          <  8:314841/149‑1 (MQ=255)
                                                                                                                       aCCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCt                                    >  2:289660/1‑137 (MQ=255)
                                                                                                                        ccATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTc                                      >  4:262000/1‑134 (MQ=255)
                                                                                                                              cTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAAcgcg                 <  8:122931/149‑1 (MQ=255)
                                                                                                                               tGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCgg                >  3:186964/1‑149 (MQ=255)
                                                                                                                                gCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCgg                >  6:2017/1‑148 (MQ=255)
                                                                                                                                 cTAATCATTCAAACTATTTAGTCTGTGACAGATCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCgg                <  5:40575/147‑1 (MQ=255)
                                                                                                                                  tAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGAt             >  5:334297/1‑149 (MQ=255)
                                                                                                                                   aaTCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCCACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATg            >  3:26554/1‑149 (MQ=255)
                                                                                                                                   aaTCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAAcgc                  >  1:365529/1‑143 (MQ=255)
                                                                                                                                    aTCATTCAAACTATTTAGTCTGTGACAGTGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGAt             <  8:344656/147‑1 (MQ=255)
                                                                                                                                     tCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTg                                                                           >  4:28407/1‑84 (MQ=255)
                                                                                                                                       aatCAAACTATTTAGTCTGTGACAGAGCCAACACCCAGTCTGTCACTGTCAGGAAAGTGTTAAAACTGCAACTCAATTACTGAAATGCACTCGAAATTAAGTGAATATACAATATCGTCCTGTTCGTAGTGAAGAACGCGGGATGttaa        >  2:319180/3‑147 (MQ=255)
                                                                                                                                       aTTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATAACGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCa        >  7:236717/1‑149 (MQ=255)
                                                                                                                                        ttCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCAt       >  6:116526/1‑149 (MQ=255)
                                                                                                                                           aaaCTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCgg                <  1:133811/137‑1 (MQ=255)
                                                                                                                                           aaaCTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCAttct    >  7:96811/1‑149 (MQ=255)
                                                                                                                                           aaaCTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCAttct    >  1:308167/1‑149 (MQ=255)
                                                                                                                                            aaCTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATc                                         >  8:225596/1‑111 (MQ=255)
                                                                                                                                            aaCTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCAttctt   >  4:91391/1‑149 (MQ=255)
                                                                                                                                             actattTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCAttcttc  <  2:197918/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTC  >  NC_000913/579970‑580259

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