Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,777,617 T→C A55A (GCT→GCC ypjF → CP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,777,6170TC93.8% 319.9 / 7.3 97A55A (GCT→GCCypjFCP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base C (47/44);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

ACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGGCGCTCTACCGGCTTGC  >  NC_000913/2777471‑2777762
                                                                                                                                                  |                                                                                                                                                 
aCAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGt                                                                                                                                                >  1:572081/1‑150 (MQ=255)
  aaTTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTAt                                                                                                                                              >  2:305781/1‑150 (MQ=255)
  aaTTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATc                                                                                                                                             >  2:1123021/1‑151 (MQ=255)
  aaTTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATc                                                                                                                                             <  2:961956/151‑1 (MQ=255)
     ttCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCAc                                                                                                                                           >  2:970455/1‑150 (MQ=255)
      tCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCActc                                                                                                                                         >  1:92157/1‑151 (MQ=255)
      tCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCActc                                                                                                                                         >  1:816326/1‑151 (MQ=255)
        gCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTg                                                                                                                                       >  1:571374/1‑151 (MQ=255)
        gCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTg                                                                                                                                       >  2:250510/1‑151 (MQ=255)
           ggcggcGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGcc                                                                                                                                    <  1:1451205/151‑1 (MQ=255)
           ggcggcGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGcc                                                                                                                                    >  2:1252948/1‑151 (MQ=255)
            gcggcgAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCg                                                                                                                                   <  2:200888/151‑1 (MQ=255)
                  aaGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAg                                                                                                                             <  1:561268/151‑1 (MQ=255)
                      cctgcctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCaa                                                                                                                         <  2:387951/151‑1 (MQ=255)
                         gcctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAAt                                                                                                                        <  2:642449/149‑1 (MQ=255)
                          cctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAAtttt                                                                                                                     >  1:1263580/1‑151 (MQ=255)
                           ctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACCCGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAAtttt                                                                                                                     <  2:678279/150‑1 (MQ=255)
                                gCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCtggtg                                                                                                               <  1:194940/151‑1 (MQ=255)
                                  cAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCtggtgg                                                                                                              >  2:1089203/1‑150 (MQ=255)
                                  cAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGa                                                                                                             <  2:933760/151‑1 (MQ=255)
                                      ggCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCTGTCAATTTTCTGGTGGaaaa                                                                                                          <  2:572015/150‑1 (MQ=255)
                                      ggCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGaaaa                                                                                                          <  2:45827/150‑1 (MQ=255)
                                      ggCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGaaa                                                                                                           >  1:204736/1‑149 (MQ=255)
                                      ggCTGTCTGGCAAATGCTACTGACACGCCTGCGGGAACACCACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAg                                                                                                         <  2:434773/151‑1 (MQ=255)
                                       gCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACGACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAg                                                                                                         >  2:1089940/1‑150 (MQ=255)
                                         tGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTAc                                                                                                      >  2:648782/1‑151 (MQ=255)
                                          gTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTa                                                                                                       <  1:229078/149‑1 (MQ=255)
                                              ggCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGa                                                                                                    <  1:1170409/148‑1 (MQ=255)
                                              ggCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGa                                                                                                    <  2:864282/148‑1 (MQ=255)
                                                cAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTgg                                                                                               <  2:631790/151‑1 (MQ=255)
                                                       tACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTAt                                                                                        >  2:754580/1‑151 (MQ=255)
                                                        aCTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATc                                                                                       <  2:1247401/151‑1 (MQ=255)
                                                        aCTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATc                                                                                       >  2:448395/1‑151 (MQ=255)
                                                           gACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGa                                                                                     <  1:615458/150‑1 (MQ=255)
                                                             cacGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATca                                                                                  <  1:252673/151‑1 (MQ=255)
                                                                   tgctgcAACAACACTATGGCCTGACACTGGACGCAACGCCGTACAGTGATGAAACTGTTATTAAGGAACATATCGACGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGCACGACCTGGTACGTATCGATCACAGAgg                                                                            >  2:1389541/1‑151 (MQ=255)
                                                                    gctgGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGa                                                                           <  2:1465287/151‑1 (MQ=255)
                                                                      tgGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGAtt                                                                         >  2:1194348/1‑151 (MQ=255)
                                                                       ggAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGAt                                                                          <  1:301918/149‑1 (MQ=255)
                                                                           caacaCTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGt                                                                    <  1:74471/151‑1 (MQ=255)
                                                                             acacTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTgg                                                                  <  1:256456/151‑1 (MQ=255)
                                                                              cacTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGc                                                                 <  1:488696/151‑1 (MQ=255)
                                                                                  aTGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGcaaca                                                             <  1:844830/151‑1 (MQ=255)
                                                                                  aTGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGcaa                                                               <  2:123058/149‑1 (MQ=255)
                                                                                   tGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGcaaca                                                             >  2:663954/1‑150 (MQ=255)
                                                                                    ggCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCaacaac                                                           >  2:1265673/1‑151 (MQ=255)
                                                                                     gCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAacaaca                                                          <  1:1226305/151‑1 (MQ=255)
                                                                                        tGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGtc                                                       >  2:1211800/1‑151 (MQ=255)
                                                                                         gACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGCATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCTTGGCAACAACAGtc                                                       >  2:135301/1‑150 (MQ=255)
                                                                                           cacTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTcc                                                    >  1:1349280/1‑151 (MQ=255)
                                                                                               gAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCatat                                                <  1:1253047/151‑1 (MQ=255)
                                                                                               gAACGACACGCCGTTCAGTGATGAAAATGTCATTAAGGAACATATCGCTGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGAACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCatat                                                >  2:1119810/1‑151 (MQ=255)
                                                                                                aaCGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCatat                                                >  1:8199/1‑150 (MQ=255)
                                                                                                  cGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATAtt                                             <  2:452159/151‑1 (MQ=255)
                                                                                                   gACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATAttt                                            <  1:978727/151‑1 (MQ=255)
                                                                                                      acGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCg                                         <  2:1124592/151‑1 (MQ=255)
                                                                                                       cGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCgt                                        <  2:354921/151‑1 (MQ=255)
                                                                                                       cGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCg                                         <  1:891658/150‑1 (MQ=255)
                                                                                                        gCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCgta                                       <  2:378636/151‑1 (MQ=255)
                                                                                                          cGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGa                                                                           >  1:615422/1‑113 (MQ=255)
                                                                                                          cGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGa                                                                           <  2:615350/113‑1 (MQ=255)
                                                                                                          cGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCgtagt                                     <  2:161325/151‑1 (MQ=255)
                                                                                                           gTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCgtagta                                    <  2:221525/151‑1 (MQ=255)
                                                                                                           gTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCgtagt                                     >  2:287347/1‑150 (MQ=255)
                                                                                                            ttCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCgtagt                                     >  1:756188/1‑149 (MQ=255)
                                                                                                            ttCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCgtagt                                     >  1:1031494/1‑149 (MQ=255)
                                                                                                             tCAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCgtagt                                     >  2:808405/1‑148 (MQ=255)
                                                                                                              cAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGat                                 <  1:878432/151‑1 (MQ=255)
                                                                                                              cAGTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGa                                  <  2:1083061/150‑1 (MQ=255)
                                                                                                              cAGTGATGAAACTGTTATTAAGGAAAATATGGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCgta                                       >  2:353739/1‑145 (MQ=255)
                                                                                                                gTGATGAAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGatat                               <  1:112936/151‑1 (MQ=255)
                                                                                                                 tgatgaAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATAtt                              >  2:537079/1‑151 (MQ=255)
                                                                                                                  gatgaAACTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTc                             <  1:393829/151‑1 (MQ=255)
                                                                                                                     gaAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTg                           <  1:948207/150‑1 (MQ=255)
                                                                                                                       aaCTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGa                        <  1:451345/151‑1 (MQ=255)
                                                                                                                         cTGTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGc                      >  1:596811/1‑151 (MQ=255)
                                                                                                                           gTTATTAAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCaa                    <  2:639777/151‑1 (MQ=255)
                                                                                                                            ttattaAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCa                     >  2:1028516/1‑149 (MQ=255)
                                                                                                                            ttattaAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGAACGAACTGGTACGTATCGATCACAGAGGTTTTTCGTGGCAACAACAGTCTCAATATATTTCCGTAGTAGATATTCTGCGAGCa                     >  2:730587/1‑149 (MQ=255)
                                                                                                                             tattaAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAgg                  >  1:60147/1‑151 (MQ=255)
                                                                                                                             tattaAGGAACATATCGACGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGAATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAgg                  >  2:1067246/1‑151 (MQ=255)
                                                                                                                              attaAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAgg                  <  1:410483/150‑1 (MQ=255)
                                                                                                                              attaAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGgc                 >  1:246917/1‑151 (MQ=255)
                                                                                                                              attaAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGgc                 >  1:382483/1‑151 (MQ=255)
                                                                                                                               ttaAGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGgcg                >  1:388605/1‑151 (MQ=255)
                                                                                                                                 aaGGAACATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGgcgc               >  2:736899/1‑150 (MQ=255)
                                                                                                                                    gAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGGCGctct            >  2:1145875/1‑150 (MQ=255)
                                                                                                                                       cATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGGCGCTCTACCg        >  2:775728/1‑151 (MQ=255)
                                                                                                                                       cATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGGCGCTCTACCg        <  2:816245/151‑1 (MQ=255)
                                                                                                                                       cATATCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGGCGCTCTACCg        >  1:229101/1‑151 (MQ=255)
                                                                                                                                        atatCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGGCGCTCTACCgg       >  2:1261437/1‑151 (MQ=255)
                                                                                                                                           tCGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGGCGCTCTACCGGCt     >  1:784752/1‑150 (MQ=255)
                                                                                                                                            cGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGGCGCTCTACCGGCtt    >  1:258470/1‑150 (MQ=255)
                                                                                                                                            cGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGGCGCTCTACCGGCtt    <  1:929395/150‑1 (MQ=255)
                                                                                                                                            cGATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGGCGCTCTACCGGCTTg   >  1:1110893/1‑151 (MQ=255)
                                                                                                                                             gATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGGCGCTCTACCGGCTTg   >  2:232325/1‑150 (MQ=255)
                                                                                                                                             gATGCCGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGGCGCTCTACCGGCTTGc  >  2:554317/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTTCCGTAGTAGATATTCTGCGAGCAAGGCGCTCTACCGGCTTGC  >  NC_000913/2777471‑2777762

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: