Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,750,398 C→T 7.1% V389M (GTG→ATG)  thiI ← tRNA sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,750,3980CT 7.1% 249.0 / 6.5 84V389M (GTG→ATG) thiItRNA sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (39/39);  new base T (2/4);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

GAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTG  >  NC_002947/5750252‑5750546
                                                                                                                                                   |                                                                                                                                                    
gAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTCGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATcaca                                                                                                                                                     <  7:253981/148‑1 (MQ=255)
  aCGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGa                                                                                                                                                  >  2:114612/1‑149 (MQ=255)
      cATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCCCATCGATCATGTGATCGa                                                                                                                                              <  5:312458/149‑1 (MQ=255)
       aTGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGAt                                                                                                                                             <  8:244212/149‑1 (MQ=255)
        tGGCCTGGAATTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCATGTGATCGATg                                                                                                                                            <  3:70780/149‑1 (MQ=255)
          gCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATgg                                                                                                                                           <  8:100500/148‑1 (MQ=255)
               gACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATg                                                                                                                                      <  3:221309/148‑1 (MQ=255)
               gACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATg                                                                                                                                      >  1:111379/1‑148 (MQ=255)
                aCTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGaa                                                                                                                                    >  5:309708/1‑149 (MQ=255)
                  ttCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACt                                                                                                                                  >  1:205227/1‑149 (MQ=255)
                     gATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTgg                                                                                                                               >  6:72863/1‑149 (MQ=255)
                     gATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTg                                                                                                                                <  6:313841/148‑1 (MQ=255)
                     gATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTg                                                                                                                                <  5:124048/148‑1 (MQ=255)
                     gATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCACGGCCACTTCCTCGATTTCGATGACCTTGCCAAGCTCAACGATCACATGATCGATGGAAATGAACtagt                                                                                                                               >  7:291964/1‑146 (MQ=255)
                      aTACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGgc                                                                                                                              <  3:260703/149‑1 (MQ=255)
                          cTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTg                                                                                                                                >  8:251879/1‑143 (MQ=255)
                          cTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTg                                                                                                                                <  7:251879/143‑1 (MQ=255)
                          cTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGgcgcgc                                                                                                                          >  3:79668/1‑149 (MQ=255)
                             ccAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTc                                                                                                                        >  6:240352/1‑148 (MQ=255)
                               aGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTc                                                                                                                        >  2:124606/1‑146 (MQ=255)
                               aGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTc                                                                                                                        >  3:63816/1‑146 (MQ=255)
                                gCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCATGTGATCGAt                                                                                                                                             >  4:25779/1‑124 (MQ=255)
                                gCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCATGTGATCGAt                                                                                                                                             <  3:25779/124‑1 (MQ=255)
                                gCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAaga                                                                                                                    <  5:72863/149‑1 (MQ=255)
                                    cAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAgcgc                                                                                                                >  3:216059/1‑149 (MQ=255)
                                      gagGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCACTCAAGAgcgcgc                                                                                                              >  1:128229/1‑149 (MQ=255)
                                         gTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTcc                                                                                                           >  7:106478/1‑149 (MQ=255)
                                          ttCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTcc                                                                                                           <  2:111379/148‑1 (MQ=255)
                                             gtattCCTTGGCTTCGGGGTGTCGAATGTCTATGACGATCTGGCCTGGGAGCGCCTCTGCCAATTCCTCTATTTCGATGTCCTTGCCCAGCTCATCTATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGc                                                                                                        <  6:42583/145‑1 (MQ=255)
                                                  ccTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGc                                                                                                   >  1:90543/1‑148 (MQ=255)
                                                     gggCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATg                                                                                               >  7:181219/1‑149 (MQ=255)
                                                      ggCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATg                                                                                               <  1:124606/148‑1 (MQ=255)
                                                         aTCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGaa                                                                                                                                    >  6:273978/1‑108 (MQ=255)
                                                         aTCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGaa                                                                                                                                    <  5:273978/108‑1 (MQ=255)
                                                             gggTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTg                                                                                       >  2:154656/1‑149 (MQ=255)
                                                               gTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCt                                                                                     <  2:205227/149‑1 (MQ=255)
                                                                     aaTGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCAGGTCGAACTGCTTttctt                                                                               >  8:311352/1‑149 (MQ=255)
                                                                      aTGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTttcttc                                                                              <  8:106478/149‑1 (MQ=255)
                                                                      aTGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTttcttc                                                                              <  5:189907/149‑1 (MQ=255)
                                                                       tGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGTTCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCg                                                                             >  7:195272/1‑149 (MQ=255)
                                                                        gTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCTCCCCATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTttcttc                                                                              <  1:114612/148‑1 (MQ=255)
                                                                          cGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACg                                                                                                     <  3:133461/122‑1 (MQ=255)
                                                                           gATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCt                                                                         <  8:181219/149‑1 (MQ=255)
                                                                              gACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCa                                                                      >  3:20666/1‑149 (MQ=255)
                                                                                       gcctggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCg                                                                  <  3:34070/144‑1 (MQ=255)
                                                                                       gcctggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGt                                                                <  4:63816/146‑1 (MQ=255)
                                                                                       gcctggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGt                                                                <  6:19748/146‑1 (MQ=255)
                                                                                         ctggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCACTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGc                                                           <  2:128229/149‑1 (MQ=255)
                                                                                            gcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTg                                                        <  5:240352/149‑1 (MQ=255)
                                                                                               gcgcCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCa                                                     <  6:207546/149‑1 (MQ=255)
                                                                                                  ccTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATgg                                                  <  3:75625/149‑1 (MQ=255)
                                                                                                  ccTCGGCCACTTCCTAGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATgg                                                  <  7:267921/149‑1 (MQ=255)
                                                                                                   cTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCgg                                                                 >  4:294923/1‑133 (MQ=255)
                                                                                                   cTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCgg                                                                 <  3:294923/133‑1 (MQ=255)
                                                                                                      ggCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGggtgg                                              >  2:139302/1‑149 (MQ=255)
                                                                                                          aCTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCgg                                          <  4:20666/149‑1 (MQ=255)
                                                                                                          aCTTACTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATAGATGGAAATGAACTTGGCGGGCTGAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCgg                                          <  8:291964/149‑1 (MQ=255)
                                                                                                             tCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGtt                                       >  3:322672/1‑149 (MQ=255)
                                                                                                                cTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATgg                                                  <  6:88530/136‑1 (MQ=255)
                                                                                                                cTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGAGATCggaag                                       >  5:88530/1‑144 (MQ=255)
                                                                                                                 tCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTAc                                    >  3:137420/1‑149 (MQ=255)
                                                                                                                   gATTTCGATGTCCTTGCCCAGCTCATCGTTCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTAcc                                   <  8:195272/148‑1 (MQ=255)
                                                                                                                   gATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTAcc                                   >  1:238430/1‑148 (MQ=255)
                                                                                                                    aTTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGggtgg                                              <  4:159927/136‑1 (MQ=255)
                                                                                                                    aTTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGggtgg                                              >  3:159927/1‑136 (MQ=255)
                                                                                                                    aTTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCga                                 <  1:98509/149‑1 (MQ=255)
                                                                                                                         gATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTc                                                                        <  2:184983/105‑1 (MQ=255)
                                                                                                                         gATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTc                                                                        >  1:184983/1‑105 (MQ=255)
                                                                                                                         gATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCa                            >  4:282835/1‑149 (MQ=255)
                                                                                                                         gATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCa                            >  6:76131/1‑149 (MQ=255)
                                                                                                                          aTGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCgg                                          <  8:57382/134‑1 (MQ=255)
                                                                                                                          aTGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCgg                                          >  7:57382/1‑134 (MQ=255)
                                                                                                                          aTGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATCGGTGGTCGGGTTTACCGAGATCAc                           >  6:240870/1‑149 (MQ=255)
                                                                                                                               cTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCAcgccg                       >  7:10339/1‑148 (MQ=255)
                                                                                                                                   cccAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTa                  <  8:10339/149‑1 (MQ=255)
                                                                                                                                       gCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTcc              >  3:310831/1‑149 (MQ=255)
                                                                                                                                         tCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCGTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCgg            >  4:20525/1‑149 (MQ=255)
                                                                                                                                          cATCGATCATGTGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCgg            >  4:255682/1‑148 (MQ=255)
                                                                                                                                           aTCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCa          >  7:68113/1‑149 (MQ=255)
                                                                                                                                           aTCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCa          >  3:324935/1‑149 (MQ=255)
                                                                                                                                            tCGATCATGTGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCAt         <  3:255682/149‑1 (MQ=255)
                                                                                                                                            tCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCAt         <  7:96636/149‑1 (MQ=255)
                                                                                                                                            tCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCAt         <  6:309708/149‑1 (MQ=255)
                                                                                                                                            tCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCAt         <  1:119519/149‑1 (MQ=255)
                                                                                                                                            tCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCAt         >  7:100774/1‑149 (MQ=255)
                                                                                                                                                  acaTGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCtt   <  2:90543/149‑1 (MQ=255)
                                                                                                                                                   catcatCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTg  <  4:216059/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTG  >  NC_002947/5750252‑5750546

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