Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,104 G→C 6.0% W204S (TGG→TCG) ‡ sseA → 3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,653,1040GC 6.0% 271.1 / 6.9 83W204S (TGG→TCG) ‡sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (39/37);  new base C (3/2);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

GGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGAC  >  NC_000913/2652955‑2653251
                                                                                                                                                     |                                                                                                                                                   
ggtgAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTg                                                                                                                                                     >  1:167305/1‑150 (MQ=255)
 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)
                                                                                                                                                     |                                                                                                                                                   
GGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGAC  >  NC_000913/2652955‑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: