Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,768,615 G→T 9.0% intergenic (+42/‑50) yfjP → / → yfjQ CP4‑57 prophage; 50S ribosome‑binding GTPase family protein/CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,768,6150GT 9.0% 406.0 / 22.7 120intergenic (+42/‑50)yfjP/yfjQCP4‑57 prophage; 50S ribosome‑binding GTPase family protein/CP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (57/52);  new base T (4/7);  total (61/59)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 120 total aligned reads displayed.

GGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTA  >  NC_000913/2768468‑2768760
                                                                                                                                                   |                                                                                                                                                 
ggTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGt                                                                                                                                                  >  4:49039/1‑149 (MQ=255)
       aCCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAg                                                                                                                                           <  4:261945/149‑1 (MQ=255)
        ccTTTCCCCTTATTCCCGCCCCGGTGCTGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCGGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAgg                                                                                                                                          <  2:65231/149‑1 (MQ=255)
        ccTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAgg                                                                                                                                          <  2:69577/149‑1 (MQ=255)
        ccTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAgg                                                                                                                                          <  3:134385/149‑1 (MQ=255)
          tttCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTACTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAGGGc                                                                                                                                        <  4:883807/149‑1 (MQ=25)
               ccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAGGGctcgga                                                                                                                                   <  1:674614/149‑5 (MQ=25)
                cgtATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttctt                                                                                                                                  <  4:821947/147‑1 (MQ=255)
                 ttATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttcttt                                                                                                                                 >  3:979486/1‑149 (MQ=255)
                 ttATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttcttt                                                                                                                                 <  3:914946/149‑1 (MQ=255)
                   aTTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTAt                                                                                                                               <  2:936831/149‑1 (MQ=255)
                   aTTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTAt                                                                                                                               <  2:235417/149‑1 (MQ=255)
                    ttCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTAtt                                                                                                                              <  2:480366/149‑1 (MQ=255)
                    ttCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGGTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTAtt                                                                                                                              <  1:599374/149‑1 (MQ=255)
                     tCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTAttt                                                                                                                             <  3:222100/149‑1 (MQ=255)
                       ccGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGt                                                                                                                           >  2:928532/1‑149 (MQ=255)
                       ccGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGt                                                                                                                           <  2:315842/149‑1 (MQ=255)
                        cGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGtt                                                                                                                          <  3:105404/149‑1 (MQ=255)
                        cGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGtt                                                                                                                          <  3:32411/149‑1 (MQ=255)
                         gCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGt                                                                                                                           <  1:337873/147‑1 (MQ=255)
                         gCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGt                                                                                                                           >  2:337873/1‑147 (MQ=255)
                           cccGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTa                                                                                                                       <  3:179580/149‑1 (MQ=255)
                                   ggACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACt                                                                                                               >  1:912093/1‑149 (MQ=255)
                                      ccATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgag                                                                                                            <  4:1016734/149‑1 (MQ=255)
                                       cATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGt                                                                                                                                                  >  3:543920/1‑110 (MQ=255)
                                       cATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGt                                                                                                                                                  <  4:543920/110‑1 (MQ=255)
                                        attattCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgagga                                                                                                          <  2:123857/149‑1 (MQ=255)
                                         ttattCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgaggag                                                                                                         >  1:731407/1‑149 (MQ=255)
                                              cAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGc                                                                                                    <  4:484049/149‑1 (MQ=255)
                                              cAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGc                                                                                                    <  1:67881/149‑1 (MQ=255)
                                                    gTGCGTTCCTCGGTGGTGTCAGTTGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTAGTCTCCGAGCCCTGTCGCAATTGCGGCAGTGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGa                                                                                              <  4:125928/149‑1 (MQ=255)
                                                     tGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGAc                                                                                             >  1:586187/1‑149 (MQ=255)
                                                        gTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGcttt                                                                                                                                     <  2:721850/106‑1 (MQ=255)
                                                        gTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGcttt                                                                                                                                     >  1:721850/1‑106 (MQ=255)
                                                        gTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCg                                                                                          >  4:446596/1‑149 (MQ=255)
                                                           ccTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgag                                                                                                            >  3:359905/1‑128 (MQ=255)
                                                           ccTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgag                                                                                                            <  4:359905/128‑1 (MQ=255)
                                                                     gtCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCgcgc                                                                             >  3:511242/1‑149 (MQ=255)
                                                                      tCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCt                                                                            >  4:555574/1‑149 (MQ=255)
                                                                      tCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCt                                                                            >  1:141891/1‑149 (MQ=255)
                                                                      tCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCt                                                                            >  4:297630/1‑149 (MQ=255)
                                                                             cccGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGgcgc                                                                     <  1:718840/149‑1 (MQ=255)
                                                                              ccGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCt                                                                    >  4:250380/1‑149 (MQ=255)
                                                                              ccGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCt                                                                    >  1:952787/1‑149 (MQ=255)
                                                                                gcgcTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGc                                                                  <  4:291725/149‑1 (MQ=255)
                                                                                  gcTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCaa                                                                >  1:916536/1‑149 (MQ=255)
                                                                                   cTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCaaa                                                               <  1:591415/149‑1 (MQ=255)
                                                                                      tCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGc                                                                  <  4:23873/143‑1 (MQ=255)
                                                                                      tCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGc                                                                  >  3:23873/1‑143 (MQ=255)
                                                                                      tCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCt                                                            >  1:480222/1‑149 (MQ=255)
                                                                                        tGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTa                                                          >  4:784317/1‑149 (MQ=255)
                                                                                        tGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTa                                                          >  4:180623/1‑149 (MQ=255)
                                                                                         gggATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTAt                                                         <  4:511242/149‑1 (MQ=255)
                                                                                             tttCTTCTTCTGAGCGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGc                                                     >  2:856947/1‑149 (MQ=255)
                                                                                              ttcttcttcTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGTCGCTGCAAATCTTAtt                                                        >  3:839830/1‑145 (MQ=255)
                                                                                              ttcttcttcTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTAtt                                                        <  4:839830/145‑1 (MQ=255)
                                                                                                 ttcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTg                                                 >  2:408095/1‑149 (MQ=255)
                                                                                                    ttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAcc                                              >  2:1008138/1‑149 (MQ=255)
                                                                                                      ctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgt                                            <  1:566948/149‑1 (MQ=255)
                                                                                                      ctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgt                                            >  4:270827/1‑149 (MQ=255)
                                                                                                      ctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgt                                            >  1:404398/1‑149 (MQ=255)
                                                                                                       tGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtc                                           >  1:168405/1‑149 (MQ=255)
                                                                                                        gAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtcc                                          >  2:649099/1‑149 (MQ=255)
                                                                                                        gAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtcc                                          <  3:777586/149‑1 (MQ=255)
                                                                                                          gtgtTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtccgt                                        <  3:297630/149‑1 (MQ=255)
                                                                                                          gtgtTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtccgt                                        <  3:299752/149‑1 (MQ=255)
                                                                                                             tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTaa                                     >  4:705682/1‑149 (MQ=255)
                                                                                                             tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTaa                                     >  1:657238/1‑149 (MQ=255)
                                                                                                                 aTTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCg                                 >  4:545359/1‑149 (MQ=255)
                                                                                                                    cTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCt                                                                                                   <  4:201538/80‑1 (MQ=255)
                                                                                                                       cccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAaga                           <  3:687629/149‑1 (MQ=255)
                                                                                                                       cccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAaga                           <  3:965149/149‑1 (MQ=255)
                                                                                                                         cTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGt                         >  3:260222/1‑149 (MQ=255)
                                                                                                                         cTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGt                         >  3:47054/1‑149 (MQ=255)
                                                                                                                          tattatTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGTCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGtt                        >  3:635971/1‑149 (MQ=255)
                                                                                                                           attattCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTa                       >  1:265501/1‑149 (MQ=255)
                                                                                                                            ttattCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTAt                      >  2:122179/1‑149 (MQ=255)
                                                                                                                               ttCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTc                   >  3:259032/1‑149 (MQ=255)
                                                                                                                                 ctcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTcgc                 <  3:742018/149‑1 (MQ=255)
                                                                                                                                 ctcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTcgc                 <  3:258095/149‑1 (MQ=255)
                                                                                                                                  tcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTg                               >  2:182037/1‑134 (MQ=255)
                                                                                                                                  tcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTg                               <  1:182037/134‑1 (MQ=255)
                                                                                                                                  tcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTcgcg                >  1:907107/1‑149 (MQ=255)
                                                                                                                                  agaaaGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTcgcg                >  3:28084/5‑149 (MQ=25)
                                                                                                                                   gaaaGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTAtt                                                        <  3:763286/105‑1 (MQ=21)
                                                                                                                                   gaaaGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTAtt                                                        >  4:763286/4‑108 (MQ=21)
                                                                                                                                   gaaaGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAaga                           >  2:937826/4‑137 (MQ=21)
                                                                                                                                   gaaaGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAaga                           <  1:937826/134‑1 (MQ=21)
                                                                                                                                   gaaaGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgt               <  4:28084/146‑1 (MQ=21)
                                                                                                                                   ccGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgt               >  3:125735/1‑149 (MQ=255)
                                                                                                                                    cGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgtg              <  3:1041423/149‑1 (MQ=255)
                                                                                                                                    cGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgtg              <  1:122179/149‑1 (MQ=255)
                                                                                                                                      aGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgtggt            <  1:408095/149‑1 (MQ=255)
                                                                                                                                      aGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgtggt            <  4:223685/149‑1 (MQ=255)
                                                                                                                                      aGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgtggt            >  3:805135/1‑149 (MQ=25)
                                                                                                                                         ccTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGcc         >  4:89055/1‑149 (MQ=255)
                                                                                                                                           tGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCa       >  2:931362/1‑149 (MQ=255)
                                                                                                                                            gTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAg      >  1:853135/1‑149 (MQ=255)
                                                                                                                                              cGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCg    <  2:731407/149‑1 (MQ=255)
                                                                                                                                               gCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGt   >  2:174527/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTA  >  NC_000913/2768468‑2768760

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: