Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,768,615 G→T 5.3% 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 5.3% 594.4 / 17.2 171intergenic (+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 (88/74);  new base T (4/5);  total (92/79)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

Only 100 of 171 total aligned reads displayed.

GGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATT  >  NC_000913/2768468‑2768762
                                                                                                                                                   |                                                                                                                                                   
ggTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGt                                                                                                                                                    <  2:897363/149‑1 (MQ=255)
   aaGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCg                                                                                                                                                 >  1:614633/1‑149 (MQ=255)
      cACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCa                                                                                                                                              <  1:193060/149‑1 (MQ=255)
       aCCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAg                                                                                                                                             <  4:135604/149‑1 (MQ=255)
       aCCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGACCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAg                                                                                                                                             <  3:50533/149‑1 (MQ=25)
        ccTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAgg                                                                                                                                            <  3:988174/149‑1 (MQ=255)
        ccTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAgg                                                                                                                                            <  2:953334/149‑1 (MQ=255)
         cTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAggg                                                                                                                                           <  2:654805/149‑1 (MQ=255)
         cTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAggg                                                                                                                                           <  4:505546/149‑1 (MQ=255)
          tttCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGc                                                                                                                                          <  1:720782/149‑1 (MQ=255)
          tttCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAGGGc                                                                                                                                          <  2:710753/149‑1 (MQ=255)
             ccccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTCCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGcttt                                                                                                                                       >  3:935164/1‑149 (MQ=255)
               ccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttct                                                                                                                                     >  1:443769/1‑149 (MQ=255)
               ccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttct                                                                                                                                     >  4:310176/1‑149 (MQ=255)
               ccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttct                                                                                                                                     >  3:982036/1‑149 (MQ=255)
               ccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttct                                                                                                                                     <  3:15859/149‑1 (MQ=255)
                      cccGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCCAGTGCGGCAGGGCTTTCTTTATTTg                                                                                                                              >  4:598727/1‑149 (MQ=255)
                       ccGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGt                                                                                                                             <  1:63110/149‑1 (MQ=255)
                       ccGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGt                                                                                                                             >  1:817322/1‑149 (MQ=255)
                        cGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCAAGCCCTGCCGCACTTGCGACAg                                                                                                                                             >  1:994909/1‑132 (MQ=25)
                        cGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCAAGCCCTGCCGCACTTGCGACAg                                                                                                                                             <  2:994909/132‑1 (MQ=25)
                            ccGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAc                                                                                                                        <  3:290741/149‑1 (MQ=255)
                             cGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAcc                                                                                                                       <  4:921380/149‑1 (MQ=255)
                                 gCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATa                                                                                                                   <  1:1003655/149‑1 (MQ=255)
                                    gACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTg                                                                                                                >  2:229138/1‑149 (MQ=255)
                                    gACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTg                                                                                                                >  3:1333468/1‑149 (MQ=255)
                                       cATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgagg                                                                                                             >  1:1188989/1‑149 (MQ=255)
                                        attattCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgagga                                                                                                            <  1:923751/149‑1 (MQ=255)
                                        attattCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgagga                                                                                                            <  2:860687/149‑1 (MQ=255)
                                         ttattCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgagg                                                                                                             <  4:261664/147‑1 (MQ=255)
                                          tattCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGt                                                                                                          >  3:482029/1‑149 (MQ=255)
                                             tCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTg                                                                                                       <  2:1464940/149‑1 (MQ=255)
                                              cAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTc                                                                                                         <  3:291113/146‑1 (MQ=255)
                                              cAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTc                                                                                                         >  4:291113/1‑146 (MQ=255)
                                                  ccGTGCGTTGCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGTAGTCTGCTTAt                                                                                                  >  1:538060/1‑149 (MQ=255)
                                                     tGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGAc                                                                                               >  4:615588/1‑149 (MQ=255)
                                                     tGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGAc                                                                                               >  4:36949/1‑149 (MQ=255)
                                                             tCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTgg                                                                                       >  2:1416940/1‑149 (MQ=255)
                                                             tCGGTGGGGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTgg                                                                                       <  4:1407592/149‑1 (MQ=255)
                                                                gtggtgTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCtt                                                                                    >  1:1279368/1‑149 (MQ=255)
                                                                   gtgtCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTcgc                                                                                 >  1:1137653/1‑149 (MQ=255)
                                                                     gtCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGa                                                                                                <  4:664358/132‑1 (MQ=255)
                                                                     gtCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGa                                                                                                >  3:664358/1‑132 (MQ=255)
                                                                         gTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTg                                                                           >  1:485626/1‑149 (MQ=255)
                                                                              ccGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttc                                                                                                                                      <  2:1277146/85‑1 (MQ=255)
                                                                                gcgcTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGc                                                                    <  2:78138/149‑1 (MQ=255)
                                                                                  gcTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGc                                                                                                      <  3:1375232/113‑1 (MQ=255)
                                                                                     gTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATc                                                               <  4:1033048/149‑1 (MQ=255)
                                                                                        tGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTa                                                            >  4:1272424/1‑149 (MQ=255)
                                                                                           gATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTgg                                                                                       >  1:845067/1‑119 (MQ=255)
                                                                                            aTTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAccc                                                                                                                      <  3:735290/87‑1 (MQ=255)
                                                                                            aTTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAccc                                                                                                                      >  4:735290/1‑87 (MQ=255)
                                                                                             tttCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGCCGCTGCAAATCTTATTCGc                                                       >  3:677709/1‑149 (MQ=255)
                                                                                              ttcttcttcTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGcc                                                      <  4:86797/149‑1 (MQ=255)
                                                                                              ttcttcttcTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGcc                                                      >  3:86797/1‑149 (MQ=255)
                                                                                              ttcttcttcTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGcc                                                      <  3:860025/149‑1 (MQ=255)
                                                                                                cttcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGt                                                    >  2:614433/1‑149 (MQ=255)
                                                                                                  tcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGa                                                  <  3:310176/149‑1 (MQ=255)
                                                                                                       tGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtc                                             >  1:264679/1‑149 (MQ=255)
                                                                                                        gAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtcc                                            <  3:615588/149‑1 (MQ=255)
                                                                                                        gAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtcc                                            >  3:1313312/1‑149 (MQ=255)
                                                                                                        gAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTCCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtcc                                            <  3:581341/149‑1 (MQ=255)
                                                                                                         aGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtccg                                           >  3:143101/1‑149 (MQ=255)
                                                                                                             tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTaa                                       >  1:1394446/1‑149 (MQ=255)
                                                                                                             tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTaa                                       >  3:775233/1‑149 (MQ=255)
                                                                                                             tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTaa                                       >  1:1243222/1‑149 (MQ=255)
                                                                                                             tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTaa                                       >  4:165475/1‑149 (MQ=255)
                                                                                                             tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTaa                                       >  2:1470833/1‑149 (MQ=255)
                                                                                                             tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTaa                                       >  3:234298/1‑149 (MQ=255)
                                                                                                             tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTATTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCTCTGCAAATCTTATTCGCCGTGAccgtcc                                            <  2:1056974/144‑1 (MQ=255)
                                                                                                                     tGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGaa                               <  2:53399/149‑1 (MQ=255)
                                                                                                                      gCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAag                              >  3:47304/1‑149 (MQ=255)
                                                                                                                        ccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttct                                                                                                                                     <  2:462377/44‑1 (MQ=255)
                                                                                                                        ccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAagag                            >  1:989953/1‑149 (MQ=255)
                                                                                                                        ccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAagag                            >  2:562584/1‑149 (MQ=255)
                                                                                                                        ccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAagag                            <  3:284669/149‑1 (MQ=255)
                                                                                                                        ccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTCCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAagag                            <  4:1330174/149‑1 (MQ=255)
                                                                                                                         cTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtccgt                                          >  2:43127/1‑134 (MQ=255)
                                                                                                                          tattctTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGtt                          >  4:1181902/1‑149 (MQ=255)
                                                                                                                              attCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTa                                                            >  3:367410/1‑111 (MQ=255)
                                                                                                                              attCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTa                                                            <  4:367410/111‑1 (MQ=255)
                                                                                                                              attCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTTACCTTCCGTTAACCCGTGAAGAGTTAttt                      >  4:1195836/1‑149 (MQ=255)
                                                                                                                              attCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTAttt                      >  1:969080/1‑149 (MQ=255)
                                                                                                                                 ctcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGa                                                  >  2:1042774/1‑118 (MQ=255)
                                                                                                                                 ctcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGa                                                  <  1:1042774/118‑1 (MQ=255)
                                                                                                                                 ctcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTcgc                   <  3:1317377/149‑1 (MQ=255)
                                                                                                                                  tcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTCTTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTcgca                  >  2:1058653/1‑148 (MQ=255)
                                                                                                                                  tcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTcgcg                  >  4:777546/1‑149 (MQ=255)
                                                                                                                                      aGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgtggt              <  1:447187/149‑1 (MQ=255)
                                                                                                                                      aGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGa                                                                                                >  3:238312/1‑67 (MQ=25)
                                                                                                                                       gCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgtggtg             >  2:42850/1‑149 (MQ=255)
                                                                                                                                        cccTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAAc                                      >  4:1243193/1‑123 (MQ=255)
                                                                                                                                        cccTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAAc                                      <  3:1243193/123‑1 (MQ=255)
                                                                                                                                         ccTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGcc           >  3:587922/1‑149 (MQ=255)
                                                                                                                                         ccTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGcc           >  1:1474552/1‑149 (MQ=255)
                                                                                                                                         ccTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGaa                               >  1:1153206/1‑129 (MQ=21)
                                                                                                                                         ccTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGcc           >  2:687731/1‑149 (MQ=255)
                                                                                                                                             tCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGc       >  1:1255475/1‑149 (MQ=255)
                                                                                                                                             tCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGc       >  1:1025481/1‑149 (MQ=255)
                                                                                                                                                  aGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTAtt  <  4:482029/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATT  >  NC_000913/2768468‑2768762

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