Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,650,953 G→A 7.0% A29V (GCG→GTG) ‡ ygjP ← metal‑dependent hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,650,9530GA 7.0% 564.9 / 21.7 186A29V (GCG→GTG) ‡ygjPmetal‑dependent hydrolase
Reads supporting (aligned to +/- strand):  ref base G (80/93);  new base A (9/4);  total (89/97)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 196 total aligned reads displayed.

CAGGTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAA  >  NC_002947/5650805‑5651102
                                                                                                                                                    |                                                                                                                                                     
cAGGTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTgcgc                                                                                                                                                      <  1:4123479/150‑1 (MQ=255)
cAGGTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTgcgc                                                                                                                                                      >  1:1904016/1‑150 (MQ=255)
cAGGTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTgcgc                                                                                                                                                      <  2:1960963/150‑1 (MQ=255)
 aGGTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTgcgc                                                                                                                                                      >  2:2711253/1‑149 (MQ=255)
   gTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTg                                                                                                                                                   >  1:3944298/1‑150 (MQ=255)
     ttGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCTCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCa                                                                                                                                                 <  2:3686317/150‑1 (MQ=255)
          gcACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGc                                                                                                                                            <  2:780770/150‑1 (MQ=255)
          gcACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGc                                                                                                                                            >  2:82033/1‑150 (MQ=255)
            aCTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCgg                                                                                                                                          >  2:700792/1‑150 (MQ=255)
            aCTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCAGGCGCAGGTGCgg                                                                                                                                          <  2:5080781/150‑1 (MQ=255)
             cTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCAGGCGCAGGTGCggc                                                                                                                                         <  1:5005712/150‑2 (MQ=255)
              ttCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCgggg                                                                                                                                        >  2:181958/1‑150 (MQ=255)
               ctgCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAg                                                                                                                                                >  2:4130267/3‑141 (MQ=255)
                  cAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTa                                                                                                                                    >  1:418827/1‑150 (MQ=255)
                     gTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGc                                                                                                                                 >  1:2760977/1‑150 (MQ=255)
                       aCTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCtt                                                                                                                               <  1:4653931/150‑1 (MQ=255)
                        cTGTTGGCGCAGATCCTGGGCATGGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTg                                                                                                                              >  2:184776/1‑150 (MQ=255)
                              gcgcAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTa                                                                                                                        >  2:766196/1‑150 (MQ=255)
                               cgcAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTaa                                                                                                                       <  1:4260838/150‑1 (MQ=255)
                               cgcAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTaa                                                                                                                       <  1:766196/150‑1 (MQ=255)
                               cgcAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTaa                                                                                                                       <  2:1207821/150‑1 (MQ=255)
                                 cAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACg                                                                                                                     <  2:5231550/150‑1 (MQ=255)
                                      ccTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACgtagta                                                                                                                <  2:3285762/150‑1 (MQ=255)
                                      ccTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACgtagta                                                                                                                <  1:3751067/150‑1 (MQ=255)
                                       cTGGGCATGGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTAc                                                                                                               <  1:184776/150‑1 (MQ=255)
                                        tGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCaa                                                                                                                           <  1:1939322/137‑1 (MQ=255)
                                        tGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACg                                                                                                              >  1:4045645/1‑150 (MQ=255)
                                         gggCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACgg                                                                                                             <  1:215097/150‑1 (MQ=255)
                                           gCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTc                                                                                                           >  2:1440807/1‑150 (MQ=255)
                                               aGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATgg                                                                                                       <  1:3221307/150‑1 (MQ=255)
                                               aGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATgg                                                                                                       >  2:2102446/1‑150 (MQ=255)
                                                gCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGgc                                                                                                      >  1:238780/1‑150 (MQ=255)
                                                   aTACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCt                                                                                                   >  2:625626/1‑150 (MQ=255)
                                                     aCAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGt                                                                                                 <  2:1483089/150‑1 (MQ=255)
                                                       aGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTAt                                                                                               >  2:3477001/1‑150 (MQ=255)
                                                       aGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTCCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTAt                                                                                               >  1:1001088/1‑150 (MQ=255)
                                                             ttGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCCCTGTATTCGGCg                                                                                         >  2:2016052/1‑150 (MQ=255)
                                                                tCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTaa                                                                                      <  1:2704150/150‑1 (MQ=255)
                                                                 cGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAAc                                                                                     <  2:4045645/150‑1 (MQ=255)
                                                                  gCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACg                                                                                    >  1:2868079/1‑150 (MQ=255)
                                                                  gCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACg                                                                                    >  2:4322815/1‑150 (MQ=255)
                                                                  gCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACg                                                                                    >  1:5206598/1‑150 (MQ=255)
                                                                    tttGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGaa                                                                                  >  1:2951405/1‑150 (MQ=255)
                                                                     ttGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAAt                                                                                 <  2:1289870/150‑1 (MQ=255)
                                                                          cGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGa                                                                            <  1:3402666/150‑1 (MQ=255)
                                                                          cGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGa                                                                            <  2:2350325/150‑1 (MQ=255)
                                                                           gTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCAGGCGCAGGTGCggcaga                                                                                                                                      >  2:5408980/1‑87 (MQ=38)
                                                                           gTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCAGGCGCAGGTGCggcaga                                                                                                                                      <  1:5408980/91‑5 (MQ=38)
                                                                             cATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTc                                                                         >  1:1557616/1‑150 (MQ=255)
                                                                             cATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTc                                                                         <  2:2412197/150‑1 (MQ=255)
                                                                             cATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTc                                                                         >  2:2808172/1‑150 (MQ=255)
                                                                              aTGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCt                                                                        <  1:5364992/150‑1 (MQ=255)
                                                                               tGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTa                                                                       >  1:122819/1‑150 (MQ=255)
                                                                               tGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTa                                                                       <  1:3536063/150‑1 (MQ=255)
                                                                                 gCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGc                                                                     <  2:4861795/150‑1 (MQ=255)
                                                                                 gCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGc                                                                     <  2:3605070/150‑1 (MQ=255)
                                                                                       aGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCt                                                                                                   <  1:136701/114‑1 (MQ=255)
                                                                                       aGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCt                                                                                                   >  2:136701/1‑114 (MQ=255)
                                                                                         ggTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCCAGTAACGTAGTACGGTCCTGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACg                                                             >  1:1177881/1‑150 (MQ=255)
                                                                                           tAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGc                                                           <  1:3643304/150‑1 (MQ=255)
                                                                                                 cgcTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACCTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGc                                                     <  1:3278073/150‑1 (MQ=255)
                                                                                                    tGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCCGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAg                                                  <  1:4144799/150‑1 (MQ=255)
                                                                                                      cAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGga                                                <  2:1111768/150‑1 (MQ=255)
                                                                                                         gTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGaa                                             <  2:5390705/150‑1 (MQ=255)
                                                                                                           aTTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGaaaa                                           <  1:78281/150‑1 (MQ=255)
                                                                                                           aTTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGaaaa                                           >  1:3007911/1‑150 (MQ=255)
                                                                                                                cccAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGcc                                      >  1:4415908/1‑150 (MQ=255)
                                                                                                                cccAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGcc                                      <  1:3828839/150‑1 (MQ=255)
                                                                                                                cccAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGcc                                      <  1:3199336/150‑1 (MQ=255)
                                                                                                                cccAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGcc                                      >  2:527568/1‑150 (MQ=255)
                                                                                                                 ccAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCg                                     <  2:1007318/150‑1 (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:353654/1‑150 (MQ=255)
                                                                                                                    ggCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGa                                  <  2:3129808/150‑1 (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)
                                                                                                                                         cccGCACCTGCAGGCGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCgg             >  2:3327229/2‑150 (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:1904016/150‑1 (MQ=255)
                                                                                                                                               ccTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCg       <  1:82033/150‑1 (MQ=255)
                                                                                                                                                 tGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGa     <  1:2788101/150‑1 (MQ=255)
                                                                                                                                                  gcgcCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCACTAGCGTCAGTCGCGGTCATCGGAt    >  2:1197187/1‑150 (MQ=255)
                                                                                                                                                    gcCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATaa  <  2:3053192/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGGTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAA  >  NC_002947/5650805‑5651102

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: