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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,919,9030CG 7.3% 262.6 / 9.2 82G282R (GGC→CGC) ‡gudXglucarate dehydratase‑related protein
Reads supporting (aligned to +/- strand):  ref base C (35/41);  new base G (4/2);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

CGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGG  >  NC_000913/2919759‑2920052
                                                                                                                                                |                                                                                                                                                     
cgtcgCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCg                                                                                                                                                  >  1:913793/1‑150 (MQ=255)
 gtcgCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGc                                                                                                                                                 >  2:1369841/1‑150 (MQ=255)
   cgCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGccc                                                                                                                                               >  1:1022262/1‑150 (MQ=255)
     cAAAGCTGCGCCCCACGGACTGCACCGGAGAGCGTCCAGAAGTGAGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGt                                                                                                                                             >  2:74847/1‑150 (MQ=255)
       aaGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTCCCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGGCGCGACGGGCAAGCCGGTCGCCCGTCg                                                                                                                                           <  2:691427/150‑1 (MQ=255)
       aaGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGTTGGCGATCATTTTAGTCTCGACGGGCGACCGGCTTGCCCGTCg                                                                                                                                           <  2:754835/150‑1 (MQ=255)
       aaGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAATCCGGGCGCCCGTCg                                                                                                                                           <  2:696960/150‑1 (MQ=255)
         gCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGaa                                                                                                                                         >  2:295605/1‑150 (MQ=255)
            gcgcCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCGACCGGCTTGCCCGTCgcgact                                                                                                                                      <  2:549882/150‑6 (MQ=255)
              gcCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTcc                                                                                                                                    <  2:789767/150‑1 (MQ=255)
              gcCACACGGACGGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGTCGGGCAAGCCGTTCGCCCGTCGATATTcg                                                                                                                                    <  2:1296325/150‑2 (MQ=255)
               ccACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCg                                                                                                                                   >  1:884488/1‑150 (MQ=255)
                cacaCGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGc                                                                                                                                  >  2:446078/1‑150 (MQ=255)
                      gACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCAc                                                                                                                            <  2:734755/150‑1 (MQ=255)
                         tGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTc                                                                                                                         <  1:1307829/150‑1 (MQ=255)
                          gCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCa                                                                                                                        <  2:685949/150‑1 (MQ=255)
                          gCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCa                                                                                                                        <  2:345641/150‑1 (MQ=255)
                               ggcaaGCGTCCAGAAGTGCGGATCGGCAAGTCGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTcc                                                                                                                   <  1:1289948/147‑1 (MQ=255)
                                gAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCg                                                                                                                  <  1:81857/150‑1 (MQ=255)
                                gAAAGCGGCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCg                                                                                                                  <  2:216889/150‑1 (MQ=255)
                                         cAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGccc                                                                                                         <  1:680735/150‑1 (MQ=255)
                                         cAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGccc                                                                                                         >  2:688847/1‑150 (MQ=255)
                                                   gATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGcc                                                                                                                                                       <  1:861582/94‑1 (MQ=255)
                                                   gATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGcc                                                                                                                                                       >  2:861582/1‑94 (MQ=255)
                                                     tCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGcgccg                                                                                             >  2:1113917/1‑150 (MQ=255)
                                                      cGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCa                                                                                                                             >  2:1017015/1‑117 (MQ=255)
                                                      cGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCa                                                                                                                             <  1:1017015/117‑1 (MQ=255)
                                                      cGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGcgccgc                                                                                            >  1:1244554/1‑150 (MQ=255)
                                                       ggCCAGTGGAATATCTACCGCATTGCGACTCACCGCATGACCCAGTTCGCGCCAGTTGGTGGGGAGCAGGTTAGTCGCGACGGGCAAGCCGGGCGCCCGTCGAAATTCCGCCCTCCCTTCACCTCCGGAGACGCCCTGTTCTGcgccgcc                                                                                           >  2:609114/1‑149 (MQ=255)
                                                       ggCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCa                                                                                                                        >  2:906192/1‑121 (MQ=255)
                                                       ggCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCa                                                                                                                        <  1:906192/121‑1 (MQ=255)
                                                       ggCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCa                                                                                           >  1:1109598/1‑150 (MQ=255)
                                                        gCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCAt                                                                                          >  1:994958/1‑150 (MQ=255)
                                                           aGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGa                                                                                       <  1:1073525/150‑1 (MQ=255)
                                                              ggAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCt                                                                                    >  1:1113315/1‑150 (MQ=255)
                                                               gAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCtt                                                                                   >  1:311767/1‑150 (MQ=255)
                                                               gAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCtt                                                                                   >  2:478148/1‑150 (MQ=255)
                                                                aaTATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTc                                                                                  <  2:461396/150‑1 (MQ=255)
                                                                 atatCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCg                                                                                 <  2:1113315/150‑1 (MQ=255)
                                                                   atCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGc                                                                               <  1:1216026/150‑1 (MQ=255)
                                                                   atCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGc                                                                               >  1:1230726/1‑150 (MQ=255)
                                                                   atCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGc                                                                               <  2:331919/150‑1 (MQ=255)
                                                                       aCCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAg                                                                           >  1:1095017/1‑150 (MQ=255)
                                                                        ccGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAgg                                                                          <  1:365329/150‑1 (MQ=255)
                                                                        ccGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAgg                                                                          <  1:1169483/150‑1 (MQ=255)
                                                                          gCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTa                                                                        >  1:1495982/1‑150 (MQ=255)
                                                                          gCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTa                                                                        >  1:1280589/1‑150 (MQ=255)
                                                                          gCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTa                                                                        >  2:99458/1‑150 (MQ=255)
                                                                                 gCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTg                                                                                                  >  2:1191001/1‑117 (MQ=255)
                                                                                 gCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTg                                                                                                  <  1:1191001/117‑1 (MQ=255)
                                                                                   aTCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAAcatca                                                               <  2:1323575/150‑1 (MQ=255)
                                                                                   aTCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAAcatca                                                               <  1:1045038/150‑1 (MQ=255)
                                                                                     cACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCAtt                                                             >  2:639070/1‑150 (MQ=255)
                                                                                      aCCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTc                                                            <  2:27011/150‑1 (MQ=255)
                                                                                               cccTTTGCGCTCCAGTTTGTGGAGATCATGTTAGTCGCGACGTGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTTCGCCGCTTGGATCTTCGGCATAGGTAAGAAAATCATTCAGCCATTTg                                                   <  2:1479831/150‑1 (MQ=255)
                                                                                               cccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTg                                                   <  1:478148/150‑1 (MQ=255)
                                                                                                ccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGATCATCATTCAGCCCTTTGc                                                  <  2:884488/150‑1 (MQ=255)
                                                                                                ccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGc                                                  <  2:353643/150‑1 (MQ=255)
                                                                                                ccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGc                                                  >  1:400808/1‑150 (MQ=255)
                                                                                                 cATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCa                                                 <  1:639070/150‑1 (MQ=255)
                                                                                                 cATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCa                                                 <  1:214287/150‑1 (MQ=255)
                                                                                                 cATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCa                                                 <  1:1496741/150‑1 (MQ=255)
                                                                                                  aTTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAg                                                                      >  2:272115/1‑128 (MQ=255)
                                                                                                  aTTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAg                                                                      <  1:272115/128‑1 (MQ=255)
                                                                                                  aTTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCAt                                                <  2:994958/150‑1 (MQ=255)
                                                                                                  aTTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCAt                                                <  1:1494080/150‑1 (MQ=255)
                                                                                                     tCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATaaa                                             <  2:311767/150‑1 (MQ=255)
                                                                                                      cgcgCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAg                                            >  1:870773/1‑150 (MQ=255)
                                                                                                        cgcCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCa                                                 >  1:1057073/1‑143 (MQ=255)
                                                                                                        cgcCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCa                                                 <  2:1057073/143‑1 (MQ=255)
                                                                                                        cgcCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGaa                                          >  2:222272/1‑150 (MQ=255)
                                                                                                         gcCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGaaa                                         >  1:1367845/1‑150 (MQ=255)
                                                                                                          ccAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAAt                                        <  1:446078/150‑1 (MQ=255)
                                                                                                             gTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGc                                     >  1:302543/1‑150 (MQ=255)
                                                                                                                 gtggCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTtca                                 >  2:986325/1‑150 (MQ=255)
                                                                                                                   ggCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCGCTTTGCATAAAGAAATGGCTtcatc                               >  1:495204/1‑150 (MQ=255)
                                                                                                                   ggCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTtcatc                               <  2:116215/150‑1 (MQ=255)
                                                                                                                   ggCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTtcatc                               <  2:1487020/150‑1 (MQ=255)
                                                                                                                              gaatttcgACGGGCGACCGGCTTGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGAACTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATAACGCAGCCATg                    >  2:1033408/7‑150 (MQ=255)
                                                                                                                                 tttcgACGGGCGACCGGCTTGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCAc                 >  2:1241876/4‑150 (MQ=255)
                                                                                                                                  ttcgACGGGCGACCGGCTTGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCAcc                >  1:468324/3‑150 (MQ=255)
                                                                                                                                   gcgACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCg               >  2:589770/1‑150 (MQ=255)
                                                                                                                                    cgACGGGCGACCGGCTTGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGt              >  1:1132434/1‑150 (MQ=255)
                                                                                                                                        gggCAAGCCGGTCGCCCTTCGAAATTCCGCCATCACTTCACGTCCGGAGATGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTTCATAAAGAATTGGCTTCATCAAGCAGCCATTCTCCGTTggg          <  2:306006/150‑1 (MQ=255)
                                                                                                                                                cGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACgg  <  2:1022262/150‑1 (MQ=255)
                                                                                                                                                cGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACgg  <  1:295605/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGG  >  NC_000913/2919759‑2920052

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