Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,203,009 G→C 6.3% intergenic (‑76/‑15) 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,203,0090GC 6.3% 357.9 / 10.8 108intergenic (‑76/‑15)cohE/croEe14 prophage; repressor protein phage e14/e14 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (52/49);  new base C (2/5);  total (54/54)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

Only 100 of 110 total aligned reads displayed.

GGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGCTGGCTAAAAA  >  NC_000913/1202862‑1203154
                                                                                                                                                   |                                                                                                                                                 
ggAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCgg                                                                                                                                                  >  2:389853/1‑149 (MQ=255)
ggAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATCTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCgg                                                                                                                                                  >  1:976740/1‑149 (MQ=255)
       ttCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACa                                                                                                                                           >  1:1216522/1‑149 (MQ=255)
        tCCCGCCAACACCATCTCTCAGTTTTCTGGCGGTAGACCGCCGGATGTCATGGATTGTTATCATAACGAAATTAAAACCCTTTTACCGTTAAGGTACAAGTATCTTGAAGTTTCATTTCAATCATGTAATATGTACACCGGAGGTACat                                                                                                                                          >  2:329147/1‑149 (MQ=255)
          ccGCCAACACCTTCTCTCTGTTTTCTGGCGTTAGACAGCCGGATGTCATGGATTTGTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACCTCCGGTGTACatat                                                                                                                                        <  1:298434/149‑1 (MQ=255)
          ccGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACatat                                                                                                                                        >  2:889946/1‑149 (MQ=255)
            gCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTg                                                                                                                                      <  2:221247/149‑1 (MQ=255)
             ccAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACCTCCGGTGTACATAttac                                                                                                                                     <  2:62262/149‑3 (MQ=255)
             ccAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGt                                                                                                                                     >  1:458506/1‑149 (MQ=255)
             ccAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGt                                                                                                                                     >  2:327147/1‑149 (MQ=255)
             ccAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGt                                                                                                                                     >  2:1370777/1‑149 (MQ=255)
              cAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACCTCCGGTGTACATAttaca                                                                                                                                    <  1:1517454/149‑4 (MQ=255)
              cAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACCTCCGGTGTACATAttaca                                                                                                                                    <  1:917799/149‑4 (MQ=255)
              cAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTa                                                                                                                                    >  1:525623/1‑149 (MQ=255)
               aaCACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACCTCCGGTGTACATAttacat                                                                                                                                   <  2:50283/149‑5 (MQ=255)
               aaCACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTAt                                                                                                                                   <  2:1018452/149‑1 (MQ=255)
               aaCACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGATAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTAt                                                                                                                                   <  1:306391/149‑1 (MQ=255)
                   ccATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGaaa                                                                                                                               <  2:107597/149‑1 (MQ=255)
                    cATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAg                                                                                                                              >  1:1536714/1‑149 (MQ=255)
                     aTCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGc                                                                                                                             >  1:1532616/1‑149 (MQ=255)
                         ctcAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTAt                                                                                                                         >  2:1113523/1‑149 (MQ=255)
                         ctcAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTAt                                                                                                                         >  2:297091/1‑149 (MQ=255)
                            aGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTgg                                                                                                                      <  2:482057/149‑1 (MQ=255)
                            aGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTgg                                                                                                                      <  1:335326/149‑1 (MQ=255)
                            aGTTTTCTGGCGTTAGACCGCCGGATGTCACGGATTTTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTgg                                                                                                                      <  1:698691/149‑1 (MQ=255)
                              ttttCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGa                                                                                                                    <  1:1591947/149‑1 (MQ=255)
                                   tGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCtt                                                                                                               <  2:973125/149‑1 (MQ=255)
                                   tGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCtt                                                                                                               >  1:823289/1‑149 (MQ=255)
                                      cGTTAGACGGCTGGATGGCATGGATTTGTTTCATAACGAAATTAAAACCATTGTACCGTTAAGGTACAAGTATCTTGAAGGTGCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTaa                                                                                                            <  2:1176555/149‑1 (MQ=255)
                                      cGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTaa                                                                                                            <  2:1384258/149‑1 (MQ=255)
                                       gTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAAc                                                                                                           >  2:925556/1‑149 (MQ=255)
                                           gACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCaaa                                                                                                       <  1:345999/149‑1 (MQ=255)
                                             ccgccgGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGa                                                                                                     >  1:980575/1‑149 (MQ=255)
                                                ccgGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAAca                                                                                                  >  2:794313/1‑149 (MQ=255)
                                                   gATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAAcagca                                                                                               <  2:1280591/149‑1 (MQ=255)
                                                     tGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAgg                                                                                             <  1:289890/149‑1 (MQ=255)
                                                         aTGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACAGGAGGTATATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGa                                                                                         >  2:1236715/1‑149 (MQ=255)
                                                          tGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAg                                                                                        >  1:1337074/1‑149 (MQ=255)
                                                          tGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAg                                                                                        <  2:823289/149‑1 (MQ=255)
                                                          tGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAg                                                                                        <  1:889946/149‑1 (MQ=255)
                                                          tGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAg                                                                                        <  2:857917/149‑1 (MQ=255)
                                                           ggATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGt                                                                                       <  1:1113523/149‑1 (MQ=255)
                                                            gATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGtt                                                                                      >  2:121759/1‑149 (MQ=255)
                                                            gATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTAGAAGGTTCATTTCAATCATGTAATATGTACACCGGCGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAATAACAGCAGGGCGAGtt                                                                                      <  2:618100/149‑1 (MQ=255)
                                                               tGTTTTCATAACGAAATTAAAACCCTTTTACCGGTAAGGTAAAAGTAGCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGc                                                                                   <  2:262810/149‑1 (MQ=255)
                                                                gTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGcc                                                                                  >  2:1205383/1‑149 (MQ=255)
                                                                 ttttCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCg                                                                                 >  2:1020704/1‑149 (MQ=255)
                                                                     cATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGaaa                                                                             >  1:1584042/1‑149 (MQ=255)
                                                                      aTAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGaaaa                                                                            >  2:59630/1‑149 (MQ=255)
                                                                       tAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAg                                                                           <  1:503273/149‑1 (MQ=255)
                                                                       tAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAg                                                                           >  1:1467734/1‑149 (MQ=255)
                                                                       tAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGATCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAg                                                                           <  1:26066/149‑1 (MQ=255)
                                                                              aTTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCta                                                                    >  2:1046485/1‑148 (MQ=255)
                                                                              aTTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTc                                                                    <  1:1463270/149‑1 (MQ=255)
                                                                                  aaaCCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAaca                                                                <  1:414856/149‑1 (MQ=255)
                                                                                   aaCCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCCacac                                                               >  2:1206818/1‑149 (MQ=255)
                                                                                   aaCCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAacac                                                               >  2:950424/1‑149 (MQ=255)
                                                                                    aCCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACAcc                                                              <  2:841877/149‑1 (MQ=255)
                                                                                    aCCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATAACAATCATGTAATATGTTCACCGGAGGTACATTTTCTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGACTGAAAAGTTGTCTCAAAAcc                                                              >  2:1230645/1‑149 (MQ=255)
                                                                                     cccTTGTACCGTTAAGGTACAAGTATCTTTAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCt                                                             <  1:560579/149‑1 (MQ=255)
                                                                                      ccTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCAAGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTg                                                            >  2:957266/1‑149 (MQ=255)
                                                                                           tACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTACTATGTACACCGGAGGTATATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTAc                                                       <  1:1342188/149‑1 (MQ=255)
                                                                                           tACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTAc                                                       <  2:1536714/149‑1 (MQ=255)
                                                                                            aCCGTTAAGGAACAAGTAACTTGAAGGTTCATATCAATCATGTAATATGTACACCCGAGGTACATATTGTATGAAAGCGAATTGGGACACTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGCAAAAGTTGGCGCAACCCCAGCCTACt                                                      >  2:242268/1‑149 (MQ=255)
                                                                                              cttaacggTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAAACAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTa                                                    >  1:1290653/7‑149 (MQ=255)
                                                                                                ttAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACg                                                  <  1:1370777/149‑1 (MQ=255)
                                                                                                ttAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACg                                                  <  1:327147/149‑1 (MQ=255)
                                                                                                ttAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACg                                                  >  2:1085306/1‑149 (MQ=255)
                                                                                                   aGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCt                                               >  1:863108/1‑149 (MQ=255)
                                                                                                        cAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTAAACACCTGGCTACTTACGGCTGGttt                                          >  1:169082/1‑149 (MQ=255)
                                                                                                         aaGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGtttt                                         <  1:1067923/149‑1 (MQ=255)
                                                                                                         aaGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGtttt                                         <  1:1073480/149‑1 (MQ=255)
                                                                                                           gTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCa                                       >  2:76239/1‑149 (MQ=255)
                                                                                                            tATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCaa                                      <  2:268027/149‑1 (MQ=255)
                                                                                                            tATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCaa                                      <  1:19596/149‑1 (MQ=255)
                                                                                                            tATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAAt                                     >  1:1211311/1‑149 (MQ=255)
                                                                                                                  gAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGTCCGGAAAAGTTGTCTCAACACCTGGCTACATACGGCTGGTTTTCAATGGCta                                >  1:559607/1‑149 (MQ=255)
                                                                                                                   aaGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAACAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCt                                                             >  1:484075/1‑119 (MQ=255)
                                                                                                                   aaGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAACAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCt                                                             <  2:484075/119‑1 (MQ=255)
                                                                                                                    aGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCtata                              >  1:9021/1‑149 (MQ=255)
                                                                                                                        tCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGtttt                                         <  1:34264/134‑1 (MQ=255)
                                                                                                                        tCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGtttt                                         >  2:34264/1‑134 (MQ=255)
                                                                                                                        tCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATaaaaa                          >  1:515863/1‑149 (MQ=255)
                                                                                                                         cATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAAACAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATaaaaaa                         <  1:1548096/149‑1 (MQ=255)
                                                                                                                            ttCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCa                     <  2:1211311/149‑1 (MQ=255)
                                                                                                                              cAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAg                    <  2:568528/149‑1 (MQ=255)
                                                                                                                               aaTCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGt                   >  1:875711/1‑149 (MQ=255)
                                                                                                                                 tCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGttt                 <  2:1532616/149‑1 (MQ=255)
                                                                                                                                   atacaaTATGTACCTCCGGTGTACATATTATATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGttttt               >  1:297645/5‑148 (MQ=255)
                                                                                                                                   aTGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTtg               <  1:1205383/149‑1 (MQ=255)
                                                                                                                                   aTGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGAATCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTAAACACCTGGCTACTTACGCCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTtg               >  2:1111611/1‑149 (MQ=255)
                                                                                                                                      tAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGTCTACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGc            >  2:22910/1‑149 (MQ=255)
                                                                                                                                      tAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGc            >  2:1503276/1‑149 (MQ=255)
                                                                                                                                      tAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGc            >  2:30187/1‑149 (MQ=255)
                                                                                                                                       aaTATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTTGCTACTTACGTCTGGTTTTCAATGGCTATAAAAAAGCCAGGTTTTTGCt           >  2:315452/1‑149 (MQ=255)
                                                                                                                                       aaTATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGCt           >  2:1452628/1‑149 (MQ=255)
                                                                                                                                       aaTATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGCt           >  1:646309/1‑149 (MQ=255)
                                                                                                                                        atatGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGCTg          >  2:911722/1‑149 (MQ=255)
                                                                                                                                         tatGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGCTgg         <  1:1020704/149‑1 (MQ=255)
                                                                                                                                                acCGGAGGTACATATTGTATGAAAGCGTTTTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGCTGGCTaaaaa  >  2:1388061/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAAGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTTGGCTCAACACCTGGCTACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGCTGGCTAAAAA  >  NC_000913/1202862‑1203154

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