Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,970,421 C→A 5.2% G43G (GGG→GGT PP_5723 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,970,4210CA 5.2% 327.4 / 6.4 96G43G (GGG→GGTPP_5723hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/53);  new base A (3/2);  total (41/55)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 103 total aligned reads displayed.

CTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAG  >  NC_002947/5970272‑5970569
                                                                                                                                                     |                                                                                                                                                     
cTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGcc                                                                                                                                                       <  1:870605/150‑1 (MQ=255)
 tCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGccc                                                                                                                                                      >  2:501291/1‑150 (MQ=255)
  cGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGcccc                                                                                                                                                     <  1:4056275/150‑1 (MQ=255)
     gATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCa                                                                                                                                                  <  1:803600/150‑1 (MQ=255)
      aTTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCaa                                                                                                                                                 <  1:3411442/150‑1 (MQ=255)
      aTTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCaa                                                                                                                                                 <  1:3411418/150‑1 (MQ=255)
      aTTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCaa                                                                                                                                                 <  1:1935296/150‑1 (MQ=255)
        tGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGa                                                                                                                                               <  2:3793665/150‑1 (MQ=255)
        tGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCGTGAAAACTGCACGCGATCGTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGa                                                                                                                                               <  2:1686204/150‑1 (MQ=255)
         ggCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGAt                                                                                                                                              <  2:1178182/150‑1 (MQ=255)
         ggCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGAt                                                                                                                                              <  2:1178224/150‑1 (MQ=255)
           ccAGTCGCCCGCACGGGAAACCTGCACACTTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGAttt                                                                                                                                            <  1:2591045/150‑1 (MQ=255)
           ccAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGAttt                                                                                                                                            <  2:4286466/150‑1 (MQ=255)
           ccAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGGCGGGGTTTGGGTGGGATGGGCCCCCCCAGGttt                                                                                                                                            >  1:4105499/1‑150 (MQ=255)
             aGTCGCCAGCCCGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCACGCTGCGCGAAAAATGGGCAACAGCTTGAAAACGGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTtgt                                                                                                                                          <  2:2346173/150‑1 (MQ=255)
             aGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTtgt                                                                                                                                          <  1:220592/150‑1 (MQ=255)
             aGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTtgt                                                                                                                                          <  2:3698848/150‑1 (MQ=255)
             aGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTtgt                                                                                                                                          <  2:3607952/150‑1 (MQ=255)
             aGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTtgt                                                                                                                                          <  2:1068511/150‑1 (MQ=255)
             aGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGACTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTtgt                                                                                                                                          <  1:331987/150‑1 (MQ=255)
               tCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGc                                                                                                                                        <  1:1930780/150‑1 (MQ=255)
                 gCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCt                                                                                                                                      <  1:2664550/150‑1 (MQ=255)
                  ccAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCtt                                                                                                                                     <  1:960202/150‑1 (MQ=255)
                  ccAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCtt                                                                                                                                     <  2:1402872/150‑1 (MQ=255)
                  ccAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCtt                                                                                                                                     <  2:1886651/150‑1 (MQ=255)
                   cAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTg                                                                                                                                    <  1:45622/150‑1 (MQ=255)
                          ggAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGa                                                                                                                             <  1:3748552/150‑1 (MQ=255)
                            aaaCATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGata                                                                                                                          <  1:4071094/150‑1 (MQ=255)
                              aCATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATAtt                                                                                                                        <  1:206149/150‑1 (MQ=255)
                                aTGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCCGCGCGCTTGGGGTGCGAAGCCCCCCCCCACCTTTGTGCCGTGTTGCTGTGATTgg                                                                                                                      >  2:2941872/1‑149 (MQ=255)
                                 tGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGcc                                                                                                                     <  1:822788/150‑1 (MQ=255)
                                   cacaATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCgg                                                                                                                   <  2:3146637/150‑1 (MQ=255)
                                   cacaATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCgg                                                                                                                   <  2:1209592/150‑1 (MQ=255)
                                   cacaATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCgg                                                                                                                   <  1:1210079/150‑1 (MQ=255)
                                   cacaATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCgg                                                                                                                   <  2:3851859/150‑1 (MQ=255)
                                   cacaATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCgg                                                                                                                   <  2:3851901/150‑1 (MQ=255)
                                    acaATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCggg                                                                                                                  <  1:536592/150‑1 (MQ=255)
                                          tCGGTGCTCAGCTGCTGGGCATTATTCGCGTAAAATGTGCTAAAGCTTGAAATCTGCTGGATATCTTTGTTGGAGGGGCCTTTCGTCGCGATTGGGTTTCGAAGCGGCCCCCTTGATTTGTTGCTTGATGCTGATATTGCCGTGTCCGCt                                                                                                            <  1:1210856/150‑1 (MQ=255)
                                               gcgcAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCa                                                                                                       <  1:1862274/150‑1 (MQ=255)
                                                cgcAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAg                                                                                                      <  1:501291/150‑1 (MQ=255)
                                                 gcagcaGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGc                                                                                                     <  2:2069609/150‑1 (MQ=255)
                                                   agcagcGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGGTTGGGGTGGGAAGCGCCCCCCAAGATTTGTGGGGTGTTTGTGGTATTTGCGGGGGCGCGGTGCAGgcc                                                                                                   >  2:3709076/1‑150 (MQ=255)
                                                       gcGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATc                                                                                               <  1:1244889/150‑1 (MQ=255)
                                                       gcGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATc                                                                                               <  1:2514035/150‑1 (MQ=255)
                                                       gcGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATc                                                                                               <  2:1312872/150‑1 (MQ=255)
                                                         ggggCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATcgc                                                                                             <  2:559063/150‑1 (MQ=255)
                                                            gCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGAc                                                                                          <  1:1965385/150‑1 (MQ=255)
                                                              aaGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGc                                                                                        <  2:1593232/150‑1 (MQ=255)
                                                                 cTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAg                                                                                     <  2:1517274/150‑1 (MQ=255)
                                                                     gcgAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCGTGCGGCGCGAAGGGGGTTCCAAGCGACCCCCCAGAATTGTGGCCTGTTGGTGATAATGGCCGGGCCCCCGGGCCGCCCCAACCCGCCGCAAGGCCg                                                                                 >  2:1898748/1‑150 (MQ=255)
                                                                       gAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGGGGGGCGCATGGGGTGGGCATCCCCCCCCCACCATTTTGCCCGGGTGCGGGTAATTGACGGGCCGCAGGGCAGCCCCAACGCCACGCAACGccccc                                                                               >  2:2865736/1‑147 (MQ=255)
                                                                       gAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGGGTGGGGGGGGAAGCGCCCCCCAACAATTGTTGCTTGTAGCTGATAATTGCGGGGCCCCGGGTGAGCCCAATAGCCACGCCAAGCCCca                                                                               >  2:1818164/1‑150 (MQ=255)
                                                                                            aaaCTGCACGCGATCTTTGTGGGAGCGGCCTTGGGGCGGCCTTGGGGTGGGGAACCCCCCCCAACAATTTTGGCGTGATGGTGGTGTTTGCGGGGCCGCGGTGCCGCCCAAACGCCAAGCCAGGCCCCACACAAAAGCCCAggtggaaga                                                          >  2:2707450/1‑145 (MQ=255)
                                                                                                   ccgctttCTTTGTGGGCGGGGCTTTGGGTGGGTTTGGGGGTGCGACCCCCCCCCCAATTGTGGGGCCTGTATGCTTATTTGGCGGGGGCCGGTGTGCCCCCCATCCGCCCCCCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGcc                                                   <  2:3595484/144‑1 (MQ=255)
                                                                                                    cgcgCTGTTGGTGGGCGCGGCCTGGGGGCGCGTTGGGGGGGCGACGCCCCCCCCAAGTTTGGTGCCGTGATGCTGATATCGCGGGGGCCGCGGTGCCCCCCACCCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCa                                                  <  2:4105499/150‑1 (MQ=255)
                                                                                                     gcgATCTTTGTGGGAGCGGCCTTGCGGCGCGTTTGGGGTGCGAACCGCCCCCCCAGATTTGTTGCGGGGTGCTGGTAATGCCGGGGCCCCAGGGGAGCCCAAAGGCGACGCAAGGCCCCCCACCAAGGCCAATGGGGACAAAGGCGcccg                                                 >  2:313241/1‑148 (MQ=255)
                                                                                                           ttggtgGGCGCGGCCTGGCGTGGCGTTGGGGGTGCGATGCCCCCCCCATTTTTGGTGCTGTGTGGCTTTTTTGGCGGGGGCGGCTGGGCCGCCCAATCGCGCCGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTtgt                                           <  2:3333926/147‑1 (MQ=255)
                                                                                                             tgtgGGTGCGGCCTTGCTTCGCTTTTGGGGCGCGACGCCCCCCCCATGTTTTGTGCGTTGTGGCTGTTTTGGCCGGGGCCGGTGTGCCGCCCACTCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTgtgc                                         <  2:4467382/150‑1 (MQ=255)
                                                                                                                   aGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGcc                                                                                                    <  1:2453104/85‑1 (MQ=255)
                                                                                                                        ccTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTg                              >  1:1779021/1‑150 (MQ=255)
                                                                                                                         cTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTgg                             >  1:287916/1‑150 (MQ=255)
                                                                                                                          ttGCGTCGCTTTTGGGGTGCGACGCGCCCCAAAAGATTTGTGCATTTTTGGTTTTTTTGCGGGGGCGGGTGTGCAGCACAATCGCGCCGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTggc                            <  2:2575931/150‑1 (MQ=255)
                                                                                                                          ttGCGTCGCAATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTggc                            >  1:148076/1‑150 (MQ=255)
                                                                                                                             cGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTggcggc                         >  2:3053194/1‑150 (MQ=255)
                                                                                                                              gTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGcc                        >  2:4426096/1‑150 (MQ=255)
                                                                                                                              gTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGcc                        >  2:4425912/1‑150 (MQ=255)
                                                                                                                              gTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGcc                        >  2:631565/1‑150 (MQ=255)
                                                                                                                              gTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGcc                        >  1:929827/1‑150 (MQ=255)
                                                                                                                              acaaatcttGGGGGGCGCTTCGCACCCCAATATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGcc                        >  1:2499420/8‑150 (MQ=255)
                                                                                                                               tCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCCGGTTTGTGCGGGCGACACTGGCGGCCa                       >  1:1575453/1‑150 (MQ=255)
                                                                                                                               tCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCa                       >  1:3466522/1‑150 (MQ=255)
                                                                                                                               tCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCa                       >  1:625315/1‑150 (MQ=255)
                                                                                                                                  atcttGGGGGGCGCTTCGCACCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTggc                            >  2:679836/4‑142 (MQ=255)
                                                                                                                                  atcttGGGGGGCGCTTCGCACCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTggc                            >  2:679907/4‑142 (MQ=255)
                                                                                                                                  atcttGGGGGGCGCTTCGCACCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTggc                            <  1:679836/139‑1 (MQ=255)
                                                                                                                                  atcttGGGGGGCGCTTCGCACCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTggc                            <  1:679907/139‑1 (MQ=255)
                                                                                                                                   gATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTaca                                 >  1:1148636/1‑136 (MQ=255)
                                                                                                                                   gATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTaca                                 <  2:1148636/136‑1 (MQ=255)
                                                                                                                                   gATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTtc                   >  2:613138/1‑150 (MQ=255)
                                                                                                                                        gggTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGt              >  2:3551430/1‑150 (MQ=255)
                                                                                                                                         ggTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGtt             >  1:2794808/1‑150 (MQ=255)
                                                                                                                                         ggTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGtt             >  1:2613434/1‑150 (MQ=255)
                                                                                                                                          gTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGttt            >  2:1592061/1‑150 (MQ=255)
                                                                                                                                          gTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGttt            >  1:3609113/1‑150 (MQ=255)
                                                                                                                                          gTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGttt            >  2:1435162/1‑150 (MQ=255)
                                                                                                                                          gTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGttt            >  2:1470703/1‑150 (MQ=255)
                                                                                                                                           tGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAgg                                                                       <  1:452558/90‑1 (MQ=255)
                                                                                                                                           tGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAgg                                                                       >  2:452558/1‑90 (MQ=255)
                                                                                                                                           tGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTc           >  2:3989412/1‑150 (MQ=255)
                                                                                                                                            gCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTACCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCa          >  1:3096666/1‑150 (MQ=255)
                                                                                                                                             cccaGCGCCCCCCATGTTTGGTGCTTTGCTGCTGATTTGGCGGGGGCCGCTGGGCCCCCCACCCGCGCCGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAg         <  2:2191581/147‑1 (MQ=255)
                                                                                                                                              gAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGt        >  1:2986984/1‑150 (MQ=255)
                                                                                                                                              gAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGt        >  2:2849274/1‑150 (MQ=255)
                                                                                                                                              gAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGt        >  1:86954/1‑150 (MQ=255)
                                                                                                                                              gAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTGTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGt        >  2:861569/1‑150 (MQ=255)
                                                                                                                                               aaGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTg       >  1:2348277/1‑150 (MQ=255)
                                                                                                                                                 gccccccccAAGTTTGGTGCCTTGCTGCTGTTTTTGCGGGGGCCGGTGTGCCGCCCACCCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGa     <  2:1832802/147‑1 (MQ=255)
                                                                                                                                                   gcCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACa   >  1:1381956/1‑150 (MQ=255)
                                                                                                                                                    ccccccAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAg  >  2:2983198/1‑150 (MQ=255)
                                                                                                                                                    ccccccAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAg  >  2:263884/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCATTTCTCCGTTTCAGTGGACAG  >  NC_002947/5970272‑5970569

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