Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,105 G→C 7.0% W204C (TGG→TGC) ‡ sseA → 3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,653,1050GC 7.0% 278.3 / 7.8 83W204C (TGG→TGC) ‡sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (38/39);  new base C (4/2);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

GTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGAC  >  NC_000913/2652956‑2653251
                                                                                                                                                     |                                                                                                                                                  
gtgAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTgg                                                                                                                                                    <  1:56978/150‑1 (MQ=255)
gtgAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTgg                                                                                                                                                    >  1:57245/1‑150 (MQ=255)
    aaGTAACCGATGAATTATTGGCAAGCCATGAAAATACGGAGCAAATAATTGATGCCCGCAATGCAGCACGTATTAACGCACTAGTTGATGAACCTCGCCCAGGATTACGTCGCGGACATAATCCCGGAGCACAGAATCTTCAGTGCACgg                                                                                                                                                >  1:1838142/1‑150 (MQ=255)
     aGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGa                                                                                                                                               >  1:307351/1‑150 (MQ=255)
             aTGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattcag                                                                                                                                       <  2:1717155/150‑7 (MQ=255)
              tGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattcagt                                                                                                                                      <  1:815506/150‑8 (MQ=255)
                tattatTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGATTTAACGAAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGcgcg                                                                                                                                    >  2:856947/1‑150 (MQ=255)
                 attattGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGa                                                                                                                                   >  2:97681/1‑150 (MQ=255)
                    attGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAgg                                                                                                                                >  1:357729/1‑150 (MQ=255)
                     ttGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGc                                                                                                                               <  2:901496/150‑1 (MQ=255)
                      tGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCg                                                                                                                              >  1:282589/1‑150 (MQ=255)
                      tGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCACGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCg                                                                                                                              >  1:1979091/1‑150 (MQ=255)
                       ggCAAGCCATGAAAATATGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACACATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGa                                                                                                                             >  2:534177/1‑150 (MQ=255)
                       ggCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGa                                                                                                                             >  1:1513862/1‑150 (MQ=255)
                       ggCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTGGATGAACCTCGCCCAGGTTTACGGCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGa                                                                                                                             <  1:615482/150‑1 (MQ=255)
                        gCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGaa                                                                                                                            <  2:167305/150‑1 (MQ=255)
                          aaGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAAGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACt                                                                                                                          <  1:304335/150‑1 (MQ=255)
                                   aaaTACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGAc                                                                                                                 <  2:771207/150‑1 (MQ=255)
                                    aaTACGTCTAAAATTATTGATTAACACCCTGATGCACTATTTAAAGCAGAAATTAATTAAAATAAAACAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGAcc                                                                                                                >  2:840231/1‑150 (MQ=255)
                                    aaTACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGAcc                                                                                                                <  1:1628800/150‑1 (MQ=255)
                                         ggCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGa                                                                                                           >  1:1403904/1‑150 (MQ=255)
                                             cAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTg                                                                                                       >  2:864378/1‑150 (MQ=255)
                                              aaaTTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCGCGGAGCACTGAATGTGCCGTGGACGGAACGGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTgg                                                                                                      <  2:1120477/150‑1 (MQ=255)
                                                     tGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGat                                                                                               >  1:1964170/1‑150 (MQ=255)
                                                         gcccgcccGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAttt                                                                                           >  2:1437722/1‑150 (MQ=255)
                                                         gcccgcccGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAttt                                                                                           >  2:505405/1‑150 (MQ=255)
                                                          cccgcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAtttt                                                                                          >  1:806694/1‑150 (MQ=255)
                                                            cgcccgGCTGCACGTTTTAACGCATTAGTTGATGAACCTCGCCCAGGTTTAGGTCGCGGACTTATTCCCGGAGCACTGAATGTTCCGTGGAGGGATCTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAtttttt                                                                                        <  1:1689472/150‑1 (MQ=255)
                                                            cgcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAtttttt                                                                                        <  2:978086/150‑1 (MQ=255)
                                                            cgcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAtttttt                                                                                        >  2:613173/1‑150 (MQ=255)
                                                             gcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTg                                                                                       <  2:903599/150‑1 (MQ=255)
                                                             gcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTg                                                                                       <  2:1652282/150‑1 (MQ=255)
                                                             gcccgGCTGCACGTTTTAACGAAGAATTTTATGAACCTCGCCCAGGTTTACTTCTCGTAAATATTCCCGTAGCACTGAATGTTCCGTTGACGGAACTGTTGCGCTAAGGCGAACTAAAAACTACCTATGAACTGGATGCGATATTTTTTg                                                                                       <  1:471964/150‑1 (MQ=255)
                                                               ccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGGTCCGTTGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGt                                                                                     <  1:1831715/150‑1 (MQ=255)
                                                                     gCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCgcggc                                                                               >  2:355945/1‑150 (MQ=255)
                                                                      cACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCgcggcg                                                                              >  2:333042/1‑150 (MQ=255)
                                                                       aCGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGt                                                                             <  2:1065462/150‑1 (MQ=255)
                                                                        cGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTc                                                                            >  1:850391/1‑150 (MQ=255)
                                                                        cGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTc                                                                            >  1:743248/1‑150 (MQ=255)
                                                                         gTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCa                                                                           >  1:822973/1‑150 (MQ=255)
                                                                                 gCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACgcc                                                                   >  1:1813151/1‑148 (MQ=255)
                                                                                 gCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGAc                                                                   <  1:260868/150‑1 (MQ=255)
                                                                                      aGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGc                                                                                                  >  2:222371/1‑114 (MQ=255)
                                                                                      aGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGc                                                                                                  <  1:222371/114‑1 (MQ=255)
                                                                                      aGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAAcc                                                              <  1:972100/150‑1 (MQ=255)
                                                                                        ttGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCaa                                                            <  1:1477636/150‑1 (MQ=255)
                                                                                         tgatgaACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAAt                                                           >  1:1300591/1‑150 (MQ=255)
                                                                                              aaCCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCg                                                      <  1:1866426/150‑1 (MQ=255)
                                                                                              aaCCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCg                                                      <  2:1964170/150‑1 (MQ=255)
                                                                                               aCCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGt                                                     >  2:1882062/1‑150 (MQ=255)
                                                                                                ccTCGCCCAGGTTTACGTCGCGGACACATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTc                                                    <  1:534177/150‑1 (MQ=255)
                                                                                                         ggTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTgg                                                                                      <  1:336737/107‑1 (MQ=255)
                                                                                                         ggTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTgg                                                                                      >  2:336737/1‑107 (MQ=255)
                                                                                                         ggTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGc                                           <  1:1580880/150‑1 (MQ=255)
                                                                                                          gTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCTTCAGCTGCGGct                                          >  1:818260/1‑150 (MQ=255)
                                                                                                             tACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTg                                       >  2:1213643/1‑150 (MQ=255)
                                                                                                              aCGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTgg                                      <  1:840231/150‑1 (MQ=255)
                                                                                                                gTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGgtg                                    <  2:307351/150‑1 (MQ=255)
                                                                                                                 tCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGgtgt                                   >  2:1215988/1‑150 (MQ=255)
                                                                                                                      gACATATTCCCGGAGCACTGAATGATCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGGCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACgg                              <  1:881413/150‑1 (MQ=255)
                                                                                                                        cATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCa                            <  2:1300591/150‑1 (MQ=255)
                                                                                                                          tatTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGc                          >  1:166844/1‑150 (MQ=255)
                                                                                                                          tatTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGc                          >  2:1516544/1‑150 (MQ=255)
                                                                                                                            ttCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGAGCGAAGGCGAACTAAAAACGACCGATGAAGTGGATCCGATATTTTTTGGTCGCGGCGACAGCTACGACAAACAAATTATCGTCAGCTGCGGCTCGGGTGTAACGGCAGCCg                        >  1:1461463/1‑150 (MQ=255)
                                                                                                                              cccGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGgtg                                    <  2:992219/136‑1 (MQ=255)
                                                                                                                              cccGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGgtg                                    >  1:992219/1‑136 (MQ=255)
                                                                                                                                cGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGt                    <  1:613173/150‑1 (MQ=255)
                                                                                                                                 ggAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGt                    >  1:812046/1‑149 (MQ=255)
                                                                                                                                 ggAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGt                    <  2:812046/149‑1 (MQ=255)
                                                                                                                                    gCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTg                <  2:822973/150‑1 (MQ=255)
                                                                                                                                    gCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTg                <  2:882516/150‑1 (MQ=255)
                                                                                                                                     cACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGt               <  1:97681/150‑1 (MQ=255)
                                                                                                                                      aCTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGtt              >  2:536379/1‑150 (MQ=255)
                                                                                                                                      aCTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGtt              <  2:266508/150‑1 (MQ=255)
                                                                                                                                      aCTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGtt              >  1:267733/1‑150 (MQ=255)
                                                                                                                                      aCTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGtt              <  2:1403904/150‑1 (MQ=255)
                                                                                                                                      aCTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGtt              <  2:166844/150‑1 (MQ=255)
                                                                                                                                       gcaccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTa             >  2:1989402/7‑150 (MQ=255)
                                                                                                                                        caccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAg            >  2:458489/6‑150 (MQ=255)
                                                                                                                                        caccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAg            >  1:1350685/6‑150 (MQ=255)
                                                                                                                                             gTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTc       >  1:1315951/1‑150 (MQ=255)
                                                                                                                                               tCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTcgc     <  2:1042480/150‑1 (MQ=255)
                                                                                                                                                 cGTGGACGGAACTGTTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGa   >  2:1584366/1‑150 (MQ=255)
                                                                                                                                                 cGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGa   <  1:1213643/150‑1 (MQ=255)
                                                                                                                                                  gTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGAc  <  1:355945/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGAC  >  NC_000913/2652956‑2653251

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: