Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,171,904 C→G 5.4% E324Q (GAA→CAA)  gatD ← galactitol‑1‑phosphate 5‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,171,9040CG 5.4% 312.2 / 5.3 91E324Q (GAA→CAA) gatDgalactitol‑1‑phosphate 5‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (38/48);  new base G (2/3);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TAGTGTAAACGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTACGCAATATTTTGC  >  NC_000913/2171755‑2172053
                                                                                                                                                     |                                                                                                                                                     
tAGTGTAAACGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTc                                                                                                                                                       <  1:87700/150‑1 (MQ=255)
tAGTGTAAACGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTc                                                                                                                                                       <  1:689837/150‑1 (MQ=255)
tAGTGTAAACGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTc                                                                                                                                                       <  2:529792/150‑1 (MQ=255)
  gtgtAAACGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCaa                                                                                                                                                     >  2:70110/1‑150 (MQ=255)
   tgtAAACGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCaaa                                                                                                                                                    <  1:990347/150‑1 (MQ=255)
    gtAAACGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAg                                                                                                                                                   >  2:35367/1‑150 (MQ=255)
    gtAAACGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAg                                                                                                                                                   <  1:1326786/150‑1 (MQ=255)
     tAAACGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGc                                                                                                                                                  >  2:810477/1‑150 (MQ=255)
      aaaCGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCt                                                                                                                                                 >  2:292957/1‑150 (MQ=255)
      aaaCGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCt                                                                                                                                                 >  2:1434771/1‑150 (MQ=255)
       aaCGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCtt                                                                                                                                                <  2:402562/150‑1 (MQ=255)
        aCGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCAAAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTc                                                                                                                                               >  1:305005/1‑150 (MQ=255)
         cGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTcc                                                                                                                                              <  1:1464107/150‑1 (MQ=255)
           aaGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTGAAAGCTTcgcc                                                                                                                                            <  2:1438969/150‑4 (MQ=255)
             gggggACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACgg                                                                                                                                          >  2:11028/1‑150 (MQ=255)
                ggACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGa                                                                                                                                       <  1:22682/150‑1 (MQ=255)
                 gACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAg                                                                                                                                      <  1:1535035/150‑1 (MQ=255)
                      ttGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGAtt                                                                                                                                 <  1:70110/150‑1 (MQ=255)
                      ttGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGAtt                                                                                                                                 <  1:1930132/150‑1 (MQ=255)
                        ggTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCATGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTGAAAGCTTcgccc                                                                                                                                           >  2:1458/1‑134 (MQ=255)
                        ggTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCATGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTGAAAGCTTcgccc                                                                                                                                           <  1:1458/138‑5 (MQ=255)
                         gTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAAt                                                                                                                              <  1:35367/150‑1 (MQ=255)
                          tGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATg                                                                                                                             >  1:617854/1‑150 (MQ=255)
                            aaTGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATTTCACGCACCGCCTGGGAGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGc                                                                                                                           <  1:1516546/150‑1 (MQ=255)
                            aaTGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGc                                                                                                                           <  2:1752120/150‑1 (MQ=255)
                               gATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTcc                                                                                                                        >  1:1777063/1‑150 (MQ=255)
                               gATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTcc                                                                                                                        <  2:1224291/150‑1 (MQ=255)
                                   tGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGc                                                                                                                    >  2:1258415/1‑150 (MQ=255)
                                    tttttGCACAATACCGGGCCAGCAGCACGCTGTCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATATCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTGAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCt                                                                                                                   <  1:347986/149‑1 (MQ=255)
                                      tttGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTa                                                                                                                 <  2:542173/150‑1 (MQ=255)
                                         gCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACt                                                                                                              >  1:897379/1‑150 (MQ=255)
                                         gCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACt                                                                                                              <  2:1261509/150‑1 (MQ=255)
                                           acaATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTa                                                                                                            <  1:1470464/150‑1 (MQ=255)
                                              aTACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGt                                                                                                         <  1:579488/150‑1 (MQ=255)
                                                aCCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTc                                                                                                       <  2:384454/150‑1 (MQ=255)
                                                aCCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTc                                                                                                       <  2:1767844/150‑1 (MQ=255)
                                                 ccGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCATCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCt                                                                                                      >  2:1800508/1‑150 (MQ=255)
                                                 ccGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCt                                                                                                      >  2:1755974/1‑150 (MQ=255)
                                                 ccGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCt                                                                                                      >  2:1925827/1‑150 (MQ=255)
                                                  cGGGCCAGCGTGATGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTg                                                                                                     >  2:1322317/1‑150 (MQ=255)
                                                  cGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTg                                                                                                     >  1:852/1‑150 (MQ=255)
                                                      ccAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAg                                                                                                 <  1:1288886/150‑1 (MQ=255)
                                                      ccAGCATCAAGCTGGCCCACCGTTTCAGCGAATGAGCAACACTTTGCCCCCCATATCATTAAGAACCATGTCACGCACCGCCAGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAg                                                                                                 >  2:1851639/1‑150 (MQ=255)
                                                       cAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGc                                                                                                >  2:787901/1‑150 (MQ=255)
                                                        aGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCa                                                                                               <  1:292957/150‑1 (MQ=255)
                                                           aTCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAAcc                                                                                            <  2:193239/150‑1 (MQ=255)
                                                                 cTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGa                                                                                                                                       >  1:1445424/1‑101 (MQ=255)
                                                                 cTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGa                                                                                                                                       <  2:1445424/101‑1 (MQ=255)
                                                                    gCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCg                                                                                   >  2:423342/1‑150 (MQ=255)
                                                                        gcgGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTc                                                                               <  2:1097070/150‑1 (MQ=255)
                                                                            tttCAGGGAATGAGCAACACTTTGCCCGGCATAGCGTTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCAc                                                                           >  2:1464630/1‑150 (MQ=255)
                                                                               cAGGGTATGAGCACCACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTcc                                                                        >  2:1351441/1‑150 (MQ=255)
                                                                               cAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTcc                                                                        >  2:1801882/1‑150 (MQ=255)
                                                                               cAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTcc                                                                        <  1:344344/150‑1 (MQ=255)
                                                                               cAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTcc                                                                        >  2:1809127/1‑150 (MQ=255)
                                                                                 gggAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACt                                                                                                              >  1:1551851/1‑110 (MQ=255)
                                                                                 gggAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACt                                                                                                              <  2:1551851/110‑1 (MQ=255)
                                                                                        aGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGc                                                               <  1:305823/150‑1 (MQ=255)
                                                                                        aGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGc                                                               <  1:298611/150‑1 (MQ=255)
                                                                                         gCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGcc                                                              >  1:205219/1‑150 (MQ=255)
                                                                                          cAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCa                                                             >  2:330602/1‑150 (MQ=255)
                                                                                           aaCACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGTCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTCTCAGCAACCGGCTCGCCGTTTCCCACTCCTCACCCGGCCaa                                                            >  2:1967432/1‑150 (MQ=255)
                                                                                               cTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGc                                                        >  2:1854077/1‑150 (MQ=255)
                                                                                                  tGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTgg                                                     <  1:1211982/150‑1 (MQ=255)
                                                                                                   gCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGa                                                    >  2:239358/1‑150 (MQ=255)
                                                                                                   gCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGa                                                    <  1:810477/150‑1 (MQ=255)
                                                                                                    cccGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGag                                                   <  1:423342/150‑1 (MQ=255)
                                                                                                        gCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGagtagt                                               >  1:329544/1‑150 (MQ=255)
                                                                                                          aTAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTc                                             <  2:329544/150‑1 (MQ=255)
                                                                                                           tAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCACAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCa                                            <  1:1002624/150‑1 (MQ=255)
                                                                                                             gCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATc                                          <  1:365562/150‑1 (MQ=255)
                                                                                                              cATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATcc                                         <  1:787901/150‑1 (MQ=255)
                                                                                                              cATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATcc                                         <  2:207715/150‑1 (MQ=255)
                                                                                                               aTTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCa                                        <  1:261537/150‑1 (MQ=255)
                                                                                                                  aCGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACt                                     >  1:1141162/1‑150 (MQ=255)
                                                                                                                   cGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTg                                    >  2:1532452/1‑150 (MQ=255)
                                                                                                                     aGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGcc                                  <  1:1322317/150‑1 (MQ=255)
                                                                                                                     aGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGcc                                  <  1:1755974/150‑1 (MQ=255)
                                                                                                                     aGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGcc                                  <  1:239358/150‑1 (MQ=255)
                                                                                                                        gATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGAt                               <  2:1141162/150‑1 (MQ=255)
                                                                                                                         aTGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATa                              <  2:384066/150‑1 (MQ=255)
                                                                                                                           gTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAAc                            <  2:897379/150‑1 (MQ=255)
                                                                                                                            tCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTc                                             >  1:34772/1‑132 (MQ=255)
                                                                                                                            tCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTc                                             <  2:34772/132‑1 (MQ=255)
                                                                                                                               cGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTc                        <  2:1777063/150‑1 (MQ=255)
                                                                                                                               cGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTc                        <  2:205219/150‑1 (MQ=255)
                                                                                                                                gCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCa                       <  1:1434771/150‑1 (MQ=255)
                                                                                                                                  aCCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGc                     >  2:63173/1‑150 (MQ=255)
                                                                                                                                     gCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGctct                  >  2:221945/1‑150 (MQ=255)
                                                                                                                                         gggCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTAc              >  2:1224634/1‑150 (MQ=255)
                                                                                                                                          gtggAAGCTTTGAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTACg             >  2:200174/4‑150 (MQ=255)
                                                                                                                                           gCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTACGc            <  1:65009/150‑1 (MQ=255)
                                                                                                                                                 cTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTACGCAATAtt      >  2:577864/1‑150 (MQ=255)
                                                                                                                                                    tCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTACGCAATATTTTg   <  1:1258415/150‑1 (MQ=255)
                                                                                                                                                     cAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTACGCAATATTTTGc  <  1:577864/150‑1 (MQ=255)
                                                                                                                                                     cAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTACGCAATATTTTGc  <  1:1925827/150‑1 (MQ=255)
                                                                                                                                                     cAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTACGCAATATTTTGc  <  1:1854077/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAGTGTAAACGAAGGGGGACCATTGGTGAATGATCTGTTTTGCACAATACCGGGCCAGCATCACGCTGGCCCGCGGTTTCAGGGAATGAGCAACACTTTGCCCGGCATAGCATTACGAGCGATGTCACGCACCGCCTGGGCGAAGCTTTCAAAGCTTCCACGGTGAGCGATTAATGGCTCCAGGCTTAACTTACGTTCTGTCAGCAACCGGCTCGCCGTTTCCCACTCCTGCCCCGGCCAAGGGCTGGAGTAGTTCATCCAACTGCCGATAACCGTCAGCTCTTTACGCAATATTTTGC  >  NC_000913/2171755‑2172053

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