Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,202,961 A→C 5.1% intergenic (‑28/‑63) cohE ← / → croE e14 prophage; repressor protein phage e14/e14 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,202,9610AC 5.1% 483.1 / 15.1 137intergenic (‑28/‑63)cohE/croEe14 prophage; repressor protein phage e14/e14 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (65/65);  new base C (3/4);  total (68/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.12e-01

Only 100 of 138 total aligned reads displayed.

CCATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACG  >  NC_000913/1202813‑1203106
                                                                                                                                                    |                                                                                                                                                  
ccATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTaa                                                                                                                                                    >  2:435878/1‑149 (MQ=255)
 cATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑g                                                                                                                                                  <  4:326137/149‑1 (MQ=255)
  aTAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑gg                                                                                                                                                 >  3:167041/1‑149 (MQ=255)
  aTAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑gg                                                                                                                                                 >  4:402023/1‑149 (MQ=255)
       ccTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAATGAAGAATTCCCGACAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGAGTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACaa                                                                                                                                            <  2:692342/149‑1 (MQ=255)
       ccTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCTTAAC‑GGTACaa                                                                                                                                            <  2:92393/149‑1 (MQ=255)
       ccTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACaa                                                                                                                                            <  2:145053/149‑1 (MQ=255)
       ccTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACaa                                                                                                                                            <  1:57023/149‑1 (MQ=255)
       ccTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACaa                                                                                                                                            <  1:303668/149‑1 (MQ=255)
        ctgctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCCACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTATAACCCTTGTACCGTTAA‑GGTACAAg                                                                                                                                           >  2:105878/1‑149 (MQ=255)
         tgctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCTTAAC‑GGTACAAgg                                                                                                                                          <  2:522962/149‑2 (MQ=255)
         tgctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCTTAAC‑GGTACAAgg                                                                                                                                          <  1:368950/149‑2 (MQ=255)
         tgctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGt                                                                                                                                          >  2:10085/1‑149 (MQ=255)
           ctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTAt                                                                                                                                        >  1:119767/1‑149 (MQ=255)
           ctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACC‑TTAACGGTACAAgggt                                                                                                                                        <  3:136808/149‑4 (MQ=255)
            tgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATc                                                                                                                                       >  1:632665/1‑149 (MQ=255)
              gTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACaa                                                                                                                                            <  1:164176/142‑1 (MQ=255)
              gTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACaa                                                                                                                                            >  2:164176/1‑142 (MQ=255)
                 tGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGaa                                                                                                                                  <  1:263140/149‑1 (MQ=255)
                 tGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGaa                                                                                                                                  <  4:89523/149‑1 (MQ=255)
                 tGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGaa                                                                                                                                  <  4:565561/149‑1 (MQ=255)
                 tGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGaa                                                                                                                                  <  4:467917/149‑1 (MQ=255)
                     ttGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGtt                                                                                                                              <  1:677884/149‑1 (MQ=255)
                          tCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCAttt                                                                                                                         >  2:654960/1‑149 (MQ=255)
                          tCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCAttt                                                                                                                         >  4:81142/1‑149 (MQ=255)
                             cgATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCaa                                                                                                                      <  1:125893/149‑1 (MQ=255)
                             cgATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCaa                                                                                                                      <  2:877104/149‑1 (MQ=255)
                              gATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAAt                                                                                                                     >  2:849387/1‑149 (MQ=255)
                                tcaatcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCa                                                                                                                   >  2:68022/1‑149 (MQ=255)
                                tcaatcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCa                                                                                                                   >  4:92348/1‑149 (MQ=255)
                                 caatcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCAt                                                                                                                  <  1:849387/149‑1 (MQ=255)
                                      aTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAat                                                                                                             <  3:402023/149‑1 (MQ=255)
                                       tAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAata                                                                                                            <  1:240891/149‑1 (MQ=255)
                                        aGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAatat                                                                                                           <  4:209109/149‑1 (MQ=255)
                                         gTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATg                                                                                                          <  1:123983/149‑1 (MQ=255)
                                            gCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGAGCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTac                                                                                                       >  2:264418/1‑149 (MQ=255)
                                               aaGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTCGAAGGTTCATTTCAATCATGTAATATGTACAcc                                                                                                    >  1:97460/1‑149 (MQ=255)
                                                  gaagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGa                                                                                                 >  1:753040/1‑149 (MQ=255)
                                                  gaagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGa                                                                                                 >  2:486198/1‑149 (MQ=255)
                                                  gaagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGa                                                                                                 >  3:352178/1‑149 (MQ=255)
                                                  gaagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGa                                                                                                 >  1:587632/1‑149 (MQ=255)
                                                   aagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAg                                                                                                <  4:152078/149‑1 (MQ=255)
                                                   aagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAg                                                                                                >  3:508166/1‑149 (MQ=255)
                                                    agaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAgg                                                                                               >  1:761792/1‑149 (MQ=255)
                                                        ttCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACa                                                                                           >  3:647560/1‑149 (MQ=255)
                                                        ttCCCGCCAACACCATCTCTCAGTGTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACAGGAGGTACa                                                                                           >  3:744797/1‑149 (MQ=255)
                                                          cccGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACata                                                                                         <  1:736592/149‑1 (MQ=255)
                                                           ccGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACatat                                                                                        >  1:375809/1‑149 (MQ=255)
                                                            cGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATAtt                                                                                       <  2:528949/149‑1 (MQ=255)
                                                               cAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTa                                                                                    <  1:654960/149‑1 (MQ=255)
                                                                aaCACCATCTCTCAGTTTTCTGTCGTTAGTGCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTAt                                                                                   <  2:838514/149‑1 (MQ=255)
                                                                    ccATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGaaa                                                                               <  3:571194/149‑1 (MQ=255)
                                                                    ccATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGaaa                                                                               <  2:587632/149‑1 (MQ=255)
                                                                    ccATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGaaa                                                                               <  1:105878/149‑1 (MQ=255)
                                                                    ccATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGaaa                                                                               <  4:473975/149‑1 (MQ=255)
                                                                    ccATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGaaa                                                                               <  4:281149/149‑1 (MQ=255)
                                                                    ccATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGaaa                                                                               <  1:358116/149‑1 (MQ=255)
                                                                      aTCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGc                                                                             >  3:640130/1‑149 (MQ=255)
                                                                      aTCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGc                                                                             <  1:486198/149‑1 (MQ=255)
                                                                         tctcAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTa                                                                          <  4:840911/149‑1 (MQ=255)
                                                                           ctgcgtttCTGGCTTTAGACCGCCGGGTGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTGCCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTAtt                                                                        >  3:740478/6‑149 (MQ=255)
                                                                            cAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTg                                                                       <  1:518196/149‑1 (MQ=255)
                                                                                    tGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCtt                                                               >  2:811144/1‑149 (MQ=255)
                                                                                       cGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTaa                                                            <  3:81142/149‑1 (MQ=255)
                                                                                         ttAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAAcc                                                          >  2:774567/1‑149 (MQ=255)
                                                                                            gACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCaaa                                                       >  4:214912/1‑149 (MQ=255)
                                                                                            gACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCaaa                                                       <  4:515211/149‑1 (MQ=255)
                                                                                                gccgGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGa                                                                    >  1:654443/1‑132 (MQ=255)
                                                                                                gccgGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAAc                                                   <  1:435878/149‑1 (MQ=255)
                                                                                                gccgGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAAc                                                   <  4:370138/149‑1 (MQ=255)
                                                                                                      tGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAgg                                             <  2:703967/149‑1 (MQ=255)
                                                                                                          aTGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCGTTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGa                                         >  3:542492/1‑149 (MQ=255)
                                                                                                          aTGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGa                                         >  3:769893/1‑149 (MQ=255)
                                                                                                          aTGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGAGCGa                                         >  4:139956/1‑149 (MQ=255)
                                                                                                               ttGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTgg                                    >  1:617286/1‑149 (MQ=255)
                                                                                                                      cATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGaaa                             >  3:100547/1‑149 (MQ=255)
                                                                                                                          aCGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGtt                         <  4:769893/149‑1 (MQ=255)
                                                                                                                            gAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTgg                       >  1:408172/1‑149 (MQ=255)
                                                                                                                              gaTTGAAACCCTTGTACCGTTAA‑GGTACAAGGATCTTGAAGGTTGATTTCAATCATGTAATATGTACAGCGGAGGTACATATTCTATGAAAGCGTATTGGGACTCTTTAAACAGAGAACAGCAGGGGGAGTTGGCCGGAAAAGTTGGCt                     <  4:508166/148‑1 (MQ=255)
                                                                                                                              aaTTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCt                     >  1:286726/1‑149 (MQ=255)
                                                                                                                               aTTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCt                     >  1:738392/1‑148 (MQ=255)
                                                                                                                               aTTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCt                     <  2:738392/148‑1 (MQ=255)
                                                                                                                               aTTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTc                    >  2:553769/1‑149 (MQ=255)
                                                                                                                                  aaaaCCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAac                 >  3:287839/1‑149 (MQ=255)
                                                                                                                                    aaCCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAacac               >  1:559241/1‑149 (MQ=255)
                                                                                                                                     aCCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACAcc              >  3:635466/1‑149 (MQ=255)
                                                                                                                                     aCCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACAcc              <  2:559241/149‑1 (MQ=255)
                                                                                                                                      cccTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCt             <  2:439715/149‑1 (MQ=255)
                                                                                                                                      cccTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCt             >  1:307894/1‑149 (MQ=255)
                                                                                                                                      cccTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCt             <  1:474957/149‑1 (MQ=255)
                                                                                                                                       ccTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGaa                                                    >  3:879174/1‑109 (MQ=255)
                                                                                                                                       ccTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTg            >  4:190448/1‑149 (MQ=255)
                                                                                                                                       ccTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTg            >  4:711664/1‑149 (MQ=255)
                                                                                                                                       acTTGTACCTTAAC‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTg            >  2:245023/2‑149 (MQ=255)
                                                                                                                                        cTTGTACCTTAAC‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTgg           >  4:683583/1‑149 (MQ=255)
                                                                                                                                        cTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACAGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTACCCAAAGAACAGCGGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTgg           <  4:280147/149‑1 (MQ=255)
                                                                                                                                         ttGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCaaa                                                       <  3:430692/104‑1 (MQ=255)
                                                                                                                                         ttGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCaaa                                                       >  4:430692/1‑104 (MQ=255)
                                                                                                                                             aCCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTGAATCATGTAATATGTACACCGGAGGTACATAGTGTATGAAGGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCAGGAAAAGTTGGCTCTACACCTTGCTACt      >  2:256233/1‑149 (MQ=255)
                                                                                                                                                 ttAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACg  >  1:104128/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACG  >  NC_000913/1202813‑1203106

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