Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,941,651 C→T 6.4% V142M (GTG→ATG)  PP_3477 ← type II secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,941,6510CT 6.4% 182.8 / 2.5 63V142M (GTG→ATG) PP_3477type II secretion protein
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base T (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

GCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACTTGGCGCCCTTGGCACTGC  >  NC_002947/3941505‑3941800
                                                                                                                                                  |                                                                                                                                                     
gcgGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTc                                                                                                                                                    <  1:483590/150‑1 (MQ=255)
gcgGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTc                                                                                                                                                    <  1:811175/150‑1 (MQ=255)
   gCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGaa                                                                                                                                                 >  2:797651/1‑150 (MQ=255)
         cGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTc                                                                                                                                           >  1:1084838/1‑150 (MQ=255)
                       atGGCGATCATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTcgcgcg                                                                                                                             >  2:585174/2‑150 (MQ=255)
                        tGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCgg                                                                                                                            <  2:1143795/150‑1 (MQ=255)
                            gatgatCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTg                                                                                                                        <  2:1084838/150‑1 (MQ=255)
                                atCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGt                                                                                                                       >  1:825637/1‑147 (MQ=255)
                                atCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGt                                                                                                                       <  2:825637/147‑1 (MQ=255)
                                atCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCgg                                                                                                                    <  2:22138/150‑1 (MQ=255)
                                 tCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCggg                                                                                                                   >  2:485046/1‑150 (MQ=255)
                                 tCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACAAGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCggg                                                                                                                   <  1:742507/150‑1 (MQ=255)
                                 tCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTACGATTCAGCTCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGGCTTTCACACGTGCTCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCggg                                                                                                                   >  1:874207/1‑150 (MQ=255)
                                       tCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGgctgct                                                                                                             >  1:664183/1‑150 (MQ=255)
                                        caccacCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGctgctg                                                                                                            <  2:27851/150‑1 (MQ=255)
                                         accaccAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCgg                                                                                                                    <  1:873438/141‑1 (MQ=255)
                                         accaccAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCgg                                                                                                                    >  2:873438/1‑141 (MQ=255)
                                          ccaccaGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTCCCGGCTTTTGCCGCGGCTCCTGGATTACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTc                                                                                                          <  1:797651/150‑1 (MQ=255)
                                            accaGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGt                                                                                                        <  1:717777/150‑1 (MQ=255)
                                              cagcagTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACATGTCGAACACGTCTTCGCCCTCGCGCggg                                                                                                                           <  1:240667/129‑1 (MQ=255)
                                              cagcagTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACATGTCGAACACGTCTTCGCCCTCGCGCggg                                                                                                                           <  1:241076/129‑1 (MQ=255)
                                              cagcagTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACATGTCGAACACGTCTTCGCCCTCGCGCggg                                                                                                                           >  2:240667/1‑129 (MQ=255)
                                              cagcagTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACATGTCGAACACGTCTTCGCCCTCGCGCggg                                                                                                                           >  2:241076/1‑129 (MQ=255)
                                                gcagTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGCCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGgc                                                                                                    >  1:2350/1‑150 (MQ=255)
                                                       gATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCgggg                                                                                                                  >  1:445371/1‑129 (MQ=255)
                                                       gATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCgggg                                                                                                                  <  2:445371/129‑1 (MQ=255)
                                                       gATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAgg                                                                                             >  2:190851/1‑150 (MQ=255)
                                                             ggTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAg                                                                                       >  2:723038/1‑150 (MQ=255)
                                                                 aaGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAgccgc                                                                                   >  2:669815/1‑150 (MQ=255)
                                                                   gCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCt                                                                                 <  2:664183/150‑1 (MQ=255)
                                                                     cTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGCCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGTGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTc                                                                               <  2:2350/150‑1 (MQ=255)
                                                                          ctgcGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAgc                                                                                      <  1:697283/138‑1 (MQ=255)
                                                                          ctgcGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAgc                                                                                      >  2:697283/1‑138 (MQ=255)
                                                                           tgcGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGgcgc                                                                                                  >  1:1187770/1‑125 (MQ=255)
                                                                           tgcGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGgcgc                                                                                                  <  2:1187770/125‑1 (MQ=255)
                                                                              gTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGCAAAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATc                                                                      <  1:976918/150‑1 (MQ=255)
                                                                                          tCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGccac                                                          >  2:665312/1‑150 (MQ=255)
                                                                                          tCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGccac                                                          >  2:203410/1‑150 (MQ=255)
                                                                                            aCCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGccacca                                                        >  2:505120/1‑150 (MQ=255)
                                                                                             ccATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCTGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAg                                                       <  1:1038800/150‑1 (MQ=255)
                                                                                               aTTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCg                                                     <  1:585174/150‑1 (MQ=255)
                                                                                                       tAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGtt                                                                                                                         >  2:1248823/1‑74 (MQ=255)
                                                                                                       tAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGtt                                                                                                                         <  1:1248823/74‑1 (MQ=255)
                                                                                                       tAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGc                                             <  1:485046/150‑1 (MQ=255)
                                                                                                        aaGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAGTACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCg                                            <  1:1136061/150‑1 (MQ=255)
                                                                                                           ggATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGa                                         >  2:236216/1‑150 (MQ=255)
                                                                                                             atatTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATa                                       <  1:723038/150‑1 (MQ=255)
                                                                                                               atTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATAcg                                     >  2:341380/1‑150 (MQ=255)
                                                                                                                ttGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATAcgc                                    >  2:1177610/1‑150 (MQ=255)
                                                                                                                   ttCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGccacca                                                        >  1:32430/1‑127 (MQ=255)
                                                                                                                   ttCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGccacca                                                        <  2:32430/127‑1 (MQ=255)
                                                                                                                      aGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCAACAGCGGGTCGAGCAGGATACGCCGGAgg                              >  2:1037352/1‑150 (MQ=255)
                                                                                                                        gCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGaa                            <  1:505120/150‑1 (MQ=255)
                                                                                                                        gCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGGCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGGTACGCCGCCGGaa                            >  2:1265424/1‑150 (MQ=255)
                                                                                                                           tttGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTa                         >  1:396368/1‑150 (MQ=255)
                                                                                                                           tttGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTa                         <  2:419972/150‑1 (MQ=255)
                                                                                                                              gCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGaa                      <  1:1177610/150‑1 (MQ=255)
                                                                                                                              gCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCACCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGaa                      <  2:827620/150‑1 (MQ=255)
                                                                                                                                 aagagCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACtt                   >  1:44405/5‑150 (MQ=255)
                                                                                                                                       ccTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACTTGGCGcc             <  1:1265424/150‑1 (MQ=255)
                                                                                                                                        cTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACTTGGCGccc            >  1:256958/1‑150 (MQ=255)
                                                                                                                                          ggAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACTTGGCGCCCtt          <  1:669815/150‑1 (MQ=255)
                                                                                                                                            aaTACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACTTGGCGCCCTTgg        <  1:236216/150‑1 (MQ=255)
                                                                                                                                                 acGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACTTGGCGCCCTTGGCACTg   >  2:1079828/1‑150 (MQ=255)
                                                                                                                                                  cGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGTGGCTGCGGTCGTCGGCGCGCAGGCGCAAGCCTAGTGCTTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACTTGGCGCCCTTGGCACTgc  <  1:665312/150‑1 (MQ=255)
                                                                                                                                                  cGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACTTGGCGCCCTTg         <  1:715252/143‑1 (MQ=255)
                                                                                                                                                  cGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACTTGGCGCCCTTg         >  2:715252/1‑143 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACTTGGCGCCCTTGGCACTGC  >  NC_002947/3941505‑3941800

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