Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,176,701 A→G 7.7% D461G (GAT→GGT)  PP_3676 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,176,7010AG 7.7% 184.3 / 3.2 64D461G (GAT→GGT) PP_3676hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base G (2/3);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGT  >  NC_002947/4176555‑4176849
                                                                                                                                                  |                                                                                                                                                    
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                                                                                                                                         <  6:242601/145‑5 (MQ=255)
               aCTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTATCACGGTAAATTTGatgta                                                                                                                                         <  2:153427/145‑5 (MQ=255)
               aCTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTATCACGGTAAATTTGatgta                                                                                                                                         >  1:153427/1‑141 (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                                                                                                                                             >  2:286170/1‑127 (MQ=255)
                             aaTCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGa                                                                                                                                             <  1:286170/127‑1 (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                                                                                              >  1:319555/1‑137 (MQ=255)
                                                                  aaaGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGAt                                                                                              <  2:319555/137‑1 (MQ=255)
                                                                    aGGGAACGGCGCGTTTCTCAAATTTCATAGTTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCAtt                                                                                <  4:44001/149‑1 (MQ=255)
                                                                    aGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCt                                                                                                            >  4:170065/1‑121 (MQ=255)
                                                                    aGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCt                                                                                                            <  3:170065/121‑1 (MQ=255)
                                                                      ggAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTc                                                                               >  1:255111/1‑148 (MQ=255)
                                                                      ggAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTc                                                                               >  5:244385/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                                                                                        >  5:294885/1‑95 (MQ=255)
                                                                                                                  cAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGcacac                                                                                        <  6:294885/95‑1 (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)
                                                                                                                                                  aTAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGt  <  6:162816/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ACCAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGT  >  NC_002947/4176555‑4176849

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: ‑