Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,176,696 C→T 7.8% Y459Y (TAC→TAT)  PP_3676 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,176,6960CT 7.8% 180.5 / 4.5 65Y459Y (TAC→TAT) PP_3676hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base T (2/3);  total (32/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.17e-01

CTGAAATACCAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGT  >  NC_002947/4176548‑4176836
                                                                                                                                                    |                                                                                                                                            
cTGAAATACCAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTAc                                                                                                                                              <  8:263253/149‑1 (MQ=255)
   aaaTACCAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCg                                                                                                                                            <  6:3713/148‑1 (MQ=255)
       aCCAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAAATATTTCTAACATCGACATCAAAGGGAACGGCGAGTTTCTCAAATTTCATAGGTAAACTTAATTCTACCAAGGGTGTTACGCTACATCAAATTTACCGTGAt                                                                                                                                        >  5:296316/1‑148 (MQ=255)
            gagCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAAttt                                                                                                                                  <  4:232839/149‑1 (MQ=255)
              gcgcGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGa                                                                                                                                <  5:226180/149‑1 (MQ=255)
              gcgcGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGa                                                                                                                                <  7:121968/149‑1 (MQ=255)
                  ggCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTATCACGGTAAATTTGatgta                                                                                                                            <  8:299149/149‑5 (MQ=255)
                      aCTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTATCACGGTAAATTTGatgta                                                                                                                            >  5:242601/1‑141 (MQ=255)
                      aCTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTATCACGGTAAATTTGatgta                                                                                                                            >  1:153427/1‑141 (MQ=255)
                      aCTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTATCACGGTAAATTTGatgta                                                                                                                            <  2:153427/145‑5 (MQ=255)
                      aCTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTATCACGGTAAATTTGatgta                                                                                                                            <  6:242601/145‑5 (MQ=255)
                        tGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTgg                                                                                                                       >  4:55313/1‑148 (MQ=255)
                          cAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCGAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCAt                                                                                                                    >  7:69774/1‑149 (MQ=255)
                             gACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCAc                                                                                                                 >  5:39312/1‑149 (MQ=255)
                               ccGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCAccc                                                                                                               <  1:228730/149‑1 (MQ=255)
                                    actCAGGATGGCGTAAATTTTTTCTATCTTCGACATCAACGGGAACGGCGCGTTTCTCATATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATATATTTGAGGCCGCTGGCATCACCCGCGGa                                                                                                          <  2:89221/147‑1 (MQ=255)
                                    aaTCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGa                                                                                                                                <  1:286170/127‑1 (MQ=255)
                                    aaTCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGa                                                                                                                                >  2:286170/1‑127 (MQ=255)
                                    aaTCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGa                                                                                                          <  2:54857/149‑1 (MQ=255)
                                           aTGGCGTAAATTATTTCTAACATCGACATCAAAGGGTACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCa                                                                                                   <  6:296316/149‑1 (MQ=255)
                                             ggCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCAc                                                                                                  >  3:44001/1‑148 (MQ=255)
                                               cGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCt                                                                                               <  7:241739/149‑1 (MQ=255)
                                                           cTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTc                                                                                    >  4:110895/1‑148 (MQ=255)
                                                             aaCATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTTATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGAt                                                                                 <  6:39312/149‑1 (MQ=255)
                                                                 tCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGca                                                                              >  8:169400/1‑148 (MQ=255)
                                                                    aCATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGcacaca                                                                          <  7:169400/149‑1 (MQ=255)
                                                                       tCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAggg                                                                       >  2:262141/1‑149 (MQ=255)
                                                                         aaaGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGAt                                                                                 <  2:319555/137‑1 (MQ=255)
                                                                         aaaGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGAt                                                                                 >  1:319555/1‑137 (MQ=255)
                                                                           aGGGAACGGCGCGTTTCTCAAATTTCATAGTTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCAtt                                                                   <  4:44001/149‑1 (MQ=255)
                                                                           aGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCt                                                                                               >  4:170065/1‑121 (MQ=255)
                                                                           aGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCt                                                                                               <  3:170065/121‑1 (MQ=255)
                                                                             ggAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTc                                                                  >  5:244385/1‑148 (MQ=255)
                                                                             ggAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTc                                                                  >  1:255111/1‑148 (MQ=255)
                                                                                  gtcgcgTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAAt                                                            <  6:281903/147‑1 (MQ=255)
                                                                                     gcgTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAAtttt                                                         >  5:165724/1‑149 (MQ=255)
                                                                                           ctaaaaTTTCATAGGTAAACTTAATTCTTAAAAGGGTGTTAAGCTACATCAAATTTAACGTGATAAATTTTAGGCCGCTGGCATTAACCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTATAGAATTTTTACgg                                                    >  2:67402/4‑148 (MQ=255)
                                                                                            tcAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGaa                                                  >  3:73351/1‑149 (MQ=255)
                                                                                            tcAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGaa                                                  >  5:162816/1‑149 (MQ=255)
                                                                                                 tttCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAAt                                                            >  3:235965/1‑134 (MQ=255)
                                                                                                 tttCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAAt                                                            <  4:235965/134‑1 (MQ=255)
                                                                                                       aGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCATACAGGGCATTCTAGAATTTTTACGGAACGGCAGCAtt                                        >  2:28362/1‑148 (MQ=255)
                                                                                                        ggTAAACTTAATTCTTCGAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTc                                       <  8:69774/148‑1 (MQ=255)
                                                                                                         gTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAActg                                               >  8:14383/1‑137 (MQ=255)
                                                                                                         gTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACgg                                               <  7:14383/139‑1 (MQ=255)
                                                                                                            aaCTTAATTCTTCCAATGTTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATTACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGTGCATTCTATAATTTTTACGGAACGGCAGCATTCAAaga                                  <  1:67402/149‑1 (MQ=255)
                                                                                                            aaCTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCATACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAag                                   <  1:28362/148‑1 (MQ=255)
                                                                                                            aaCTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCaaa                                    >  6:193315/1‑147 (MQ=255)
                                                                                                                      ttcCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAAc                                                 <  3:156720/124‑1 (MQ=255)
                                                                                                                      ttcCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAAc                                                 >  4:156720/1‑124 (MQ=255)
                                                                                                                        ccAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGAc                         <  5:193315/146‑1 (MQ=255)
                                                                                                                         cAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGcacac                                                                           <  6:294885/95‑1 (MQ=255)
                                                                                                                         cAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGcacac                                                                           >  5:294885/1‑95 (MQ=255)
                                                                                                                          aaGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACt                    >  4:239693/1‑149 (MQ=255)
                                                                                                                           aGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTg                   <  3:239693/149‑1 (MQ=255)
                                                                                                                                gtTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGt                                                                                        <  3:266853/75‑1 (MQ=255)
                                                                                                                                gtTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGt                                                                                        >  4:266853/1‑75 (MQ=255)
                                                                                                                                      cTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATg        <  3:208041/149‑1 (MQ=255)
                                                                                                                                          aTCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGcc      >  7:228645/1‑147 (MQ=255)
                                                                                                                                          aTCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCaaa    >  4:207296/1‑149 (MQ=255)
                                                                                                                                          aTCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCaa    <  3:110895/149‑1 (MQ=255)
                                                                                                                                           tCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATTCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAg   <  6:244385/149‑1 (MQ=255)
                                                                                                                                           tCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAg   >  6:114597/1‑149 (MQ=255)
                                                                                                                                            cAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGt  >  6:125948/1‑149 (MQ=255)
                                                                                                                                               aTTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCTCGGAGTGCTCACCCTATTCCGTATTCGATTCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGc       <  1:52472/141‑1 (MQ=255)
                                                                                                                                               aTTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGc       >  2:52472/1‑141 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CTGAAATACCAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGT  >  NC_002947/4176548‑4176836

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