Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,856,520 T→A 6.3% F149L (TTT→TTA)  PP_5133 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,856,5200TA 6.3% 184.4 / 3.0 63F149L (TTT→TTA) PP_5133membrane protein
Reads supporting (aligned to +/- strand):  ref base T (33/26);  new base A (2/2);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

CAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCTGGGCATGCTACTGGTCTTG  >  NC_002947/5856374‑5856667
                                                                                                                                                  |                                                                                                                                                   
cAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCAAAGCGTTCGCTTAGc                                                                                                                                                   <  2:254686/149‑1 (MQ=255)
     cccGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCACTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCcgc                                                                                                                                               >  2:238607/1‑148 (MQ=255)
      ccGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGa                                                                                                                                             >  4:350359/1‑149 (MQ=255)
          tGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGcc                                                                                                                                         <  1:382285/149‑1 (MQ=255)
          tGAAGGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGACTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGcc                                                                                                                                         <  4:325853/149‑1 (MQ=255)
           gAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGcc                                                                                                                                         >  7:52064/1‑148 (MQ=255)
              tgtgTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCAGGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAgg                                                                                                                                     <  2:255522/149‑1 (MQ=255)
                 gtTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGcc                                                                                                                                                  >  3:70355/1‑133 (MQ=255)
                 gtTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGcc                                                                                                                                                  <  4:70355/133‑1 (MQ=255)
                    cTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCcgcg                                                                                                                                >  7:53199/1‑148 (MQ=255)
                     tGGGCGCTGGGCAGCCTGCTCGGGGTGAACCTGCTGATGCTTCTTGGCTTGGGCTTGGGTTTTGCCCTGAGCGGCGAGCATGGCGGAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCa                                                                                                                              <  4:69987/149‑1 (MQ=255)
                       ggCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCAcc                                                                                                                            >  6:342797/1‑149 (MQ=255)
                        gcgcTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCt                                                                                                                           >  1:326500/1‑149 (MQ=255)
                                cAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAg                                                                                                                   <  4:5293/149‑1 (MQ=255)
                                      gCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCgctgc                                                                                                             <  8:52064/149‑1 (MQ=255)
                                          gggCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGt                                                                                                         >  7:95588/1‑149 (MQ=255)
                                                   ctgctgATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGtgctgctgc                                                                                                 >  4:211697/1‑148 (MQ=255)
                                                     gctgATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAAc                                                                                              >  3:381002/1‑149 (MQ=255)
                                                      ctgATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGt                                                                                                         >  2:349795/1‑137 (MQ=255)
                                                        gATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAAGGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGgcgc                                                                                                                      >  2:93860/1‑122 (MQ=255)
                                                        gATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAAGGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGgcgc                                                                                                                      <  1:93860/122‑1 (MQ=255)
                                                           gCTACATGGCTGGGGCATGGGTATTGCCCTGACCGGCGAGCATGGCGCAGCGCTGGGCCTCTAGTGGCTGTGGCTAAGCGAACGCTATGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTCCTGGTGCTGCTGCAACGTCAAc                                                                                        >  6:114209/1‑149 (MQ=255)
                                                              aCTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGc                                                                                      >  8:115769/1‑148 (MQ=255)
                                                               cTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGcc                                                                                     <  8:95588/148‑1 (MQ=255)
                                                                   gCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAAc                                                                                <  3:350359/149‑1 (MQ=255)
                                                                           aTGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGgct                                                                        >  4:42387/1‑149 (MQ=255)
                                                                              ggTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCCACCGGTGgctgct                                                                     >  5:322269/1‑149 (MQ=255)
                                                                                   tGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAATGTCAACGCCTCAACCGGTGGCTGCTCGGc                                                                 <  1:237043/148‑1 (MQ=255)
                                                                                      ccTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCa                                                                                                                    >  1:333750/1‑94 (MQ=255)
                                                                                      ccTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCa                                                                                                                    <  2:333750/94‑1 (MQ=255)
                                                                                      ccTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGt                                                             >  1:16004/1‑149 (MQ=255)
                                                                                       cTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGtt                                                            >  5:33027/1‑149 (MQ=255)
                                                                                           gcggcgAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTa                                                        <  6:33027/149‑1 (MQ=255)
                                                                                            cggggAGCATGGGGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCTAGGCGGCGCACCTGGCGCCAGCGCTGGTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTac                                                       <  5:114209/149‑1 (MQ=255)
                                                                                            cggcgAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTac                                                       <  7:115769/149‑1 (MQ=255)
                                                                                             ggcgAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTaca                                                      >  6:200483/1‑149 (MQ=255)
                                                                                               cgAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACg                                                    <  4:381002/149‑1 (MQ=255)
                                                                                                     tGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCtgt                                              <  8:262271/149‑1 (MQ=255)
                                                                                                       gcgcAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCtgtg                                             <  2:130723/148‑1 (MQ=255)
                                                                                                        cgcAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACgg                                                   <  7:286037/141‑1 (MQ=255)
                                                                                                        cgcAGGGCTGGGCCGCTTGTGGCAGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACgg                                                   >  8:286037/1‑141 (MQ=255)
                                                                                                             ggCTGGGCCGCTTGTAGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGctgctg                                      >  7:241391/1‑149 (MQ=255)
                                                                                                              gCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTgg                                     <  2:326500/149‑1 (MQ=255)
                                                                                                                tGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGctgctg                                      >  6:311478/1‑146 (MQ=255)
                                                                                                                tGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGctgctg                                      <  5:311478/146‑1 (MQ=255)
                                                                                                                tGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGc                                    >  5:23004/1‑148 (MQ=255)
                                                                                                                   gCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATg                                <  3:34384/149‑1 (MQ=255)
                                                                                                                   gCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATg                                <  3:42387/149‑1 (MQ=255)
                                                                                                                   gCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGATGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGATGGTACACGGCCTGTGGCTGCTGGCGAt                                 >  3:396895/1‑148 (MQ=255)
                                                                                                                    ccGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGc                               >  2:277795/1‑149 (MQ=255)
                                                                                                                        ttGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCAc                           >  5:113851/1‑149 (MQ=255)
                                                                                                                             tcgcgggcaAAGCGTTCGCTTAGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCgcgc                      >  6:13293/6‑149 (MQ=255)
                                                                                                                             gctgtggctAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCgcgc                      >  3:102156/1‑149 (MQ=255)
                                                                                                                              ctgtggctAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAGCGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCt                     >  3:43190/1‑149 (MQ=255)
                                                                                                                                cgggcaAAGCGTTCGCTTAGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACACCTCAACCGGTGGCTGCTCGGCCTGTTGGTa                                                        >  6:276773/3‑112 (MQ=255)
                                                                                                                                cgggcaAAGCGTTCGCTTAGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACACCTCAACCGGTGGCTGCTCGGCCTGTTGGTa                                                        <  5:276773/110‑1 (MQ=255)
                                                                                                                                   gctAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCTGGGc                 <  6:113851/148‑1 (MQ=255)
                                                                                                                                     tAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCTGGGCAt               >  5:299273/1‑148 (MQ=255)
                                                                                                                                       aGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCTGGGCAt               <  6:299273/146‑1 (MQ=255)
                                                                                                                                       aGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCTGGGCATGc             >  6:233722/1‑148 (MQ=255)
                                                                                                                                        gCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAGCGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCTGGGCATGCTa           <  4:43190/149‑1 (MQ=255)
                                                                                                                                        gCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACgg                                                   >  8:324388/1‑109 (MQ=255)
                                                                                                                                              gCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCTGGGCATGCTACTGGTc     >  7:9738/1‑149 (MQ=255)
                                                                                                                                                  tGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCTGGGCATGCTACTGGTCTTg  <  8:241391/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTTGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTCGGCCTGTTGGTACACGGCCTGTGGCTGCTGGCGATGCTCACCGCGCTGGGCATGCTACTGGTCTTG  >  NC_002947/5856374‑5856667

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