Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 265,181 A→T 5.1% L9H (CTC→CAC)  yafX ← CP4‑6 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913265,1810AT 5.1% 680.8 / 8.3 199L9H (CTC→CAC) yafXCP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (96/90);  new base T (2/11);  total (98/101)
Fisher's exact test for biased strand distribution p-value = 1.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

CACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGAT  >  NC_000913/265032‑265330
                                                                                                                                                     |                                                                                                                                                     
cacaGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGa                                                                                                                                                       <  2:1233694/150‑1 (MQ=255)
 acaGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAg                                                                                                                                                      <  1:1192447/150‑1 (MQ=255)
 acaGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAg                                                                                                                                                      <  1:2820657/150‑1 (MQ=255)
   aGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGt                                                                                                                                                    >  2:640797/1‑150 (MQ=255)
      aGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGt                                                                                                                                                 >  1:1883615/1‑150 (MQ=255)
        cGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGtgc                                                                                                                                               >  2:145475/1‑150 (MQ=255)
        cGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGtgc                                                                                                                                               >  2:2444011/1‑150 (MQ=255)
          gTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGtgctg                                                                                                                                             >  1:2425798/1‑150 (MQ=255)
               tccaGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCt                                                                                                                                        >  1:1063795/1‑150 (MQ=255)
               tccaGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCt                                                                                                                                        >  1:1544138/1‑150 (MQ=255)
                 caGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCtgt                                                                                                                                      <  2:1258800/150‑1 (MQ=255)
                 caGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCtgt                                                                                                                                      >  1:2906641/1‑150 (MQ=255)
                 caGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCtgt                                                                                                                                      <  1:2070198/150‑1 (MQ=255)
                    gCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGtgtt                                                                                                                                   >  2:2462893/1‑150 (MQ=255)
                     cGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGtgttg                                                                                                                                  <  1:851307/150‑1 (MQ=255)
                     cGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGtgttg                                                                                                                                  <  2:3039252/150‑1 (MQ=255)
                        aGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCa                                                                                                                               <  2:2564885/150‑1 (MQ=255)
                          atatGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGTGGTCGTGCTGCGTCTGTGTTGTCATg                                                                                                                             >  2:1541716/1‑150 (MQ=255)
                          ataaGCGGCTCCAGGAGTAGCCACAGCGGGATAGCGCCAAAGTACTGTGGACAGAAATAGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTATCGGGAGCGTGGTCGGGCTGCGTCTGTGTTGTCATg                                                                                                                             >  2:458208/1‑150 (MQ=255)
                           tatGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATgg                                                                                                                            <  2:1471688/150‑1 (MQ=255)
                           tatGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATgg                                                                                                                            <  2:2406517/150‑1 (MQ=255)
                           tatGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATgg                                                                                                                            <  1:1819063/150‑1 (MQ=255)
                           tatGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAAGTCGGGGTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATgg                                                                                                                            >  2:898865/1‑150 (MQ=255)
                              gCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTa                                                                                                                         >  1:2918475/1‑150 (MQ=255)
                               cGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTAt                                                                                                                        >  1:244186/1‑150 (MQ=255)
                                ggCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATg                                                                                                                       >  2:75042/1‑150 (MQ=255)
                                 gCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGt                                                                                                                      >  1:7327/1‑150 (MQ=255)
                                 gCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGt                                                                                                                      >  1:2513594/1‑150 (MQ=255)
                                  cTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGtc                                                                                                                     <  1:2116618/150‑1 (MQ=255)
                                    ccAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGtctc                                                                                                                   <  1:403085/150‑1 (MQ=255)
                                      aGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCg                                                                                                                 >  2:1364434/1‑150 (MQ=255)
                                           tAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAAt                                                                                                            <  1:1351284/150‑1 (MQ=255)
                                            aGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATa                                                                                                           >  1:1172600/1‑150 (MQ=255)
                                            aGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATa                                                                                                           >  1:1853376/1‑150 (MQ=255)
                                              ccACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATaaa                                                                                                         >  1:151131/1‑150 (MQ=255)
                                                aCTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATaaaaa                                                                                                       <  2:1608794/150‑1 (MQ=255)
                                                            cgcCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCTTGACGGTCTGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGTTCTTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGa                                                                                           <  2:2608986/150‑1 (MQ=255)
                                                            cgcCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGa                                                                                           <  2:1812392/150‑1 (MQ=255)
                                                             gcCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGAt                                                                                          >  2:2159179/1‑150 (MQ=255)
                                                               cAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGAtgg                                                                                        >  2:527022/1‑150 (MQ=255)
                                                                 aaGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGAtggct                                                                                      >  2:1285884/1‑150 (MQ=255)
                                                                  aGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGAtggctg                                                                                     <  2:1829912/150‑1 (MQ=255)
                                                                       tgtgGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTa                                                                                >  2:224595/1‑150 (MQ=255)
                                                                        gtgGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGtct                                                                               >  1:722921/1‑148 (MQ=255)
                                                                          ggCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATgg                                                                             <  2:2918475/150‑1 (MQ=255)
                                                                            ccAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATgggg                                                                           >  2:667392/1‑150 (MQ=255)
                                                                            ccAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATgggg                                                                           <  2:1086868/150‑1 (MQ=255)
                                                                              aGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGat                                                                         >  2:119612/1‑150 (MQ=255)
                                                                              aGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGat                                                                         >  2:2932474/1‑150 (MQ=255)
                                                                                aaaTCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGatat                                                                       >  2:2728568/1‑150 (MQ=255)
                                                                                aaaTCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGagat                                                                       >  2:1778501/1‑150 (MQ=255)
                                                                                   tCGATACGGTGTTCATGGGGGACGGGCGTGAGGGTCTTTTCAAATTGGGGTTGGTTAGCGGGTGCGTGGTGGTGCTGGGTCTTTGTTTTCTTGGGTATGTCTCCGAAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATaaa                                                                    <  1:165737/150‑1 (MQ=21)
                                                                                   tCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATaaa                                                                    >  1:1557221/1‑150 (MQ=255)
                                                                                     gATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAATCGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAgt                                                                  >  1:1200386/1‑150 (MQ=255)
                                                                                          tcTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCGAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGtt                                                             >  1:860585/2‑150 (MQ=255)
                                                                                          tcTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGtt                                                             >  1:610791/2‑150 (MQ=255)
                                                                                           cTGTTCATCGGGGACCGGCGTGACGGTCAGTTAAAATTGGGGTTGGTTTGGGGGTGCTAGGTGTTGCTGCTTCTTTGTTGTCTTGGGTATGTCTCCTTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTc                                                            <  1:1589408/150‑1 (MQ=25)
                                                                                             gTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTTTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCgg                                                          <  1:1599793/150‑1 (MQ=255)
                                                                                             gTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCgg                                                          >  2:269576/1‑150 (MQ=255)
                                                                                              ttCATGGGGGACCGGGGTGACGGTCTGTTCAATTTGGGGTTGGTTTGCGGGTGCTTGGTGGTGCTGGGTGTGTTTTGTCTTGGGTATGTCTCCGTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCggg                                                         <  1:1449386/150‑1 (MQ=21)
                                                                                              ttCATGGGGGACCGGGGTGACGGTCAGTTCAATTTCGGGTTGGTTCGGGGGTGGGTGGTCGTGCTGCGTCTGTGTTGTCTTGGGTTTGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCggg                                                         <  1:2990147/150‑1 (MQ=25)
                                                                                              ttCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCggg                                                         <  2:1901363/150‑1 (MQ=255)
                                                                                                cATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGa                                                       <  2:2437383/150‑1 (MQ=255)
                                                                                                 aTCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAg                                                      <  2:2335048/150‑1 (MQ=255)
                                                                                                   gggggACGGGCGTGACGGTCTGTTCAAATTCGGGTTGTTTTGGGGGTGCGTGGTCGTGGTGCGTCGGTGTTGTCTTGGTTATGTATCCGTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGt                                                    <  1:2627683/149‑1 (MQ=25)
                                                                                                   cGGGGACCGGCGTGACGTTCAGTTCAAATTCGGTTTGTTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGTTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGt                                                    <  1:932325/150‑1 (MQ=255)
                                                                                                   cGGGGACCGGCGTGACGGTCAGTTCATATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGt                                                    <  2:2100146/150‑1 (MQ=255)
                                                                                                   cGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGt                                                    <  2:2208193/150‑1 (MQ=255)
                                                                                                   cGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGt                                                    <  2:2282976/150‑1 (MQ=255)
                                                                                                     gggACGGGCTTGACGGTCTGTTCAATTTGGGGTTGGTTTGGGGGTGCGTGGTGGTGCTGCTTGTGTTTTGTCTTGGTTTTGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGa                                                  <  1:2947993/150‑1 (MQ=25)
                                                                                                      ggACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGaa                                                 >  1:2177386/1‑150 (MQ=255)
                                                                                                         ccGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGc                                              >  1:2166700/1‑150 (MQ=255)
                                                                                                          cGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCg                                             >  1:1345229/1‑150 (MQ=255)
                                                                                                           ggCGTGCCGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCATTAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCgg                                            <  1:1456543/150‑1 (MQ=255)
                                                                                                            gCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCggg                                           <  1:326092/150‑1 (MQ=255)
                                                                                                            gCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCggg                                           <  1:1794822/150‑1 (MQ=255)
                                                                                                                 aCGGTCTGTTCACATTCGGGTGGGTTCGCGGGTGGGTGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAACACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAAt                                      <  1:2653296/150‑1 (MQ=25)
                                                                                                                   ggTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGc                                    >  1:2884589/1‑150 (MQ=255)
                                                                                                                    gTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCt                                   >  2:2833828/1‑150 (MQ=255)
                                                                                                                      cAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTc                                 <  1:1017909/150‑1 (MQ=255)
                                                                                                                       aGTTCAAATTCGGGTTGGTTAGCGGGTGCGTGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTcg                                <  1:2569869/150‑1 (MQ=255)
                                                                                                                       aGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTcg                                <  1:769632/150‑1 (MQ=255)
                                                                                                                        gTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTcgc                               <  1:104401/150‑1 (MQ=255)
                                                                                                                        gTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTcgc                               <  1:3037356/150‑1 (MQ=255)
                                                                                                                        gTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTcgc                               >  1:698820/1‑150 (MQ=255)
                                                                                                                         ttCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTcgcg                              >  2:1871915/1‑150 (MQ=255)
                                                                                                                         ttCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCCGGGAGGTGAATGCGGGTAAATGCTTcgcg                              >  2:1850760/1‑150 (MQ=255)
                                                                                                                           cAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCTTGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGAt                            <  2:1429785/150‑1 (MQ=255)
                                                                                                                             atttCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCa                          <  2:2122798/148‑1 (MQ=255)
                                                                                                                                tCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCg                       <  2:1345229/150‑1 (MQ=255)
                                                                                                                                 cGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTTTATTGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTCTGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCgg                      <  2:2356879/150‑1 (MQ=255)
                                                                                                                                 cGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCgg                      <  2:1123952/150‑1 (MQ=255)
                                                                                                                                  gggTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCggg                     <  2:1522015/150‑1 (MQ=255)
                                                                                                                                    gttggttAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTAATGGGTATGTCTCCGTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTg                   <  1:2653883/150‑1 (MQ=255)
                                                                                                                                            gCGGGTGCGAGGTCGTGCTGCGTCTGTGTTATCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTc           >  2:804126/1‑150 (MQ=255)
                                                                                                                                              gggTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTg         <  2:41137/150‑1 (MQ=255)
                                                                                                                                                 tGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTAc      <  1:684001/150‑1 (MQ=255)
                                                                                                                                                    gAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGa   >  2:2854978/1‑150 (MQ=255)
                                                                                                                                                    gAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGa   >  1:1397631/1‑150 (MQ=255)
                                                                                                                                                    gAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGa   <  2:698820/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACAGAAGCGGTCCATCCAGGCGAAGATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGAT  >  NC_000913/265032‑265330

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑