Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,650,954 C→G 7.0% A29P (GCG→CCG) ‡ ygjP ← metal‑dependent hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,650,9540CG 7.0% 610.6 / 27.1 185A29P (GCG→CCG) ‡ygjPmetal‑dependent hydrolase
Reads supporting (aligned to +/- strand):  ref base C (79/93);  new base G (9/4);  total (88/97)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 197 total aligned reads displayed.

CAGGTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAAG  >  NC_002947/5650805‑5651103
                                                                                                                                                     |                                                                                                                                                     
cAGGTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTgcgc                                                                                                                                                       <  1:4123479/150‑1 (MQ=255)
cAGGTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTgcgc                                                                                                                                                       >  1:1904016/1‑150 (MQ=255)
 aGGTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTgcgc                                                                                                                                                       >  2:2711253/1‑149 (MQ=255)
   gTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTg                                                                                                                                                    >  1:81325/1‑150 (MQ=255)
   gTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTg                                                                                                                                                    >  1:3944298/1‑150 (MQ=255)
     ttGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCTCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCa                                                                                                                                                  <  2:3686317/150‑1 (MQ=255)
     ttGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCa                                                                                                                                                  <  2:1205068/150‑1 (MQ=255)
          gcACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGc                                                                                                                                             >  2:82033/1‑150 (MQ=255)
          gcACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGc                                                                                                                                             <  2:780770/150‑1 (MQ=255)
            aCTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCgg                                                                                                                                           >  2:700792/1‑150 (MQ=255)
            aCTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCgg                                                                                                                                           <  2:1790386/150‑1 (MQ=255)
            aCTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCAGGCGCAGGTGCgg                                                                                                                                           <  2:5080781/150‑1 (MQ=255)
             cTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCAGGCGCAGGTGCggc                                                                                                                                          <  1:5005712/150‑2 (MQ=255)
               ctgCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAg                                                                                                                                                 >  2:4130267/3‑141 (MQ=255)
                  cAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTa                                                                                                                                     >  1:418827/1‑150 (MQ=255)
                   aGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAg                                                                                                                                    <  2:1840133/150‑1 (MQ=255)
                     gTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGc                                                                                                                                  >  1:2760977/1‑150 (MQ=255)
                       aCTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCtt                                                                                                                                <  1:4653931/150‑1 (MQ=255)
                           ttGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCaa                                                                                                                            <  1:844453/150‑1 (MQ=255)
                              gcgcAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTa                                                                                                                         >  2:766196/1‑150 (MQ=255)
                               cgcAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTaa                                                                                                                        <  1:4260838/150‑1 (MQ=255)
                               cgcAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTaa                                                                                                                        <  1:766196/150‑1 (MQ=255)
                                 cAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACg                                                                                                                      <  2:5231550/150‑1 (MQ=255)
                                       cTGGGCATGGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTAc                                                                                                                <  1:184776/150‑1 (MQ=255)
                                        tGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCaa                                                                                                                            >  2:1939322/1‑137 (MQ=255)
                                        tGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCaa                                                                                                                            <  1:1939322/137‑1 (MQ=255)
                                        tGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACg                                                                                                               >  1:4045645/1‑150 (MQ=255)
                                         gggCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACgg                                                                                                              <  1:215097/150‑1 (MQ=255)
                                                gCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGgc                                                                                                       >  1:238780/1‑150 (MQ=255)
                                                  cATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGgcgc                                                                                                     >  2:2047265/1‑150 (MQ=255)
                                                   aTACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCt                                                                                                    >  2:3414117/1‑150 (MQ=255)
                                                   aTACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCt                                                                                                    >  2:625626/1‑150 (MQ=255)
                                                      cAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGgcgaagaa                                                                                                 <  2:3311445/150‑6 (MQ=255)
                                                       aGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTCCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTAt                                                                                                >  1:1001088/1‑150 (MQ=255)
                                                             ttGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCg                                                                                          >  2:1080733/1‑150 (MQ=255)
                                                                tCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTaa                                                                                       <  1:2704150/150‑1 (MQ=255)
                                                                 cGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAAc                                                                                      <  2:4045645/150‑1 (MQ=255)
                                                                  gCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACg                                                                                     >  2:4322815/1‑150 (MQ=255)
                                                                  gCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACg                                                                                     >  1:5206598/1‑150 (MQ=255)
                                                                  gCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACg                                                                                     >  1:2868079/1‑150 (MQ=255)
                                                                    tttGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGaa                                                                                   >  1:2951405/1‑150 (MQ=255)
                                                                         aCGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGTGTAACGAATGCCg                                                                              >  1:3859621/1‑150 (MQ=255)
                                                                          cGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGa                                                                             <  2:2350325/150‑1 (MQ=255)
                                                                           gTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCAGGCGCAGGTGCggcaga                                                                                                                                       >  2:5408980/1‑87 (MQ=38)
                                                                           gTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCAGGCGCAGGTGCggcaga                                                                                                                                       <  1:5408980/91‑5 (MQ=38)
                                                                            tCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGAtt                                                                           >  2:3199336/1‑150 (MQ=255)
                                                                             cATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTc                                                                          <  2:2412197/150‑1 (MQ=255)
                                                                             cATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTc                                                                          >  2:2808172/1‑150 (MQ=255)
                                                                             cATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTc                                                                          >  1:1557616/1‑150 (MQ=255)
                                                                              aTGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCt                                                                         <  1:5364992/150‑1 (MQ=255)
                                                                               tGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTa                                                                        >  1:122819/1‑150 (MQ=255)
                                                                               tGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTa                                                                        <  2:1283304/150‑1 (MQ=255)
                                                                               tGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTa                                                                        <  1:3541958/150‑1 (MQ=255)
                                                                               tGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTa                                                                        <  1:3665672/150‑1 (MQ=255)
                                                                                 gCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGc                                                                      <  2:3605070/150‑1 (MQ=255)
                                                                                 gCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGc                                                                      <  2:2206066/150‑1 (MQ=255)
                                                                                 gCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGc                                                                      <  2:4861795/150‑1 (MQ=255)
                                                                                      cAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATc                                                                 <  2:1475283/150‑1 (MQ=255)
                                                                                       aGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCt                                                                                                    <  1:136701/114‑1 (MQ=255)
                                                                                         ggTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCCAGTAACGTAGTACGGTCCTGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACg                                                              >  1:1177881/1‑150 (MQ=255)
                                                                                         ggTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACg                                                              >  2:1322078/1‑150 (MQ=255)
                                                                                                 cgcTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGc                                                      <  1:3414117/150‑1 (MQ=255)
                                                                                                 cgcTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACCTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGc                                                      <  1:3278081/150‑1 (MQ=255)
                                                                                                    tGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCCGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAg                                                   <  1:4144799/150‑1 (MQ=255)
                                                                                                      cAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTTCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGga                                                 <  1:3268185/150‑1 (MQ=255)
                                                                                                         gTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGaa                                              <  2:5390705/150‑1 (MQ=255)
                                                                                                           aTTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGaaaa                                            <  2:3112454/150‑1 (MQ=255)
                                                                                                           aTTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGaaaa                                            >  1:3007911/1‑150 (MQ=255)
                                                                                                               aCCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGc                                        >  2:1124920/1‑150 (MQ=255)
                                                                                                                cccAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGcc                                       >  1:4415908/1‑150 (MQ=255)
                                                                                                                cccAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGcc                                       <  1:3199336/150‑1 (MQ=255)
                                                                                                                cccAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGcc                                       >  2:527568/1‑150 (MQ=255)
                                                                                                                 ccAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCg                                      >  2:1980303/1‑150 (MQ=255)
                                                                                                                 ccAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCg                                      <  2:5105630/150‑1 (MQ=255)
                                                                                                                 ccAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCg                                      <  2:4002488/150‑1 (MQ=255)
                                                                                                                  cAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCgg                                     >  1:2216346/1‑150 (MQ=255)
                                                                                                                   aGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCggg                                    >  1:4494308/1‑150 (MQ=255)
                                                                                                                   aGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCggg                                    >  1:353637/1‑150 (MQ=255)
                                                                                                                     gCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGaa                                  <  2:3944298/150‑1 (MQ=255)
                                                                                                                             tGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAg                          >  2:679638/1‑150 (MQ=255)
                                                                                                                                cGATCATCTGCCGCCCCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGt                       <  2:4197519/150‑1 (MQ=255)
                                                                                                                                cGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCTGGGAAGCCGCTAGCGt                       <  1:773181/150‑1 (MQ=255)
                                                                                                                                 gATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTc                      >  2:3759856/1‑150 (MQ=255)
                                                                                                                                  atcatcTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCa                     <  1:1080733/150‑1 (MQ=255)
                                                                                                                                  atcatcTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCa                     <  1:625626/150‑1 (MQ=255)
                                                                                                                                     atcTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTc                  <  2:4723397/150‑1 (MQ=255)
                                                                                                                                     atcTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTc                  >  1:4362337/1‑150 (MQ=255)
                                                                                                                                     atcTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTc                  >  1:1278044/1‑150 (MQ=255)
                                                                                                                                      tcTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTcg                 >  1:1398854/1‑150 (MQ=255)
                                                                                                                                        ccccGCACCTGCAGGCGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTcgcg               >  1:2856131/3‑150 (MQ=255)
                                                                                                                                        ccccGCACCTGCAGGCGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTcgcg               >  1:1875226/3‑150 (MQ=255)
                                                                                                                                        ccccGCACCTGCAGGCGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTcgcg               >  1:1874913/3‑150 (MQ=255)
                                                                                                                                         cccGCACCTGCAGGCGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGga                                                 >  1:3755759/2‑115 (MQ=255)
                                                                                                                                         cccGCACCTGCAGGCGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCgg              >  1:3168437/2‑150 (MQ=255)
                                                                                                                                         cccGCACCTGCAGGCGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCgg              >  2:3832078/2‑150 (MQ=255)
                                                                                                                                          ccGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGt             >  1:3066175/1‑150 (MQ=255)
                                                                                                                                              aCCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATc         <  2:1424576/150‑1 (MQ=255)
                                                                                                                                                 tGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGa      <  1:2788101/150‑1 (MQ=255)
                                                                                                                                                    gcCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATaa   <  2:1878048/150‑1 (MQ=255)
                                                                                                                                                     ccTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAAg  >  2:2874576/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGGTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAAG  >  NC_002947/5650805‑5651103

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