Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,427,281 G→A 9.4% I274I (ATC→ATT)  rhlE‑I ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,427,2810GA 9.4% 180.6 / 7.0 64I274I (ATC→ATT) rhlE‑IATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base G (28/30);  new base A (3/3);  total (31/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.55e-01

GATCGGCAGGTCGAGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAG  >  NC_002947/5427134‑5427429
                                                                                                                                                   |                                                                                                                                                    
gATCGGCAGGTCGAGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGa                                                                                                                                                     >  6:240990/1‑149 (MQ=255)
 aTCGGCAGGTCGAGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGAt                                                                                                                                                    <  7:131265/149‑1 (MQ=255)
     gCAGGTCGAGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGAATGcc                                                                                                                                                 >  3:52648/1‑148 (MQ=255)
     gCAGGTCGAGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTAACGGGCCTTGAAGCTGTCCAGCGAACGCTGGCGTGTTGCCTGTGGGCGGTCACCGTGGATGCCg                                                                                                                                                >  4:260046/1‑149 (MQ=255)
      cAGGTCGAGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGt                                                                                                                                               <  6:6478/149‑1 (MQ=255)
      cAGGTCGAGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGt                                                                                                                                               >  3:85332/1‑149 (MQ=255)
                ttGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCAc                                                                                                                                     <  7:129961/149‑1 (MQ=255)
                ttGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCAc                                                                                                                                     <  1:198129/149‑1 (MQ=255)
                   acgacCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGc                                                                                                                                   <  4:64383/148‑1 (MQ=255)
                   acgacCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGAATGCCGACGGCATCCACgg                                                                                                                                   <  4:52648/148‑2 (MQ=255)
                       ccAGCGGCAAGTCATCGATATCCAGGGCACTGTCCGGCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCg                                                                                                                              <  6:239010/149‑1 (MQ=255)
                        cAGCTGCAAGTCATCGTTATGCAGGCCTCGTGCCTGCACATCTGTCGCCAGCAGCTCTTGCACCTCTCGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGTGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCg                                                                                                                              <  7:226710/148‑1 (MQ=255)
                               aaGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCaa                                                                                                                      <  8:170549/149‑1 (MQ=255)
                                  tCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCg                                                                                                                              <  4:233854/138‑1 (MQ=255)
                                  tCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCg                                                                                                                              >  3:233854/1‑138 (MQ=255)
                                  tCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAgcg                                                                                                                   >  4:232509/1‑149 (MQ=255)
                                      cGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTc                                                                                                                >  8:86885/1‑148 (MQ=255)
                                      cGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTc                                                                                                                >  2:298968/1‑148 (MQ=255)
                                        atatCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTccac                                                                                                             <  4:9892/149‑1 (MQ=255)
                                           tCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTGGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAg                                                                                                          >  5:306828/1‑149 (MQ=255)
                                           tCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCGTCGGCCAGCAAGCGCTCCACCAg                                                                                                          >  8:47439/1‑149 (MQ=255)
                                            ccAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGc                                                                                                         >  7:65418/1‑149 (MQ=255)
                                             cAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCt                                                                                                        >  5:138671/1‑149 (MQ=255)
                                              aGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCt                                                                                                        >  5:251621/1‑148 (MQ=255)
                                                gCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTccac                                                                                                             <  3:14226/141‑1 (MQ=255)
                                                gCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTccac                                                                                                             >  4:14226/1‑141 (MQ=255)
                                                   aCGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGa                                                                                                  >  3:25070/1‑149 (MQ=255)
                                                           ccACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTg                                                                                          >  8:7053/1‑149 (MQ=255)
                                                            cacaTCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGc                                                                                         <  8:65418/149‑1 (MQ=255)
                                                              caTCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCgg                                                                                       <  4:25070/149‑1 (MQ=255)
                                                               aTCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTccacc                                                                                                            <  3:185918/127‑1 (MQ=255)
                                                               aTCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTccacc                                                                                                            >  4:185918/1‑127 (MQ=255)
                                                                      gCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTgg                                                                                <  3:232509/148‑1 (MQ=255)
                                                                          ccaGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGaa                                                                            >  6:5522/1‑148 (MQ=255)
                                                                            aGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAacac                                                                         <  1:298968/149‑1 (MQ=255)
                                                                            aGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAacac                                                                         >  1:195533/1‑149 (MQ=255)
                                                                            aGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAacac                                                                         >  7:228572/1‑149 (MQ=255)
                                                                             gCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAacac                                                                         >  2:268883/1‑148 (MQ=255)
                                                                               aCTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAg                                                                      <  2:195533/149‑1 (MQ=255)
                                                                               aCTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAg                                                                      <  8:228572/149‑1 (MQ=255)
                                                                                cTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGc                                                                     <  7:86885/149‑1 (MQ=255)
                                                                                 ttGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGAATGCCGACGGCATCCACggcga                                                                                                                                >  2:192341/1‑85 (MQ=38)
                                                                                 ttGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGAATGCCGACGGCATCCACggcga                                                                                                                                <  1:192341/89‑5 (MQ=38)
                                                                                   gCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCa                                                                    <  7:7053/147‑1 (MQ=255)
                                                                                     aCCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGAc                                                                                                 >  5:153290/1‑116 (MQ=255)
                                                                                     aCCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGAc                                                                                                 <  6:153290/116‑1 (MQ=255)
                                                                                         cACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCAc                                                                   >  8:208442/1‑142 (MQ=255)
                                                                                         cACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCAc                                                                   <  7:208442/142‑1 (MQ=255)
                                                                                                 ttGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCgctg                                                    >  1:315/1‑149 (MQ=255)
                                                                                                   gAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGaa                                                                            <  4:182210/123‑1 (MQ=255)
                                                                                                   gAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGaa                                                                            >  3:182210/1‑123 (MQ=255)
                                                                                                       cTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGc                                                               >  7:97566/1‑132 (MQ=255)
                                                                                                       cTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGc                                                               <  8:97566/132‑1 (MQ=255)
                                                                                                       cTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTgcg                                              >  2:218711/1‑149 (MQ=255)
                                                                                                          tCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCgcag                                           <  6:251621/149‑1 (MQ=255)
                                                                                                            cAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCgcagc                                          <  2:315/148‑1 (MQ=255)
                                                                                                                       ggCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAAc                              <  6:138671/149‑1 (MQ=255)
                                                                                                                        gCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGaa                               >  2:165654/1‑147 (MQ=255)
                                                                                                                          gggTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCGCGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAgg                           <  7:275944/149‑1 (MQ=255)
                                                                                                                                 cTGTGGGCGCTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGAAGCCATTGCGGGTCTTGGCGCACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGAGGCAGATCAGGTCAGACt                    >  6:190849/1‑149 (MQ=255)
                                                                                                                                             aCCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTg         <  1:165654/148‑1 (MQ=255)
                                                                                                                                              ccGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTc       <  7:47439/149‑1 (MQ=255)
                                                                                                                                               cGTGAATGCCGACGGCATCCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCa                                  >  2:102757/1‑121 (MQ=255)
                                                                                                                                               cGTGAATGCCGACGGCATCCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCa                                  <  1:102757/121‑1 (MQ=255)
                                                                                                                                                  ggATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCca     >  4:258301/1‑147 (MQ=255)
                                                                                                                                                   gATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCcaca   <  5:5522/148‑1 (MQ=255)
                                                                                                                                                   gATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAg  >  6:279296/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GATCGGCAGGTCGAGGTTGACGACCAGCGGCAAGTCATCGATATCCAGGCCACGGGCCGCCACATCGGTCGCCACCAGCACTTGCACCTCACGGGCCTTGAAGCTGTCCAGCGCACGCTGGCGGGTTGCCTGTGGGCGGTCACCGTGGATGCCGTCGGCATTCACGCCTTCGGCCAGCAAGCGCTCCACCAGCTGATCGACGCCATTGCGGGTCTTGGCGAACACCAGCACCTGCTTCCAGCGCTGCTTGCGCAGCAGGTGGCAGAACAGGTCAGCCTTGCGCTTCTTGTCCACAG  >  NC_002947/5427134‑5427429

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