Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,428,065 A→T 7.9% L13Q (CTG→CAG)  rhlE‑I ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,428,0650AT 7.9% 183.9 / 5.2 63L13Q (CTG→CAG) rhlE‑IATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base A (30/28);  new base T (3/2);  total (33/30)
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.61e-01

GCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACT  >  NC_002947/5427919‑5428213
                                                                                                                                                   |                                                                                                                                                    
gcTGCAGCACAGGCAGGGCAAAGACTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGAGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCgcag                                                                                                                                                     >  6:64662/1‑148 (MQ=255)
     agcaCAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCgg                                                                                                                                               >  5:331086/1‑149 (MQ=255)
        acaGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTc                                                                                                                                            <  2:204117/149‑1 (MQ=255)
            gcaggGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGAt                                                                                                                                         >  7:26934/1‑148 (MQ=255)
                   aaaGCCTGGGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGTGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGa                                                                                                                                 <  6:76553/149‑1 (MQ=255)
                      gCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGtt                                                                                                                              >  5:2440/1‑149 (MQ=255)
                        cTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGTCCATCAGGTGGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGATCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTg                                                                                                                            <  7:220133/149‑1 (MQ=255)
                        cTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTg                                                                                                                            <  5:269324/149‑1 (MQ=255)
                         tGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCTGGTCGCGGGCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTgg                                                                                                                           <  5:165719/149‑1 (MQ=255)
                              tCTTCCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGaaa                                                                                                                      <  5:119023/149‑1 (MQ=255)
                                ttGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAAtt                                                                                                                    >  2:251495/1‑149 (MQ=255)
                                ttGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAAtt                                                                                                                    >  6:362265/1‑149 (MQ=255)
                                ttGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAAtt                                                                                                                    >  3:361448/1‑149 (MQ=255)
                                  gCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGGGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGACTTTGGCGAAAt                                                                                                                     <  8:286972/146‑1 (MQ=255)
                                   ccGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCAt                                                                                                                 <  4:348704/149‑1 (MQ=255)
                                     ggTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATgg                                                                                                               >  8:299111/1‑149 (MQ=255)
                                          cGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATa                                                                                                          >  4:352062/1‑149 (MQ=255)
                                                     gCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAggggg                                                                                               >  6:6334/1‑149 (MQ=255)
                                                          cATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAg                                                                                          <  5:6334/149‑1 (MQ=255)
                                                              aGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGaa                                                                                      <  7:299111/149‑1 (MQ=255)
                                                               ggTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAg                                                                                     <  4:74684/149‑1 (MQ=255)
                                                                    cggccggcCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAggcggc                                                                                <  4:333155/149‑1 (MQ=255)
                                                                       ccggccAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAg                                                                             <  2:61449/149‑1 (MQ=255)
                                                                        cggccAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGt                                                                            >  3:198222/1‑149 (MQ=255)
                                                                           ccAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTa                                                                         >  4:348040/1‑149 (MQ=255)
                                                                            cAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTAc                                                                        >  8:127230/1‑149 (MQ=255)
                                                                            cAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTAc                                                                        >  8:122000/1‑149 (MQ=255)
                                                                             aGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTTGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTAcc                                                                       <  7:146630/149‑1 (MQ=255)
                                                                             aGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTAcc                                                                       <  6:2440/149‑1 (MQ=255)
                                                                              gCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCg                                                                      >  4:159131/1‑149 (MQ=255)
                                                                                    cggGGATGGCTTTGGCCTGCACCGGTGTCGTGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGc                                                                <  7:275103/149‑1 (MQ=255)
                                                                                      gggATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGccc                                                              <  3:352062/149‑1 (MQ=255)
                                                                                           ggCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGcc                                                         >  6:341619/1‑149 (MQ=255)
                                                                                                    cTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTg                                                   >  5:119234/1‑146 (MQ=255)
                                                                                                       cACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGa                                               <  6:119234/147‑1 (MQ=255)
                                                                                                        aCCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGAtt                                             >  8:276195/1‑148 (MQ=255)
                                                                                                         ccGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGAttt                                            >  7:212415/1‑148 (MQ=255)
                                                                                                          cggggtcggggtGGTGTAGTCCAGTTGCATGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGAttt                                            >  1:23342/1‑148 (MQ=255)
                                                                                                          cggggtcggggtGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCa                                          >  1:260069/1‑149 (MQ=255)
                                                                                                          cggggtcggggtGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCa                                          <  8:26934/149‑1 (MQ=255)
                                                                                                          cggggtcggggtGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCa                                          <  5:362265/149‑1 (MQ=255)
                                                                                                           ggggtcggggtGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCAc                                         <  1:251495/149‑1 (MQ=255)
                                                                                                           ggggtcggggtGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCAc                                         <  6:331086/149‑1 (MQ=255)
                                                                                                                         gtAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTg                            >  6:170616/1‑148 (MQ=255)
                                                                                                                          tAGTCCAGTTGC‑TGCAGCGTGCGCAGCGGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTa                          >  7:90345/1‑149 (MQ=255)
                                                                                                                          tAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCa                                                                                                                  >  6:241779/1‑61 (MQ=255)
                                                                                                                          tAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCa                                                                                                                  <  5:241779/61‑1 (MQ=255)
                                                                                                                            gtacAGTTTC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGGCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTAc                         <  5:64662/145‑1 (MQ=255)
                                                                                                                             tCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGTCCCTGCCTACTTGCAGATTTCAc                                         <  5:259402/131‑1 (MQ=255)
                                                                                                                             tCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGTCCCTGCCTACTTGCAGATTTCAc                                         >  6:259402/1‑131 (MQ=255)
                                                                                                                             tCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAg                       >  3:308762/1‑149 (MQ=255)
                                                                                                                             tCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAg                       >  5:171063/1‑149 (MQ=255)
                                                                                                                              ccAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAgg                      >  6:45783/1‑149 (MQ=255)
                                                                                                                               cAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAgg                      <  5:341619/148‑1 (MQ=255)
                                                                                                                                 gTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGtt                    >  3:67292/1‑148 (MQ=255)
                                                                                                                                 gTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTc                   <  7:325347/149‑1 (MQ=255)
                                                                                                                                   cgc‑tgctGCGCACGCTGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAg                 >  3:313381/2‑149 (MQ=255)
                                                                                                                                   cgc‑tgctGCGCACGCTGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAg                 >  3:318622/2‑149 (MQ=255)
                                                                                                                                       tgcAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGCTTTCACCTTTCAACGGTTGTACAGGTTCAGgcg              >  5:279696/1‑149 (MQ=255)
                                                                                                                                             gcacgcTGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTgc        <  4:318622/146‑1 (MQ=255)
                                                                                                                                             gcacgcTGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTgc        <  4:313381/146‑1 (MQ=255)
                                                                                                                                              tGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGtt                             <  3:363315/127‑1 (MQ=255)
                                                                                                                                              tGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGtt                             >  4:363315/1‑127 (MQ=255)
                                                                                                                                              cacgcTGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAgg                                                                                                  >  4:97928/3‑58 (MQ=255)
                                                                                                                                                 ccagcagCGGTTCGAACAGGGCGAGTTTGGCGAAATTCAAGGCAATTCCGTCAGGGGGTTCAGCGAAGCCGGCAAGTTTTCCGCATCGCCCCTGACTACTTAAAGATTTCACCTTTCAACCGTGGTACAGGTACAGGCGCCTAGCGCCa    >  4:306219/2‑149 (MQ=255)
                                                                                                                                                   agcagcGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACt  >  3:123875/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGC‑TGCAGCGTGCGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACT  >  NC_002947/5427919‑5428213

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