Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,202,961 A→C 5.3% 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.3% 724.1 / 23.7 203intergenic (‑28/‑63)cohE/croEe14 prophage; repressor protein phage e14/e14 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (92/100);  new base C (4/7);  total (96/107)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 212 total aligned reads displayed.

CCATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACG  >  NC_000913/1202813‑1203106
                                                                                                                                                    |                                                                                                                                                  
    aaaCCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTa                                                                                                                                               >  3:698405/1‑149 (MQ=255)
     aaCCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTAc                                                                                                                                              >  4:709118/1‑149 (MQ=255)
      aCCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACa                                                                                                                                             >  1:297703/1‑149 (MQ=255)
       ccTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACaa                                                                                                                                            >  2:350745/1‑149 (MQ=255)
       ccTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACC‑TTAACGGTACaa                                                                                                                                            <  2:1366441/149‑1 (MQ=255)
        ctgctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAg                                                                                                                                           >  4:259866/1‑149 (MQ=255)
        ctgctgGTCTGGGTTGGCGCGCGATCAAACATAGTGGCAAAGGAAGGGTTCACGCCAACACCAGCTATCGGTGTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTGTTCATAACGAAATTAAAACCCTTGTACCTTCAC‑GGTACAAg                                                                                                                                           <  4:183492/149‑1 (MQ=255)
         tgctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCTTAAC‑GGTACAAgg                                                                                                                                          <  3:1215641/149‑2 (MQ=255)
         tgctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGt                                                                                                                                          >  3:1068686/1‑149 (MQ=255)
         tgctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGt                                                                                                                                          >  3:55895/1‑149 (MQ=255)
         tgctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGt                                                                                                                                          >  1:1204403/1‑149 (MQ=255)
          gctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTa                                                                                                                                         <  2:192156/149‑1 (MQ=255)
           ctgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCTTAAC‑GGTACAAgggt                                                                                                                                        <  3:567080/149‑4 (MQ=255)
            tgGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATc                                                                                                                                       <  1:1309632/149‑1 (MQ=255)
             ggTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCt                                                                                                                                      <  2:570446/149‑1 (MQ=255)
              gTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAGGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCtt                                                                                                                                     <  1:82626/149‑1 (MQ=255)
                cTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGa                                                                                                                                   >  4:455395/1‑149 (MQ=255)
                    gTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGt                                                                                                                               >  1:1249051/1‑149 (MQ=255)
                     ttGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGtt                                                                                                                              <  2:1204403/149‑1 (MQ=255)
                       ggCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCa                                                                                                                            <  3:367015/149‑1 (MQ=255)
                       ggCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCa                                                                                                                            <  3:561539/149‑1 (MQ=255)
                        gCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCAt                                                                                                                           >  4:1153982/1‑149 (MQ=255)
                             cgATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCaa                                                                                                                      <  2:80611/149‑1 (MQ=255)
                              gATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTc                                                                                                                             >  3:1253784/1‑141 (MQ=255)
                              gATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTc                                                                                                                             <  4:1253784/141‑1 (MQ=255)
                              gATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAAt                                                                                                                     >  4:28142/1‑149 (MQ=255)
                               atcaatcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCTTAAC‑GGTACAAggg                                                                                                                                         >  4:652614/1‑126 (MQ=255)
                                    tcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACC‑TTAACGGTACAAgg                                                                                                                                          <  3:542769/122‑2 (MQ=255)
                                    tcaTAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACC‑TTAACGGTACAAgg                                                                                                                                          >  4:542769/1‑121 (MQ=255)
                                             cAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTaca                                                                                                      <  3:442733/149‑1 (MQ=255)
                                             cAAAGGAAGAAATCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTaa                                                                                                                                                    >  1:774306/1‑104 (MQ=255)
                                              aaaGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTacac                                                                                                     >  3:622370/1‑149 (MQ=255)
                                               aaGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACAcc                                                                                                    <  2:628572/149‑1 (MQ=255)
                                                  gaagaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGa                                                                                                 <  1:224015/149‑1 (MQ=255)
                                                     gaaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGt                                                                                              <  1:570320/149‑1 (MQ=255)
                                                      aaTTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTa                                                                                             <  3:192686/149‑1 (MQ=255)
                                                        ttCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTaa                                                                                                                                                    >  4:856267/1‑93 (MQ=255)
                                                           ccGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACatat                                                                                        <  1:884693/149‑1 (MQ=255)
                                                           ccGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACatat                                                                                        >  2:1390777/1‑149 (MQ=255)
                                                           ccGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACatat                                                                                        <  2:300834/149‑1 (MQ=255)
                                                             gCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTg                                                                                      >  2:1146967/1‑149 (MQ=255)
                                                             gCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTg                                                                                      <  1:1103103/149‑1 (MQ=255)
                                                               cAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTa                                                                                    <  2:89619/149‑1 (MQ=255)
                                                                aaCACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTAt                                                                                   <  3:6571/149‑1 (MQ=255)
                                                                    ccATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGaaa                                                                               <  2:396365/149‑1 (MQ=255)
                                                                     cATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACaa                                                                                                                                            <  2:1442708/87‑1 (MQ=255)
                                                                     cATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACaa                                                                                                                                            >  1:1442708/1‑87 (MQ=255)
                                                                      aTCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGc                                                                             >  1:815188/1‑149 (MQ=255)
                                                                         tctcAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTg                                                                                                                                    >  2:940673/1‑91 (MQ=255)
                                                                          ctcAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTAt                                                                         >  4:1107908/1‑149 (MQ=255)
                                                                           tcAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATAtt                                                                                       >  4:497505/1‑134 (MQ=255)
                                                                            cAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTg                                                                       <  3:1153982/149‑1 (MQ=255)
                                                                             aGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTa                                                                                                                                         <  2:122609/82‑1 (MQ=255)
                                                                             aGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTa                                                                                                                                         >  1:122609/1‑82 (MQ=255)
                                                                              gTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTggg                                                                     >  4:574844/1‑149 (MQ=255)
                                                                               ttttCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATGATGTAATATGTACACCGGAGGTAc                                                                                            >  3:884147/1‑125 (MQ=255)
                                                                               ttttCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGa                                                                    >  1:1077490/1‑149 (MQ=255)
                                                                               ttttCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGa                                                                    <  1:937858/149‑1 (MQ=255)
                                                                                 ttCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGAct                                                                  <  1:1066874/149‑1 (MQ=255)
                                                                                    tGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCtt                                                               <  3:813277/149‑1 (MQ=255)
                                                                                    tGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCtt                                                               >  3:348657/1‑149 (MQ=255)
                                                                                       cGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTaa                                                            <  1:60167/149‑1 (MQ=255)
                                                                                       cGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTaa                                                            <  2:128613/149‑1 (MQ=255)
                                                                                       cGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTaa                                                            <  2:480050/149‑1 (MQ=255)
                                                                                        gTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCCTTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGCATTGGGACTCTTTAAc                                                           >  4:875760/1‑149 (MQ=255)
                                                                                         ttAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAAcc                                                          <  1:522952/149‑1 (MQ=255)
                                                                                            gACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATc                                                                                                                                       <  1:1232382/69‑1 (MQ=255)
                                                                                            gACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCaaa                                                       >  1:582000/1‑149 (MQ=255)
                                                                                               cgccgGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGaa                                                    <  2:297703/149‑1 (MQ=255)
                                                                                               cgccgGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGaa                                                    <  3:28142/149‑1 (MQ=255)
                                                                                                     aTGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAAcagcag                                              <  1:330683/149‑1 (MQ=255)
                                                                                                        tCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGc                                           <  1:368978/149‑1 (MQ=255)
                                                                                                           tGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCaa                                                                                                                      <  4:616351/71‑1 (MQ=255)
                                                                                                           tGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAg                                        <  2:1249051/149‑1 (MQ=255)
                                                                                                             gATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGtt                                      <  3:489596/149‑1 (MQ=255)
                                                                                                             gATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGtt                                      >  3:53727/1‑149 (MQ=255)
                                                                                                             gATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGtt                                      >  3:1337670/1‑149 (MQ=255)
                                                                                                              aTTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTg                                     <  3:604922/149‑1 (MQ=255)
                                                                                                               ttGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTgg                                    >  3:1394899/1‑149 (MQ=255)
                                                                                                                tGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGc                                   <  1:161051/149‑1 (MQ=255)
                                                                                                                  ttttCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGaaa                                                                               >  1:655686/1‑103 (MQ=255)
                                                                                                                  ttttCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCg                                 >  1:1292190/1‑149 (MQ=255)
                                                                                                                     tCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGaa                              >  1:986298/1‑149 (MQ=255)
                                                                                                                        tAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAg                           <  2:591298/149‑1 (MQ=255)
                                                                                                                         aaCGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTTAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGt                          >  2:766611/1‑149 (MQ=255)
                                                                                                                            gAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTgg                       >  3:1266514/1‑149 (MQ=255)
                                                                                                                            gAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTgg                       <  3:978354/149‑1 (MQ=255)
                                                                                                                            gAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTgg                       >  3:38756/1‑149 (MQ=255)
                                                                                                                               aTTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTc                    <  1:1158479/149‑1 (MQ=255)
                                                                                                                                ttAAAACCCTTGTACCGTTAA‑GGTACAAGCATCTTTAAGGTTCATTTCAATCATGCAATATGTACACCGGAGGTTCATAGTGTATCAAAGCGTATTGGGACTCTTTAACCAAGGAACAGGAGTACGAGTTGGCCGGAAAAGTTGGCTCa                   >  4:905204/1‑149 (MQ=255)
                                                                                                                                    gatacTTGTACATTAAC‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAacac               >  2:1034977/5‑149 (MQ=255)
                                                                                                                                    aaCCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAacac               >  2:537760/1‑149 (MQ=255)
                                                                                                                                     aCCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACAcc              >  1:661338/1‑149 (MQ=255)
                                                                                                                                     aCCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACAcc              >  1:1141747/1‑149 (MQ=255)
                                                                                                                                      tacTTGTACCTTAAC‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCt             >  3:245514/3‑149 (MQ=255)
                                                                                                                                       ccTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAAcagca                                               <  4:373172/114‑1 (MQ=255)
                                                                                                                                        cTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTgg           <  3:1367998/149‑1 (MQ=255)
                                                                                                                                           gTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTa        >  2:1402375/1‑149 (MQ=255)
                                                                                                                                               cGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTa    >  4:80506/1‑149 (MQ=255)
                                                                                                                                                 ttAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACg  <  2:928280/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCATAAACCTGCTGGTCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACG  >  NC_000913/1202813‑1203106

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