Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,651,319 G→C 8.8% intergenic (+41/+29) PP_2323 → / ← aroF‑I hypothetical protein/tyr‑sensitive phospho‑2‑dehydro‑3‑deoxyheptonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,651,3190GC 8.8% 274.1 / 16.7 91intergenic (+41/+29)PP_2323/aroF‑Ihypothetical protein/tyr‑sensitive phospho‑2‑dehydro‑3‑deoxyheptonate aldolase
Reads supporting (aligned to +/- strand):  ref base G (41/42);  new base C (2/6);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

GGGCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATACT  >  NC_002947/2651175‑2651460
                                                                                                                                                |                                                                                                                                             
gggCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGccc                                                                                                                                          >  1:606255/1‑150 (MQ=255)
  gCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCgg                                                                                                                                        <  2:123970/150‑1 (MQ=255)
    tctGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCg                                                                                                                                      >  2:620711/1‑150 (MQ=255)
    tctGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCg                                                                                                                                      >  1:695898/1‑150 (MQ=255)
    tctGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCg                                                                                                                                      <  2:1594522/150‑1 (MQ=255)
     ctGGAGTACGCCGAGCAGATGGGGTACACGGTCCTGCTGGCGTGGGCGTAGGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGt                                                                                                                                     <  2:612943/150‑1 (MQ=255)
      tGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGcc                                                                                                                                           <  2:1997696/143‑1 (MQ=255)
      tGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGcc                                                                                                                                           >  1:1997696/1‑143 (MQ=255)
      tGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGtt                                                                                                                                    >  1:1107076/1‑150 (MQ=255)
      tGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGtt                                                                                                                                    <  2:695898/150‑1 (MQ=255)
      tGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGtt                                                                                                                                    <  1:620711/150‑1 (MQ=255)
      tGGAGTACGCCGAGCAGATGGAGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAAGCCTGAGCTGACGACAGGCACCAGACACAACAAGACACGCCGGGCAGTGCCCGGCGtt                                                                                                                                    <  1:53780/150‑1 (MQ=255)
       ggAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGttt                                                                                                                                   <  1:1473656/150‑1 (MQ=255)
       ggAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGttt                                                                                                                                   <  2:1065030/150‑1 (MQ=255)
       ggAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGttt                                                                                                                                   <  2:1064640/150‑1 (MQ=255)
        gAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGtttt                                                                                                                                  <  2:1131607/150‑1 (MQ=255)
              gCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGttttt                                                                                                                                 >  2:1482344/1‑145 (MQ=255)
              gCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGttttt                                                                                                                                 <  1:1482344/145‑1 (MQ=255)
               ccGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTcc                                                                                                                           >  1:2246809/1‑150 (MQ=255)
                 gAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCg                                                                                                                         >  2:1580601/1‑150 (MQ=255)
                    cAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCaa                                                                                                                      <  1:1198053/150‑1 (MQ=255)
                      gATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATg                                                                                                                    >  1:1747641/1‑150 (MQ=255)
                      gATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATg                                                                                                                    >  2:659995/1‑150 (MQ=255)
                         ggggTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTc                                                                                                                 <  2:606255/150‑1 (MQ=255)
                           ggTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAg                                                                                                               >  1:212171/1‑150 (MQ=255)
                           ggTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAg                                                                                                               <  2:1107076/150‑1 (MQ=255)
                             tACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGcc                                                                                                             <  1:1232482/150‑1 (MQ=255)
                                    tGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTa                                                                                                      <  1:2218945/150‑1 (MQ=255)
                                     gATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGTCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTAc                                                                                                     <  2:219357/150‑1 (MQ=255)
                                          cGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGc                                                                                                <  1:1580601/150‑1 (MQ=255)
                                            tCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCgg                                                                                              <  1:659995/150‑1 (MQ=255)
                                               tGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGgcag                                                                                           <  1:1600597/150‑1 (MQ=255)
                                                  tCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGcac                                                                                        >  2:1224265/1‑150 (MQ=255)
                                                       tgtgGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGgca                                                                                            <  2:1166922/141‑1 (MQ=255)
                                                       tgtgGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGgca                                                                                            >  1:1166922/1‑141 (MQ=255)
                                                         tgGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTg                                                                                 >  1:529795/1‑150 (MQ=255)
                                                           gAGCGTTACATGGAGCACCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAg                                                                               >  2:1961412/1‑150 (MQ=255)
                                                              cGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCtt                                                                            <  1:1017776/150‑1 (MQ=255)
                                                               gTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTg                                                                           >  2:1748206/1‑150 (MQ=255)
                                                                ttACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTgg                                                                          >  2:410235/1‑150 (MQ=255)
                                                                          cgcCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTg                                                                                                 >  2:790693/1‑117 (MQ=255)
                                                                          cgcCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTg                                                                                                 <  1:790693/117‑1 (MQ=255)
                                                                           gcCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTccc                                                                                                                          >  1:621308/1‑91 (MQ=255)
                                                                           gcCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTccc                                                                                                                          <  2:621308/91‑1 (MQ=255)
                                                                             cagcagCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTgcg                                                             >  1:667653/1‑150 (MQ=255)
                                                                             cagcagCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTgcg                                                             >  2:1722578/1‑150 (MQ=255)
                                                                               gcagcGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCa                                                           >  1:649401/1‑150 (MQ=255)
                                                                                       aTTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTtc                                                   >  2:1327383/1‑150 (MQ=255)
                                                                                        ttCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGtt                                                                                                                                    >  1:729934/1‑68 (MQ=255)
                                                                                        ttCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGtt                                                                                                                                    <  2:729934/68‑1 (MQ=255)
                                                                                          cGAGCAAGGCCTGAGGTGACGACAGGCACCAGACACAACAAACAACGCCCGGCAGTACCCGGCGTTTTTTATTCCCGCATTGTTAAGCCGCCCTTACGCTGCGGCAGCACCTCCTTGGCCTTGGCATGCATGCTTCGCAGGGTTTTCTCg                                                <  2:310123/150‑1 (MQ=255)
                                                                                          cGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCg                                                <  1:410235/150‑1 (MQ=255)
                                                                                           gAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCgg                                               >  2:2178681/1‑150 (MQ=255)
                                                                                              cAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAAt                                                                                                                     <  2:593328/77‑1 (MQ=255)
                                                                                              cAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAAt                                                                                                                     >  1:593328/1‑77 (MQ=255)
                                                                                               aaGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTGCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACACCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGc                                           <  2:1327045/150‑1 (MQ=255)
                                                                                                  gCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGCGCCTGGCATGCATGCTGCGAAGGGTTTCCTTGGTAGCCGa                                        <  2:1785955/150‑1 (MQ=255)
                                                                                                  gCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTCGCCGa                                        <  2:1747641/150‑1 (MQ=255)
                                                                                                     tGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCa                                     >  1:1824794/1‑150 (MQ=255)
                                                                                                       aGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGt                                   >  2:1240032/1‑150 (MQ=255)
                                                                                                          tGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTg                                                                                 >  1:2117533/1‑101 (MQ=255)
                                                                                                          tGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTg                                                                                 <  2:2117533/101‑1 (MQ=255)
                                                                                                              gacAGGCACCATACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCACTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGGCGATGCa                            >  1:587014/1‑150 (MQ=255)
                                                                                                              gacAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCa                            <  1:1722578/150‑1 (MQ=255)
                                                                                                               acAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAg                           <  2:649401/150‑1 (MQ=255)
                                                                                                               acAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAg                           <  2:529795/150‑1 (MQ=255)
                                                                                                                         ccgatctcaaACAACGCCGGGCAGTGCCCGGCTTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGATCTTGGCATGCATGCTGCGCAGGGTTTTCTCCGTTTCCTACCAGTCGATGCAGGCGTCCGTGa                 <  1:1229463/143‑1 (MQ=255)
                                                                                                                          acacAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGAt                <  1:2322984/150‑1 (MQ=255)
                                                                                                                          acacAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGAt                >  2:1647952/1‑150 (MQ=255)
                                                                                                                           cacaATAAAAAACGCCGGGCAGTGCCCGCCGTTTTTCACCTCCGCAATGTTCAGCCGCACTTACGCTGCGGCAGCACATCCTTCGTCGTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATc               <  2:587014/150‑1 (MQ=255)
                                                                                                                            acaacaAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGc                         <  2:156979/139‑1 (MQ=255)
                                                                                                                            acaacaAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGc                         >  1:156979/1‑139 (MQ=255)
                                                                                                                             caacaaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCgg                                                                                              <  2:465087/69‑1 (MQ=255)
                                                                                                                             caacaaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCgg                                                                                              >  1:465087/1‑69 (MQ=255)
                                                                                                                              aacaaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGAc                                       >  2:1345694/1‑123 (MQ=255)
                                                                                                                              aacaaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGAc                                       <  1:1345694/123‑1 (MQ=255)
                                                                                                                               acaaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGaca           >  1:979494/1‑150 (MQ=255)
                                                                                                                                 aaaaaaCGCCGGGCAGTGCCCGGCGTTTTTCATTACCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCGAAGAGTGCTCGGTAGCCGACTAGTCGATGCAGGCGTCCGTGATAGatcgg         >  2:1229463/1‑146 (MQ=255)
                                                                                                                                  aaaaaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGc        <  1:1647952/150‑1 (MQ=255)
                                                                                                                                    aaaCGCCGGGCACTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAg                           >  2:1408711/1‑129 (MQ=255)
                                                                                                                                    aaaCGCCGGGCACTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAg                           <  1:1408711/129‑1 (MQ=255)
                                                                                                                                       cGCCGGGCAGTGCCCGGCGTTTTTCATTTCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATAc   >  2:1968131/1‑150 (MQ=255)
                                                                                                                                       cGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGCCCAGTCGATGCAGGCGTCCGTGATCGACACGCCATAc   <  1:1961412/150‑1 (MQ=255)
                                                                                                                                        gCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCGGCCGCGCTTACGATGCGGCAGCAGATCCTTGAGCTTGGCATGCATGCTGCGCAAGGTTTTCTCCTTAGCCTACCACTGGATGCCGGAGTAGGTGATGGACAGGACAGACt  >  2:1032559/1‑150 (MQ=255)
                                                                                                                                        gCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATc               >  2:794460/1‑137 (MQ=255)
                                                                                                                                        gCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATc               <  1:794460/137‑1 (MQ=255)
                                                                                                                                        gCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATACt  >  2:469855/1‑150 (MQ=255)
                                                                                                                                          cGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAgg                                                         >  1:918127/1‑93 (MQ=255)
                                                                                                                                          cGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAgg                                                         <  2:918127/93‑1 (MQ=255)
                                                                                                                                             gCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGt                       >  1:1543916/1‑124 (MQ=255)
                                                                                                                                             gCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGt                       <  2:1543916/124‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
GGGCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATACT  >  NC_002947/2651175‑2651460

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