Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,694,072 T→G 5.0% R713R (CGA→CGC ligD ← DNA ligase D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,694,0720TG 5.0% 349.4 / 5.9 118R713R (CGA→CGCligDDNA ligase D
Reads supporting (aligned to +/- strand):  ref base T (52/60);  new base G (4/2);  total (56/62)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 120 total aligned reads displayed.

CGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGCAGGGTCAGGGTCCAG  >  NC_002947/3693927‑3694218
                                                                                                                                                 |                                                                                                                                                  
cgcCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGATGCCCCCGTTTACCTCGTCCCACCCGGCTCGcc                                                                                                                                                <  1:1617833/150‑1 (MQ=255)
 gcCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCg                                                                                                                                               >  1:1489774/1‑150 (MQ=255)
 gcCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCg                                                                                                                                               >  2:2790738/1‑150 (MQ=255)
         cGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAgg                                                                                                                                       <  1:1117996/150‑1 (MQ=255)
            tGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGca                                                                                                                                    <  2:2426601/150‑1 (MQ=255)
               gcAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCATGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGCGCCGTAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACaa                                                                                                                                 <  1:2859646/150‑1 (MQ=255)
               gcAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACaa                                                                                                                                 >  1:376413/1‑150 (MQ=255)
               gcAGGTAGGCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCCCCCGGCTCGCCGGGTCAAGGGCACaa                                                                                                                                 <  2:1312451/150‑1 (MQ=255)
                   gTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCgg                                                                                                                                              >  2:1699525/1‑133 (MQ=255)
                     attctaTGAAAATGCGTCCCACGCGATTTTTCGGCCCCCACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGt                                                                                                                           <  1:113039/145‑1 (MQ=255)
                      gTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTg                                                                                                                          <  1:2346637/150‑1 (MQ=255)
                       tCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGa                                                                                                                         >  1:1505518/1‑150 (MQ=255)
                        cGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGaa                                                                                                                        >  2:1984891/1‑150 (MQ=255)
                        cGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGaa                                                                                                                        <  2:282793/150‑1 (MQ=255)
                         gATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAAt                                                                                                                       >  1:1013747/1‑150 (MQ=255)
                               aaTGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCtt                                                                                                                 <  1:734638/150‑1 (MQ=255)
                               aaTGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCtt                                                                                                                 <  1:2440466/150‑1 (MQ=255)
                                 tGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGc                                                                                                               >  2:2166460/1‑150 (MQ=255)
                                   cGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGcca                                                                                                             <  2:2105774/150‑1 (MQ=255)
                                    gTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAggg                                                                                                                                      >  1:1866549/1‑124 (MQ=255)
                                    gTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAggg                                                                                                                                      <  2:1866549/124‑1 (MQ=255)
                                    gTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGccac                                                                                                            >  2:1433076/1‑150 (MQ=255)
                                    gTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGccac                                                                                                            <  2:376413/150‑1 (MQ=255)
                                     tCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGccacc                                                                                                           >  1:2731085/1‑150 (MQ=255)
                                      cccACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGccacca                                                                                                          <  2:1013747/150‑1 (MQ=255)
                                       ccACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCcaccac                                                                                                         >  1:2731974/1‑150 (MQ=255)
                                           gcgATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTg                                                                                                       <  1:481814/148‑1 (MQ=255)
                                           gcgATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTg                                                                                                       >  2:481814/1‑148 (MQ=255)
                                           gcgATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGt                                                                                                     >  2:1746114/1‑150 (MQ=255)
                                           gcgATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCCGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGt                                                                                                     <  1:863576/150‑1 (MQ=255)
                                              aTTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGc                                                                                                               <  2:1985583/137‑1 (MQ=255)
                                              aTTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGc                                                                                                               >  1:1985583/1‑137 (MQ=255)
                                              aTTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGtttt                                                                                                  >  2:1063554/1‑150 (MQ=255)
                                              aTTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGtttt                                                                                                  <  1:2032578/150‑1 (MQ=255)
                                              aTTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGtttt                                                                                                  <  2:2230931/150‑1 (MQ=255)
                                                 tttCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAg                                                                                               <  2:1674685/150‑1 (MQ=255)
                                                 tttCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAg                                                                                               <  2:1505518/150‑1 (MQ=255)
                                                    cGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCAACACTGGTTTTCAGGaa                                                                                            <  1:2581887/150‑1 (MQ=255)
                                                      gCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAaca                                                                                          >  2:1165406/1‑150 (MQ=255)
                                                       ccccGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAacac                                                                                         <  1:198137/150‑1 (MQ=255)
                                                        cccGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACAcc                                                                                        >  1:2334472/1‑150 (MQ=255)
                                                        cccGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACAcc                                                                                        >  2:756302/1‑150 (MQ=255)
                                                          cGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCtt                                                                                      >  2:2563154/1‑150 (MQ=255)
                                                           gACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCCCGGAAAGTGCCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTc                                                                                     <  1:2091508/150‑1 (MQ=255)
                                                            acacCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCcaccac                                                                                                         >  1:2705888/1‑129 (MQ=255)
                                                            acacCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCcaccac                                                                                                         <  2:2705888/129‑1 (MQ=255)
                                                            acacCCCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCa                                                                                    <  1:1433076/150‑1 (MQ=255)
                                                                cgcgCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGcc                                                                                >  2:2058288/1‑150 (MQ=255)
                                                                      cAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTc                                                                          <  1:1063554/150‑1 (MQ=255)
                                                                               cGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCgg                                                                                                                                              >  1:2928433/1‑73 (MQ=255)
                                                                               cGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCgg                                                                                                                                              <  2:2928433/73‑1 (MQ=255)
                                                                               cGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAgg                                                                                              <  2:1434536/121‑1 (MQ=255)
                                                                               cGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCa                                                                                    <  2:2264841/131‑1 (MQ=255)
                                                                               cGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCa                                                                                    >  1:2264841/1‑131 (MQ=255)
                                                                                     cAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAg                                                           <  2:2731974/150‑1 (MQ=255)
                                                                                     cAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAg                                                           <  1:1963164/150‑1 (MQ=255)
                                                                                      aGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTTAgg                                                          <  1:756302/150‑1 (MQ=255)
                                                                                      aGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAgg                                                          <  1:167602/150‑1 (MQ=255)
                                                                                         ttGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTc                                                       >  1:819872/1‑150 (MQ=255)
                                                                                             ccAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCgg                                                                                                                                              >  2:876655/1‑59 (MQ=255)
                                                                                             ccAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCgg                                                                                                                                              <  1:876655/59‑1 (MQ=255)
                                                                                             ccAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGtt                                                   <  1:2166460/150‑1 (MQ=255)
                                                                                              cAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTg                                                  >  2:1058910/1‑150 (MQ=255)
                                                                                                     ttGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTggg                                                >  1:3069752/1‑145 (MQ=255)
                                                                                                     ttGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGcc                                           <  1:2058288/150‑1 (MQ=255)
                                                                                                      tGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAGTACCCTTGCCACCACTGGtttt                                                                                                  >  1:477908/1‑94 (MQ=255)
                                                                                                      tGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAGTACCCTTGCCACCACTGGtttt                                                                                                  <  2:477908/94‑1 (MQ=255)
                                                                                                            tGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTccacca                                    <  1:1984891/150‑1 (MQ=255)
                                                                                                              gCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATg                                  <  1:2563154/150‑1 (MQ=255)
                                                                                                                gTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGa                                                                                                                         >  1:1670305/1‑61 (MQ=255)
                                                                                                                gTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGa                                                                                                                         <  2:1670305/61‑1 (MQ=255)
                                                                                                                gTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGc                                >  1:279314/1‑150 (MQ=255)
                                                                                                                    cTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGc                                            <  1:1711606/134‑1 (MQ=255)
                                                                                                                     tGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTcc                           >  1:1732770/1‑150 (MQ=255)
                                                                                                                     tGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTcc                           <  2:465802/150‑1 (MQ=255)
                                                                                                                     tGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTcc                           >  2:1719640/1‑150 (MQ=255)
                                                                                                                          tCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCATGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGc                      >  2:2029835/1‑150 (MQ=255)
                                                                                                                            cTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAg                    >  1:3049552/1‑150 (MQ=255)
                                                                                                                               gACCTCGTCCCACCCGGATCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAgcgc                 >  2:274492/1‑150 (MQ=255)
                                                                                                                                cccTCGTCCCACCCGGCTCGCCGCGTCAAGGGCACCACCACGTGCATACCCTAGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTTGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGca                <  2:1536613/149‑1 (MQ=255)
                                                                                                                                aCCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGca                <  2:1765341/150‑1 (MQ=255)
                                                                                                                                  cTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagg              >  2:2257482/1‑150 (MQ=255)
                                                                                                                                     ccttgaCCCGGCGAGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGcc                               >  1:1703605/6‑130 (MQ=255)
                                                                                                                                     ccttgaCCCGGCGAGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGcc                               <  2:1703605/125‑1 (MQ=255)
                                                                                                                                      tCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAAt                                                                                                                       <  1:2308323/41‑1 (MQ=255)
                                                                                                                                      tCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAAt                                                                                                                       >  2:2308323/1‑41 (MQ=255)
                                                                                                                                      tCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTgg                                                                                                      >  1:1760324/1‑58 (MQ=255)
                                                                                                                                      tCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTgg                                                                                                      <  2:1760324/58‑1 (MQ=255)
                                                                                                                                      tCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGc                                            >  1:2147760/1‑116 (MQ=255)
                                                                                                                                      cttgaCCCGGCGAGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGa                                                            <  1:494920/96‑1 (MQ=255)
                                                                                                                                      cttgaCCCGGCGAGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGa                                                            >  2:494920/5‑100 (MQ=255)
                                                                                                                                      cttgaCCCGGCGAGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggtca          >  2:759639/5‑150 (MQ=255)
                                                                                                                                       cccACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggtcag         >  2:2390143/1‑150 (MQ=255)
                                                                                                                                         cACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggtcaggg       >  2:1685932/1‑150 (MQ=255)
                                                                                                                                         cACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggtcaggg       <  2:616325/150‑1 (MQ=255)
                                                                                                                                         cACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggtcaggg       <  1:1719640/150‑1 (MQ=255)
                                                                                                                                          aCCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCg                                                                         <  2:319491/83‑1 (MQ=255)
                                                                                                                                            ccGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGCAGGGTCAGGGTcc    >  1:959218/1‑150 (MQ=255)
                                                                                                                                              ggCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGCAGGGTCAGGGTCCAg  <  2:2496533/150‑1 (MQ=255)
                                                                                                                                                 tCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagggtcagg        >  1:493096/1‑141 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGCAGGGTCAGGGTCCAG  >  NC_002947/3693927‑3694218

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: