Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,919,905 G→C 5.1% T281S (ACC→AGC)  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,9050GC 5.1% 272.8 / 4.9 92T281S (ACC→AGC) gudXglucarate dehydratase‑related protein, substrate unknown
Reads supporting (aligned to +/- strand):  ref base G (39/48);  new base C (3/2);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

AGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGT  >  NC_000913/2919756‑2920053
                                                                                                                                                     |                                                                                                                                                    
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:1103343/150‑1 (MQ=255)
                       cGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATc                                                                                                                               <  2:333460/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:850982/1‑150 (MQ=255)
                                                                                           cGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAg                                                           >  1:366562/1‑150 (MQ=255)
                                                                                               tGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCt                                                       <  2:33151/150‑1 (MQ=255)
                                                                                                  cccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCgttt                                                    >  1:1264023/1‑149 (MQ=255)
                                                                                                   ccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGc                                                   <  1:1079508/150‑1 (MQ=255)
                                                                                                   ccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGc                                                   <  1:113783/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)
                                                                                                                                                    ggTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATATGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGt  <  1:803131/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGT  >  NC_000913/2919756‑2920053

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: