Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,919,903 C→G 5.1% G282R (GGC→CGC)  gudX ← glucarate dehydratase‑related protein, substrate unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,919,9030CG 5.1% 275.5 / 4.8 93G282R (GGC→CGC) gudXglucarate dehydratase‑related protein, substrate unknown
Reads supporting (aligned to +/- strand):  ref base C (41/47);  new base G (2/3);  total (43/50)
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.96e-01

AGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAAC  >  NC_000913/2919756‑2920050
                                                                                                                                                   |                                                                                                                                                   
aGTCGTCGCAAAGCTGCGCCACCCGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGCAGTATCTCCCCCATTGAGCATCCCCGCATGTCCCAATTCGCGCCAGGTGGTGGCGAGGATGTGAGTCGCGTCGGGCCAGCCgg                                                                                                                                                   >  1:289364/1‑150 (MQ=255)
    gtcgCAAAGCTGCGCCACACGGACTGCACCGGATAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGc                                                                                                                                               >  1:346411/1‑150 (MQ=255)
    gtcgCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGc                                                                                                                                               >  1:811704/1‑150 (MQ=255)
     tcgCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACCGGCTTGCCCGTCGcg                                                                                                                                              <  2:813190/150‑2 (MQ=255)
     tcgCAAAGCGGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGcc                                                                                                                                              <  2:1180699/150‑1 (MQ=255)
       gCAAATCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCTTGACCGATTTCGCGACATTTGGTGGCGTTCTTGTTAGTCGCGACGGGCAAGCCGGTCGCCCg                                                                                                                                            <  2:462032/150‑1 (MQ=255)
       gCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGGATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCg                                                                                                                                            <  2:189438/150‑1 (MQ=255)
       gCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATAGCTTCCGCATGGAGCATCACCGCATGACCCATTTCGCGCCTGGTGGTGGCGGTCAGGTTAGTCGCGCCGGGCCAGCCGGTCGCCCg                                                                                                                                            >  2:472462/1‑150 (MQ=255)
         aacgCTGCGCCACCCGGACTGCACCGGAAAGCGTCCAGAAGAGCGGATCCGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTc                                                                                                                                          <  1:656127/147‑1 (MQ=255)
          aaGCTGCGCCACACGGACTGCACCGGAAACCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGAGCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCGACCGGCTTGCCCGTCg                                                                                                                                         <  1:1294697/150‑1 (MQ=255)
            gCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCAGCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCgta                                                                                                                                       >  2:17514/1‑148 (MQ=255)
            gCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCGACCGGCTTGCCCGTCgcg                                                                                                                                       <  1:1267736/150‑3 (MQ=255)
                cgcCACACGGACTGCACCGGACAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTc                                                                                                                                   <  1:624653/150‑1 (MQ=255)
                cgcCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGGGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTc                                                                                                                                   <  1:1154326/150‑1 (MQ=255)
                   cacaCGGACTGCACCGGATAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGc                                                                                                                                <  2:346411/150‑1 (MQ=255)
                      cccgACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGAGGAATATCTACCGCAATGAGCATCACCGCATGACCCTTTTCGCGCCAGTATGTGGCGCTCATGTTAGTCGCGACGGTCAAGCCGGTCGCCCGTCGAAATTCCggaat                                                                                                                             <  2:220634/147‑5 (MQ=255)
                       cGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATc                                                                                                                            <  2:333460/150‑1 (MQ=255)
                       cGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATc                                                                                                                            <  2:1103343/150‑1 (MQ=255)
                             gCACCGGAAAGCGTCCAGACGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCa                                                                                                                      <  1:295400/150‑1 (MQ=255)
                             gCACCGGAAAGCGTCCAGAAGGGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCa                                                                                                                      <  2:1249601/150‑1 (MQ=255)
                               aCCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACg                                                                                                                    >  1:522928/1‑150 (MQ=255)
                                ccGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGt                                                                                                                   >  2:1096872/1‑150 (MQ=255)
                                  ggAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTcc                                                                                                                 >  2:529674/1‑150 (MQ=255)
                                     aaGCGTCCAGAAGTTCGGATCGGCAAGTGGAATATCTACGGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCTATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGga                                                                                                              <  1:827810/150‑1 (MQ=255)
                                     aaGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGga                                                                                                              >  1:509778/1‑150 (MQ=255)
                                      aGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGgag                                                                                                             <  1:46649/150‑1 (MQ=255)
                                       gCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGgaga                                                                                                            <  2:404911/150‑1 (MQ=255)
                                          tCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGc                                                                                                         <  2:739538/150‑1 (MQ=255)
                                          tCCAGAAGTGCGGATCGGCAAGGGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTTGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGc                                                                                                         <  2:704472/150‑1 (MQ=255)
                                                 gTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCTCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCGGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTc                                                                                                  <  1:849810/150‑1 (MQ=255)
                                                 gTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTc                                                                                                  >  2:535706/1‑150 (MQ=255)
                                                   gCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTg                                                                                                <  1:1285890/150‑1 (MQ=255)
                                                   gCGGATCGGCAAGTGGAAATTCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCAGGTCCGGAGAAGCCCTGTTCTg                                                                                                <  2:1039369/150‑1 (MQ=255)
                                                   gCGGATCGGCAACTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTg                                                                                                <  2:1044734/150‑1 (MQ=255)
                                                    cGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTgc                                                                                               >  2:1121634/1‑150 (MQ=255)
                                                         cGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGACGCCCTGGTCTGcgccgc                                                                                          >  2:803131/1‑150 (MQ=255)
                                                          ggCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCa                                                                                         >  2:113783/1‑150 (MQ=255)
                                                          ggCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTAACGTCCGGAGAAGCCCTGTTCTGCGCCGCa                                                                                         >  2:584846/1‑150 (MQ=255)
                                                          ggCAAGTGGAATATCTACCGCATTGAGCATCACCGAATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCTACAGGCAAGCCGCTCCCCCGTCGAAATTCCGCCATCACATCACGTCCGGAGCAGCCAGGGTCTGCGCCGCa                                                                                         >  2:252194/1‑150 (MQ=255)
                                                             aaGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATgg                                                                                      >  2:567319/1‑150 (MQ=255)
                                                             aaGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCAGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATgg                                                                                      >  1:15086/1‑150 (MQ=255)
                                                             aaGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTGGGGGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATgg                                                                                      <  2:171325/150‑1 (MQ=255)
                                                               gTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTGGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGAt                                                                                    <  1:977444/150‑1 (MQ=255)
                                                                  gAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCtt                                                                                 >  1:594293/1‑150 (MQ=255)
                                                                     tatCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCgg                                                                              <  2:465591/150‑1 (MQ=255)
                                                                        cTACCGCATTGCGCATCACCGCATGACCCATTTCGCGCCAGTTGCTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGGCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATTGATCTTCGGCAt                                                                           <  2:251150/150‑1 (MQ=255)
                                                                         tACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGATATTCCGCCATCACTTCAGGTCCGTAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATa                                                                          <  2:733425/150‑1 (MQ=255)
                                                                             gCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTa                                                                      >  1:112059/1‑150 (MQ=255)
                                                                                ttGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTTGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGa                                                                   <  1:813645/150‑1 (MQ=255)
                                                                                ttGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGa                                                                   >  1:274897/1‑150 (MQ=255)
                                                                                   aGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCCCGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAAca                                                                >  2:506936/1‑150 (MQ=255)
                                                                                    gCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGACcat                                                               >  1:370867/1‑150 (MQ=255)
                                                                                    gCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAAcat                                                               >  1:352695/1‑150 (MQ=255)
                                                                                     cATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGGCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGCCATAGGTAAGAAcatc                                                              >  1:480059/1‑150 (MQ=255)
                                                                                      attaCCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCGCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAAcatca                                                             <  2:299667/147‑1 (MQ=255)
                                                                                          ccGCATGACCCATTTAGCGCCAGTTGGTTGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCTTGGATCTTCGGCATAGGTAAGAACATCATTCa                                                         <  1:198200/150‑1 (MQ=255)
                                                                                           cGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAg                                                        >  1:366562/1‑150 (MQ=255)
                                                                                           cGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAg                                                        >  1:850982/1‑150 (MQ=255)
                                                                                               tGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCt                                                    <  2:33151/150‑1 (MQ=255)
                                                                                                  cccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCgttt                                                 >  1:1264023/1‑149 (MQ=255)
                                                                                                   ccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGc                                                <  1:113783/150‑1 (MQ=255)
                                                                                                   ccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGc                                                <  1:1079508/150‑1 (MQ=255)
                                                                                                    cATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCa                                               >  2:1154946/1‑150 (MQ=255)
                                                                                                    cATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTAACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTAAGCCCTTTGCa                                               >  2:1162074/1‑150 (MQ=255)
                                                                                                         cgcgCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAg                                          <  1:943098/150‑1 (MQ=255)
                                                                                                         cgcgCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAg                                          <  1:529674/150‑1 (MQ=255)
                                                                                                         cgcgCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAg                                          <  1:135665/150‑1 (MQ=255)
                                                                                                           cgcCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGaa                                        >  1:452813/1‑150 (MQ=255)
                                                                                                             ccAGTTGGTGGCGATCATGTTAGTCGCGACGGTCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAAt                                      <  2:316194/150‑1 (MQ=255)
                                                                                                             ccAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAAt                                      >  1:1232296/1‑150 (MQ=255)
                                                                                                              cAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATg                                     >  1:1249730/1‑150 (MQ=255)
                                                                                                              cAGTTGGTGGCGAACATGTTAGTCCCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCAAGGATCTTCGGCATAGGTAAGACCATCATTCAGCCCTTTGCATAAAGAAATg                                     >  2:964714/1‑150 (MQ=255)
                                                                                                                   ggtggCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTtc                                >  1:387604/1‑150 (MQ=255)
                                                                                                                   ggtggCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTtc                                <  2:710496/150‑1 (MQ=255)
                                                                                                                    gtggCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTtca                               <  2:274897/150‑1 (MQ=255)
                                                                                                                       gCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGATCATCATTCGGCCCTTTGCATAAAGAATTGGCTtcatca                            >  1:90701/1‑150 (MQ=255)
                                                                                                                         gATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAag                          <  1:1000316/150‑1 (MQ=255)
                                                                                                                         gATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGGAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAag                          >  2:486223/1‑150 (MQ=255)
                                                                                                                           tCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAagca                        <  1:1069705/150‑1 (MQ=255)
                                                                                                                              tGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCTGCCCTTTGCATAAAGAAATGGCTTCATCCAGCAGcc                     >  2:853513/1‑150 (MQ=255)
                                                                                                                               gTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAACTTCCGCCCTCACTTCACGTCCGGAGACGCCCTGTTCTGCGCCGCATGGAGCTGCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGGTTCATCAAGCAGCCa                    >  1:820036/1‑150 (MQ=255)
                                                                                                                                 tAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATg                  <  1:795036/150‑1 (MQ=255)
                                                                                                                                     ttcgACGGGCGACCGGCTTGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCTTAGGTAAGAACATCATTCAGCCCGTTGCATAAAGAAATGGCTTCATCAAGCCGCCATGCAcc              >  2:113090/3‑150 (MQ=255)
                                                                                                                                     cgcgACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCAcc              >  1:1066180/1‑150 (MQ=255)
                                                                                                                                       cgACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGt            <  2:1189701/150‑1 (MQ=255)
                                                                                                                                        gACGGGCGACCGGCTTGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGtt           >  1:560640/1‑150 (MQ=255)
                                                                                                                                            ggCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGa       <  1:970637/150‑1 (MQ=255)
                                                                                                                                             gCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTGGGGAt      >  1:1209546/1‑150 (MQ=255)
                                                                                                                                              cAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATc     <  1:1096872/150‑1 (MQ=255)
                                                                                                                                              cAAGCCGGTCGCCCGTCGAAATTCCGCCATCACATCACGTCCGGAGAAGCCCTGTTCTGCGCCGCAAGGATCTTCGGCATAGGTAAGCACATCATTCCGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATc     >  1:242564/1‑150 (MQ=255)
                                                                                                                                                aGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCaa   <  1:359754/150‑1 (MQ=255)
                                                                                                                                                 gCCGGTCGCCCGTCGAACTTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAAc  <  2:625338/150‑1 (MQ=255)
                                                                                                                                                 gCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAAc  <  2:596499/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAAC  >  NC_000913/2919756‑2920050

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