Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,768,609 T→C 6.2% intergenic (+36/‑56) 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,6090TC 6.2% 575.5 / 16.4 174intergenic (+36/‑56)yfjP/yfjQCP4‑57 prophage; 50S ribosome‑binding GTPase family protein/CP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (88/75);  new base C (5/6);  total (93/81)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 177 total aligned reads displayed.

TCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGC  >  NC_000913/2768461‑2768757
                                                                                                                                                    |                                                                                                                                                    
tCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGt                                                                                                                                                      <  1:1369921/149‑1 (MQ=255)
 cGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTc                                                                                                                                                     <  3:595096/149‑1 (MQ=255)
  gATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCGGAGCCCTGTCg                                                                                                                                                    >  4:1026234/1‑149 (MQ=255)
     aCGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGGCCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCaa                                                                                                                                                 <  3:1039022/149‑1 (MQ=255)
       ggTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGt                                                                                                                                               <  1:877830/149‑1 (MQ=255)
       ggTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGt                                                                                                                                               <  2:897363/149‑1 (MQ=255)
          aaGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCg                                                                                                                                            >  1:614633/1‑149 (MQ=255)
             cACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCa                                                                                                                                         <  1:193060/149‑1 (MQ=255)
               ccTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAgg                                                                                                                                       <  3:988174/149‑1 (MQ=255)
               ccTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAgg                                                                                                                                       <  2:953334/149‑1 (MQ=255)
                cTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAggg                                                                                                                                      >  3:691967/1‑149 (MQ=255)
                cTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAggg                                                                                                                                      <  2:654805/149‑1 (MQ=255)
                cTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAggg                                                                                                                                      <  4:505546/149‑1 (MQ=255)
                cTTTCCCCTCATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAggg                                                                                                                                      <  4:1427167/149‑1 (MQ=255)
                 tttCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGc                                                                                                                                     <  1:720782/149‑1 (MQ=255)
                     cccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttc                                                                                                                                 <  3:970775/149‑1 (MQ=255)
                     cccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttc                                                                                                                                 >  4:1062521/1‑149 (MQ=255)
                     cccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttc                                                                                                                                 >  3:347019/1‑149 (MQ=255)
                      ccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCaa                                                                                                                                                 <  3:659840/132‑1 (MQ=255)
                      ccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttct                                                                                                                                <  3:15859/149‑1 (MQ=255)
                      ccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttct                                                                                                                                >  1:443769/1‑149 (MQ=255)
                      ccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttct                                                                                                                                >  2:572019/1‑149 (MQ=255)
                       cTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttctt                                                                                                                               >  1:1464650/1‑149 (MQ=255)
                       cTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCacttgc                                                                                                                                             >  2:989892/1‑130 (MQ=255)
                       cTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCacttgc                                                                                                                                             <  1:989892/135‑6 (MQ=255)
                             cccGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCCAGTGCGGCAGGGCTTTCTTTATTTg                                                                                                                         >  4:598727/1‑149 (MQ=255)
                              ccGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGt                                                                                                                        >  1:817322/1‑149 (MQ=255)
                              ccGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGt                                                                                                                        <  1:63110/149‑1 (MQ=255)
                               cGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCAAGCCCTGCCGCACTTGCGACAg                                                                                                                                        >  1:994909/1‑132 (MQ=25)
                                gCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGttt                                                                                                                      >  4:303164/1‑149 (MQ=255)
                                   ctggTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAc                                                                                                                   <  4:1251976/147‑1 (MQ=255)
                                   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)
                                             ccATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgag                                                                                                         <  4:316702/149‑1 (MQ=255)
                                              cATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgagg                                                                                                        >  1:1188989/1‑149 (MQ=255)
                                               attattCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgagga                                                                                                       <  2:860687/149‑1 (MQ=255)
                                               attattCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgagga                                                                                                       <  1:923751/149‑1 (MQ=255)
                                                ttattCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgaggag                                                                                                      >  3:1407592/1‑149 (MQ=255)
                                                ttattCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgagg                                                                                                        <  4:261664/147‑1 (MQ=255)
                                                ttattCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgagg                                                                                                        >  3:261664/1‑147 (MQ=255)
                                                    tCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTg                                                                                                  <  4:1285725/149‑1 (MQ=255)
                                                         ccGTGCGTTGCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGTAGTCTGCTTAt                                                                                             >  1:538060/1‑149 (MQ=255)
                                                           gTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGa                                                                                           <  3:303164/149‑1 (MQ=255)
                                                            tGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGAc                                                                                          <  3:1062521/149‑1 (MQ=255)
                                                            tGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGAc                                                                                          >  4:36949/1‑149 (MQ=255)
                                                              cGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGAccc                                                                                        <  2:614633/149‑1 (MQ=255)
                                                                    tCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTgg                                                                                  >  2:1416940/1‑149 (MQ=255)
                                                                          gtgtCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTcgc                                                                            >  1:1137653/1‑149 (MQ=255)
                                                                            gtCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGa                                                                                           <  4:664358/132‑1 (MQ=255)
                                                                                gTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTg                                                                      >  1:485626/1‑149 (MQ=255)
                                                                                gTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTg                                                                      <  4:702448/149‑1 (MQ=255)
                                                                                     ccGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttc                                                                                                                                 >  1:1277146/1‑85 (MQ=255)
                                                                                     ccGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttc                                                                                                                                 <  2:1277146/85‑1 (MQ=255)
                                                                                       gcgcTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGc                                                               >  1:1172015/1‑149 (MQ=255)
                                                                                       gcgcTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGc                                                               <  2:78138/149‑1 (MQ=255)
                                                                                         gcTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGc                                                                                                 >  4:1375232/1‑113 (MQ=255)
                                                                                                   aTTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAccc                                                                                                                 <  3:735290/87‑1 (MQ=255)
                                                                                                    tttCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGCCGCTGCAAATCTTATTCGc                                                  >  3:677709/1‑149 (MQ=255)
                                                                                                     ttcttcttcTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGcc                                                 <  3:860025/149‑1 (MQ=255)
                                                                                                     ttcttcttcTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGcc                                                 <  4:86797/149‑1 (MQ=255)
                                                                                                     ttcttcttcTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGcc                                                 >  3:86797/1‑149 (MQ=255)
                                                                                                        ttcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTg                                              >  3:49986/1‑149 (MQ=255)
                                                                                                           ttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAcc                                           >  2:1048922/1‑149 (MQ=255)
                                                                                                              tGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtc                                        >  1:264679/1‑149 (MQ=255)
                                                                                                               gAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtcc                                       <  3:615588/149‑1 (MQ=255)
                                                                                                               gAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtcc                                       >  3:1313312/1‑149 (MQ=255)
                                                                                                                aGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtccg                                      >  3:143101/1‑149 (MQ=255)
                                                                                                                    tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTaa                                  >  3:775233/1‑149 (MQ=255)
                                                                                                                    tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTaa                                  >  1:1243222/1‑149 (MQ=255)
                                                                                                                    tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTaa                                  >  3:234298/1‑149 (MQ=255)
                                                                                                                    tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTaa                                  >  2:1470833/1‑149 (MQ=255)
                                                                                                                         ttCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGt                             >  2:1430454/1‑149 (MQ=255)
                                                                                                                            tGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGaa                          <  2:53399/149‑1 (MQ=255)
                                                                                                                             gCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAag                         >  3:47304/1‑149 (MQ=255)
                                                                                                                             gCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAag                         >  4:483999/1‑149 (MQ=255)
                                                                                                                               ccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttct                                                                                                                                <  2:462377/44‑1 (MQ=255)
                                                                                                                               ccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAagag                       <  3:284669/149‑1 (MQ=255)
                                                                                                                               ccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTCCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAagag                       <  4:1330174/149‑1 (MQ=255)
                                                                                                                                cTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgag                                                                                                         <  1:1337307/66‑1 (MQ=255)
                                                                                                                                cTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtccgt                                     >  2:43127/1‑134 (MQ=255)
                                                                                                                                 tattctTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGtt                     >  4:1181902/1‑149 (MQ=255)
                                                                                                                                     attCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTa                                                       >  3:367410/1‑111 (MQ=255)
                                                                                                                                     attCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTTACCTTCCGTTAACCCGTGAAGAGTTAttt                 >  4:1195836/1‑149 (MQ=255)
                                                                                                                                     attCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTAttt                 >  1:969080/1‑149 (MQ=255)
                                                                                                                                      ttCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTc                >  4:626124/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                                                                                           <  4:238312/67‑1 (MQ=25)
                                                                                                                                              gCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgtggtg        <  3:562583/149‑1 (MQ=255)
                                                                                                                                              gCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgtggtg        >  2:42850/1‑149 (MQ=255)
                                                                                                                                                ccTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGcc      >  3:587922/1‑149 (MQ=255)
                                                                                                                                                ccTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGcc      >  2:687731/1‑149 (MQ=255)
                                                                                                                                                    tCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGc  >  1:1025481/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGC  >  NC_000913/2768461‑2768757

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: