Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,202,957 G→T 5.1% intergenic (‑24/‑67) 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,9570GT 5.1% 471.6 / 12.2 139intergenic (‑24/‑67)cohE/croEe14 prophage; repressor protein phage e14/e14 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (64/66);  new base T (3/4);  total (67/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.61e-01

Only 100 of 139 total aligned reads displayed.

CTCCCATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACT  >  NC_000913/1202810‑1203102
                                                                                                                                                   |                                                                                                                                                  
 tcCCATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGtt                                                                                                                                                  <  4:523114/149‑1 (MQ=255)
 tcCCATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGtt                                                                                                                                                  <  3:380621/149‑1 (MQ=255)
   ccATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTaa                                                                                                                                                >  2:435878/1‑149 (MQ=255)
     aTAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑gg                                                                                                                                             >  3:167041/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                                                                                                                                      <  1:122427/149‑2 (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)
            tgctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACC‑TTAACGGTACAAgg                                                                                                                                      <  1:346094/149‑2 (MQ=255)
              ctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTAt                                                                                                                                    >  1:119767/1‑149 (MQ=255)
               tgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATc                                                                                                                                   >  1:632665/1‑149 (MQ=255)
                 gTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACaa                                                                                                                                        >  2:164176/1‑142 (MQ=255)
                 gTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACaa                                                                                                                                        <  1:164176/142‑1 (MQ=255)
                 gTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCtt                                                                                                                                 >  3:473975/1‑149 (MQ=255)
                 gTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCtt                                                                                                                                 <  4:116848/149‑1 (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                                                                                                                  <  2:877104/149‑1 (MQ=255)
                                cgATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCaa                                                                                                                  <  1:125893/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)
                                        caTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTaa                                                                                                          <  4:348699/149‑1 (MQ=255)
                                         aTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAat                                                                                                         <  3:402023/149‑1 (MQ=255)
                                           aGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAatat                                                                                                       >  3:674696/1‑149 (MQ=255)
                                              ggCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTa                                                                                                    >  4:407396/1‑149 (MQ=255)
                                               gCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGAGCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTac                                                                                                   >  2:264418/1‑149 (MQ=255)
                                               gCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTaa                                                                                                          <  1:19928/142‑1 (MQ=255)
                                                  aaGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTCGAAGGTTCATTTCAATCATGTAATATGTACAcc                                                                                                >  1:97460/1‑149 (MQ=255)
                                                     gaagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGa                                                                                             >  2:486198/1‑149 (MQ=255)
                                                     gaagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGa                                                                                             >  3:152078/1‑149 (MQ=255)
                                                     gaagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGa                                                                                             >  1:587632/1‑149 (MQ=255)
                                                      aagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAg                                                                                            >  3:508166/1‑149 (MQ=255)
                                                      aagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAg                                                                                            <  4:152078/149‑1 (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)
                                                                      acCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAAGCATGTAATATGTACACCGGAGGTACGTATTGTATGaa                                                                            >  3:574105/1‑149 (MQ=255)
                                                                       ccATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGaaa                                                                           <  1:105878/149‑1 (MQ=255)
                                                                       ccATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGaaa                                                                           <  2:587632/149‑1 (MQ=255)
                                                                       ccATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGaaa                                                                           <  4:281149/149‑1 (MQ=255)
                                                                       ccATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGaaa                                                                           <  3:571194/149‑1 (MQ=255)
                                                                         aTCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGc                                                                         <  1:486198/149‑1 (MQ=255)
                                                                            tctcAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTa                                                                      <  4:840911/149‑1 (MQ=255)
                                                                              tcAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCTTAtt                                                                    >  1:642595/1‑149 (MQ=255)
                                                                               cAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTg                                                                   <  1:518196/149‑1 (MQ=255)
                                                                                aGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAGAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTgg                                                                  <  1:739195/149‑1 (MQ=255)
                                                                                       tGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCtt                                                           >  2:811144/1‑149 (MQ=255)
                                                                                            ttAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAAcc                                                      >  2:774567/1‑149 (MQ=255)
                                                                                               gACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCaaa                                                   >  4:214912/1‑149 (MQ=255)
                                                                                                   gccgGATGTCATGGATTGTTTTCATAACGAAATTAAAACCTTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAAGCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAAc                                               <  4:574105/149‑1 (MQ=255)
                                                                                                   gccgGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAAc                                               <  4:370138/149‑1 (MQ=255)
                                                                                                   gccgGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAAc                                               <  1:435878/149‑1 (MQ=255)
                                                                                                       gATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAAcagca                                           >  2:843465/1‑149 (MQ=255)
                                                                                                         tGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAgg                                         <  2:703967/149‑1 (MQ=255)
                                                                                                             aTGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGa                                     >  3:769893/1‑149 (MQ=255)
                                                                                                                gATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGtt                                  >  3:356497/1‑149 (MQ=255)
                                                                                                                gATTGTTCTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGtt                                  >  1:439715/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)
                                                                                                                                aaaTTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGc                  <  3:751786/149‑1 (MQ=255)
                                                                                                                                aaaTTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGc                  <  4:490546/149‑1 (MQ=255)
                                                                                                                                 gaTTGAAACCCTTGTACCGTTAA‑GGTACAAGGATCTTGAAGGTTGATTTCAATCATGTAATATGTACAGCGGAGGTACATATTCTATGAAAGCGTATTGGGACTCTTTAAACAGAGAACAGCAGGGGGAGTTGGCCGGAAAAGTTGGCt                 <  4:508166/148‑1 (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)
                                                                                                                                    tAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCaa              <  4:167041/149‑1 (MQ=255)
                                                                                                                                     aaaaCCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGTGAGTTGGCCGGAAAAGTTGGCTCAac             >  1:264315/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         >  1:307894/1‑149 (MQ=255)
                                                                                                                                         cccTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCt         <  2:439715/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                                                   >  4:430692/1‑104 (MQ=255)
                                                                                                                                                aCCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTGAATCATGTAATATGTACACCGGAGGTACATAGTGTATGAAGGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCAGGAAAAGTTGGCTCTACACCTTGCTACt  >  2:256233/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTCCCATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACT  >  NC_000913/1202810‑1203102

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: ‑