Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,138,457 C→A 6.0% R173S (CGC→AGC)  yadG → ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,138,4570CA 6.0% 321.3 / 6.7 100R173S (CGC→AGC) yadGABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (43/51);  new base A (3/3);  total (46/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

Only 100 of 101 total aligned reads displayed.

GCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAGCATGCGCCAG  >  NC_002947/2138314‑2138606
                                                                                                                                               |                                                                                                                                                     
gCCTGTGGTACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGa                                                                                                                                                 <  1:1225078/150‑1 (MQ=255)
 ccTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGAt                                                                                                                                                >  2:1938422/1‑150 (MQ=255)
 ccTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGAt                                                                                                                                                <  2:771516/150‑1 (MQ=255)
 ccTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGAt                                                                                                                                                <  1:252031/150‑1 (MQ=255)
 ccTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGAt                                                                                                                                                <  1:507386/150‑1 (MQ=255)
  cTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGAtg                                                                                                                                               <  2:244997/150‑1 (MQ=255)
  cTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGAtg                                                                                                                                               >  1:883927/1‑150 (MQ=255)
   tgtgGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACAGAGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGAtgt                                                                                                                                              >  1:1779881/1‑150 (MQ=255)
    gtgGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGAtgtg                                                                                                                                             <  1:974887/150‑1 (MQ=255)
     tgGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCGGCGCAGCTCGAtatcc                                                                                                                                            <  1:1284520/150‑5 (MQ=255)
           aaaCGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTc                                                                                                                                      <  1:663810/150‑1 (MQ=255)
           aaaCGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTc                                                                                                                                      >  2:83296/1‑150 (MQ=255)
           aaaCGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACAGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCg                                                                                                                                                  <  2:2009402/138‑1 (MQ=255)
           aaaCGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACAGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCg                                                                                                                                                  >  1:2009402/1‑138 (MQ=255)
                 gATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACt                                                                                                                                <  1:1938422/150‑1 (MQ=255)
                 gATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCCGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACt                                                                                                                                >  2:1368904/1‑150 (MQ=255)
                       cAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTg                                                                                                                          >  2:342064/1‑150 (MQ=255)
                        aGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCTCCGCTCGATGTGGAGCTTCCTCACTGAGCTGa                                                                                                                         <  1:395181/150‑1 (MQ=255)
                          tCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCATGGACGAACACCCCGCCGGTGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAAc                                                                                                                       >  1:2028817/1‑150 (MQ=255)
                          tCGCGTTACTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCAAGAACGCCTCCGTCTGATACCGTTCGAACCCACCGCCGGGGTGGCTATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAAc                                                                                                                       <  1:755925/150‑1 (MQ=255)
                            gcgTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCa                                                                                                                     >  2:735947/1‑150 (MQ=255)
                            gcgTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCa                                                                                                                     >  1:2315305/1‑150 (MQ=255)
                                   ttGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGc                                                                                                              >  2:859418/1‑150 (MQ=255)
                                   ttGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGc                                                                                                              >  2:1986243/1‑150 (MQ=255)
                                              cATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACcatcat                                                                                                   <  1:1148269/150‑1 (MQ=255)
                                               aTGAAGCGGCGCCTGATGATTGCCCGCGCATTGGTCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGc                                                                                                                            >  2:114558/1‑124 (MQ=255)
                                               aTGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGc                                                                                                                            <  1:114558/124‑1 (MQ=255)
                                               aTGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCatcatc                                                                                                  >  1:1629334/1‑150 (MQ=255)
                                                tGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATcc                                                                                                 <  1:83296/150‑1 (MQ=255)
                                                tGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATcc                                                                                                 <  1:308011/150‑1 (MQ=255)
                                                 gAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCt                                                                                                >  2:1260344/1‑150 (MQ=255)
                                                  aaGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTc                                                                                               <  2:1629334/150‑1 (MQ=255)
                                                   aGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCt                                                                                                                           >  2:1505851/1‑121 (MQ=255)
                                                   aGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCt                                                                                                                           <  1:1505851/121‑1 (MQ=255)
                                                        cgcCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCaccacc                                                                                         >  1:1638972/1‑150 (MQ=255)
                                                         gcCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCt                                                                                                                                    <  2:2161547/106‑1 (MQ=255)
                                                         gcCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCt                                                                                                                                    >  1:2161547/1‑106 (MQ=255)
                                                         gcCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCAccaccc                                                                                        <  2:883927/150‑1 (MQ=255)
                                                            tgatgaTTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGc                                                                                                              <  2:113618/125‑1 (MQ=255)
                                                            tgatgaTTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGc                                                                                                              >  1:113618/1‑125 (MQ=255)
                                                             gatgatTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCAGTGAACTGAACCAGAAGGGCATCACCATCATCCTCACCACACACTa                                                                                    >  2:1875623/1‑150 (MQ=255)
                                                              atgatTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACt                                                                                                                                <  2:2237116/105‑1 (MQ=255)
                                                              atgatTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACt                                                                                                                                >  1:2237116/1‑105 (MQ=255)
                                                               tgatTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATc                                                                                  >  1:1066793/1‑150 (MQ=255)
                                                                 atTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTg                                                                                <  2:243259/150‑1 (MQ=255)
                                                                    gCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGaa                                                                             <  1:342064/150‑1 (MQ=255)
                                                                      ccGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTAACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAaga                                                                           >  1:1985749/1‑150 (MQ=255)
                                                                        gcgcATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAgg                                                                         <  1:117498/150‑1 (MQ=255)
                                                                          gcATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCt                                                                       <  1:162366/150‑1 (MQ=255)
                                                                          gcATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCt                                                                       <  1:920245/150‑1 (MQ=255)
                                                                          gcATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCt                                                                       <  1:1260344/150‑1 (MQ=255)
                                                                                aTCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGagcag                                                                 <  1:859418/150‑1 (MQ=255)
                                                                                 tCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGc                                                                                                                            <  2:750437/90‑1 (MQ=255)
                                                                                 tCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGc                                                                                                                            >  1:750437/1‑90 (MQ=255)
                                                                                  ccACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCt                                                               <  2:2315305/150‑1 (MQ=255)
                                                                                  ccACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCt                                                               <  2:569833/150‑1 (MQ=255)
                                                                                    aCGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCtgt                                                             >  2:2051373/1‑150 (MQ=255)
                                                                                      gAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGCGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGc                                                           <  1:1268803/150‑1 (MQ=255)
                                                                                          cgcgCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAgg                                                                                                                <  2:2320830/93‑1 (MQ=255)
                                                                                          cgcgCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAgg                                                                                                                >  1:2320830/1‑93 (MQ=255)
                                                                                           gcgcCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCCCGATGTGGAGCTTCCTCACTGCGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGGCGTaa                                                      >  1:429839/1‑150 (MQ=255)
                                                                                            cgcCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCAccaccca                                                                                       <  1:1816509/116‑1 (MQ=255)
                                                                                            cgcCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCAccaccca                                                                                       >  2:1816509/1‑116 (MQ=255)
                                                                                             gcCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCAt                                                                                                            >  2:1827586/1‑94 (MQ=255)
                                                                                             gcCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCAt                                                                                                            <  1:1827586/94‑1 (MQ=255)
                                                                                             gcCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACa                                                    <  2:1488567/150‑1 (MQ=255)
                                                                                              ccTGCTGATCCTGGACGAACCCACCGCCTGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACAt                                                   <  2:1713179/150‑1 (MQ=255)
                                                                                               ctgctgATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATc                                                  >  2:1685111/1‑150 (MQ=255)
                                                                                                   tgATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCatcatc                                                                                                  >  1:1265276/1‑98 (MQ=255)
                                                                                                   tgATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCatcatc                                                                                                  <  2:1265276/98‑1 (MQ=255)
                                                                                                   tgATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGca                                              >  2:1712297/1‑150 (MQ=255)
                                                                                                       ccTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGcatcat                                          <  1:2051373/150‑1 (MQ=255)
                                                                                                       ccTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGcatcat                                          >  1:533468/1‑150 (MQ=255)
                                                                                                        cTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAaga                                                                           <  1:1877264/116‑1 (MQ=255)
                                                                                                        cTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAaga                                                                           >  2:1877264/1‑116 (MQ=255)
                                                                                                        cTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGACCAAGCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCatcatc                                         <  1:735947/150‑1 (MQ=255)
                                                                                                         tGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCtgt                                                             >  2:252663/1‑129 (MQ=255)
                                                                                                         tGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCtgt                                                             <  1:252663/129‑1 (MQ=255)
                                                                                                         tGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCg                                        >  2:2207532/1‑150 (MQ=255)
                                                                                                          ggACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAAACATAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGa                                       <  1:381140/150‑1 (MQ=255)
                                                                                                           gACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACca                                                                                                       <  2:1587881/85‑1 (MQ=255)
                                                                                                           gACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACca                                                                                                       >  1:1587881/1‑85 (MQ=255)
                                                                                                           gACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGAc                                      >  2:753085/1‑150 (MQ=255)
                                                                                                                  ccACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCa                               >  2:1323736/1‑150 (MQ=255)
                                                                                                                       gccgGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCAtc                          <  2:1066793/150‑1 (MQ=255)
                                                                                                                          ggggTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTCCCTCATTGAGCTGAACCTGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGGATCATCGACCTTGGCAGCAtcgtc                       <  2:1779881/150‑1 (MQ=255)
                                                                                                                            ggTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCga                     <  2:1638972/150‑1 (MQ=255)
                                                                                                                             gTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGGAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCgag                    >  2:379698/1‑150 (MQ=255)
                                                                                                                              tGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCgaga                   >  2:1490119/1‑150 (MQ=255)
                                                                                                                               ccacatCGAGCAGCGCAGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTg                                                                                <  2:1282528/84‑1 (MQ=255)
                                                                                                                               ccacatCGAGCAGCGCAGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTg                                                                                >  1:1282528/5‑88 (MQ=255)
                                                                                                                                gATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAac                 <  1:2207532/150‑1 (MQ=255)
                                                                                                                                cacatCGAGCGGCGCAGCTCGATATGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGcc                                                          >  2:169366/4‑109 (MQ=255)
                                                                                                                                cacatCGAGCGGCGCAGCTCGATATGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGcc                                                          <  1:169366/106‑1 (MQ=255)
                                                                                                                                  tatCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCTACCATGGCACCATCGTCGAGAacac               >  2:1821621/1‑150 (MQ=255)
                                                                                                                                   atCGAGCGGCGCAGCTCGATGTGGAGCTTCCTCACTGAGCTGGACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACAcc              >  2:866930/1‑150 (MQ=255)
                                                                                                                                    tCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCCCTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCa             <  2:429839/150‑1 (MQ=255)
                                                                                                                                        gCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAGCAt         >  2:1399788/1‑150 (MQ=255)
                                                                                                                                          tGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAGCATgc       <  2:533468/150‑1 (MQ=255)
                                                                                                                                               cgcTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAGCATGCGCCAg  >  1:1068915/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAGCATGCGCCAG  >  NC_002947/2138314‑2138606

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