Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,339,251 C→T 5.7% intergenic (‑210/+33) PP_2935 ← / ← ybiT membrane protein/ABC‑F family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,339,2510CT 5.7% 206.3 / 2.6 70intergenic (‑210/+33)PP_2935/ybiTmembrane protein/ABC‑F family ATPase
Reads supporting (aligned to +/- strand):  ref base C (32/34);  new base T (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

CAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCG  >  NC_002947/3339108‑3339397
                                                                                                                                               |                                                                                                                                                   
cagGCTAATGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCAt                                                                                                                                               <  1:660789/150‑1 (MQ=255)
    cTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCga                                                                                                                                          >  1:118288/1‑150 (MQ=255)
    cTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCga                                                                                                                                          <  1:1402257/150‑1 (MQ=255)
     tAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCTGGCCATctcgg                                                                                                                                          >  1:1246729/1‑146 (MQ=255)
     tAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCgag                                                                                                                                         >  1:925891/1‑150 (MQ=255)
      aaCGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCgaga                                                                                                                                        >  1:82005/1‑150 (MQ=255)
      aaCGAATTATTCGTCGGGCCGTGCCGGTCGGTGAGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGACCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCgaga                                                                                                                                        >  1:474280/1‑150 (MQ=255)
        cGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTATTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCCAGCCCGCGAAGAGGCCCGGCCATCCGAGATg                                                                                                                                      <  1:1313090/150‑1 (MQ=255)
               ttCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGc                                                                                                                               >  1:1161496/1‑150 (MQ=255)
                tcgtcgGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCTGGCCATCTCGGATGGCCGGGcc                                                                                                                              <  2:641432/150‑1 (MQ=38)
                     gggCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCTGGCCATCTCGGATGGCCGGGcctcttc                                                                                                                         <  1:1248120/150‑6 (MQ=38)
                              cggtcggtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATc                                                                                                                                              <  1:1100288/121‑1 (MQ=255)
                              cggtcggtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATc                                                                                                                                              >  2:1100288/1‑121 (MQ=255)
                              cggtcggtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAg                                                                                                                <  2:92957/150‑1 (MQ=255)
                                gtcggtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGAGGGCACGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGAc                                                                                                              <  2:953624/150‑1 (MQ=255)
                                gtcggtGGGATGGGGATGTACTCTTCGCGGCCTCGCTCGCGCGCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGAc                                                                                                              <  2:474280/150‑1 (MQ=39)
                                  cggtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGTGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCa                                                                                                            <  2:199080/150‑1 (MQ=255)
                                        gatggggatgTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTTCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCacac                                                                                                      >  2:695300/1‑150 (MQ=255)
                                              gatgTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGATGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGc                                                                                                <  1:1185722/150‑1 (MQ=255)
                                                 gTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTAc                                                                                             >  2:622245/1‑150 (MQ=255)
                                                  tACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACg                                                                                            >  1:897225/1‑150 (MQ=255)
                                                     tcctCGCGGGCTTCACCGCTCGCCCAAGGCCACAAAAGCGTCCAACACCTGTGCGGTCCCTGTGGGCGCGTCCAAGCCCGCGAAGAGGCCCGCCCATCCGCGATGGCCAGGCCGTCCCATCAATCCGACCAGCACACCCTGGCTACGCAg                                                                                         <  2:431855/147‑1 (MQ=255)
                                                       TTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCga                                                                                                                                          >  1:668497‑M1/1‑99 (MQ=255)
                                                       TTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGt                                                                                       >  1:464308‑M1/1‑150 (MQ=255)
                                                       TTCGCGGGCTTGCCCGCGCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCga                                                                                                                                          <  2:668497‑M1/99‑1 (MQ=255)
                                                        TCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTa                                                                                      >  1:622245‑M1/1‑150 (MQ=255)
                                                          GCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGc                                                                                                                               >  2:815347‑M1/1‑107 (MQ=255)
                                                          GCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGc                                                                                                                               <  1:815347‑M1/107‑1 (MQ=255)
                                                          GCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGt                                                                                    <  2:962933‑M1/150‑1 (MQ=255)
                                                            GGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGAt                                                                                                                                       >  2:774690‑M1/1‑97 (MQ=255)
                                                            GGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGAt                                                                                                                                       <  1:774690‑M1/97‑1 (MQ=255)
                                                                TTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGt                                                                              <  2:1268352‑M1/150‑1 (MQ=255)
                                                                TTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGt                                                                              >  1:524476‑M1/1‑150 (MQ=255)
                                                                  GCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAg                                                                            >  2:1161496‑M1/1‑150 (MQ=255)
                                                                       CTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTAcc                                                                       <  2:842267‑M1/150‑1 (MQ=255)
                                                                         cccacagggaccgcacaggtattggacgcttgtggggccttgggagagcgggcaagcccgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGc                                                                     >  1:782677‑M2/68‑150 (MQ=255)
                                                                           cccAGGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTg                                                                   <  2:932275/150‑1 (MQ=255)
                                                                            ccAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGAc                                                                                                              <  1:489567/106‑1 (MQ=255)
                                                                            ccAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGAc                                                                                                              >  2:489567/1‑106 (MQ=255)
                                                                            ccAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGa                                                                  >  1:406080/1‑150 (MQ=255)
                                                                            ccAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGa                                                                  >  2:183380/1‑150 (MQ=255)
                                                                             cAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGaa                                                                 >  2:1406319/1‑150 (MQ=255)
                                                                              aaGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAg                                                                >  2:1186145/1‑150 (MQ=255)
                                                                                  ccccACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGTGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTcca                                                            <  1:1185266/150‑1 (MQ=255)
                                                                                   cccACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTccac                                                           <  1:995813/150‑1 (MQ=255)
                                                                                   cccACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTACCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTccac                                                           >  1:1240671/1‑150 (MQ=255)
                                                                                      cacaggtattggacgcttgtggggccttgggagagcgggcaagcccgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGa                                                                                                               <  2:99140‑M2/41‑1 (MQ=255)
                                                                                      cacaggtattggacgcttgtggggccttgggagagcgggcaagcccgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGa                                                                                                               >  1:99140‑M2/55‑95 (MQ=255)
                                                                                       caAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCcaccaca                                                       <  1:599100/150‑1 (MQ=255)
                                                                                               cAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGa                                                                                                               >  2:983525/1‑86 (MQ=255)
                                                                                               cAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGa                                                                                                               <  1:983525/86‑1 (MQ=255)
                                                                                               cAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAg                                                                                                           <  2:302100/90‑1 (MQ=255)
                                                                                               cAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAg                                                                                                           >  1:302100/1‑90 (MQ=255)
                                                                                                caTACCTTTGCGGTCCCTGTGGGAGCGGGCAATCCCTCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGcc                                              <  2:82005/149‑1 (MQ=255)
                                                                                                 aTACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGCCCACCACACCATCGGCCg                                             >  2:831392/1‑150 (MQ=255)
                                                                                                       gtgCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGCCCACCACACCATCGGCCGACAGCt                                       <  1:831392/150‑1 (MQ=255)
                                                                                                          cGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCga                                                                                                                                          >  2:1402111/1‑48 (MQ=255)
                                                                                                          cGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCga                                                                                                                                          <  1:1402111/48‑1 (MQ=255)
                                                                                                          cGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGc                                                                     >  1:1354194/1‑117 (MQ=255)
                                                                                                          cGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGc                                                                     <  2:1354194/117‑1 (MQ=255)
                                                                                                              cccTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCgatgat                                <  2:118288/150‑1 (MQ=255)
                                                                                                               ccTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCTGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGatgatg                               >  2:1043265/1‑150 (MQ=255)
                                                                                                               ccTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGaa                                                                 <  1:1203319/116‑1 (MQ=255)
                                                                                                               ccTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGaa                                                                 >  2:1203319/1‑116 (MQ=255)
                                                                                                                     ggAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCATACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTg                         <  1:675682/150‑1 (MQ=255)
                                                                                                                                gcccgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTACCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGcgacga              <  2:1240671‑M2/138‑1 (MQ=255)
                                                                                                                                 cccgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCgacgac             >  1:1134999‑M2/12‑150 (MQ=255)
                                                                                                                                 cccgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGAAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCgacgac             >  1:807386‑M2/12‑150 (MQ=255)
                                                                                                                                            gCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTcgcg  <  1:695300/150‑1 (MQ=255)
                                                                                                                                            gCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTcgcg  <  2:925891/150‑1 (MQ=255)
                                                                                                                                             cccGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGc                       >  2:1098625/1‑128 (MQ=255)
                                                                                                                                             cccGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGc                       <  1:1098625/128‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCG  >  NC_002947/3339108‑3339397

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