Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,610,648 G→C 5.7% M21I (ATG→ATC emrE → small multidrug resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,610,6480GC 5.7% 299.6 / 4.4 94M21I (ATG→ATCemrEsmall multidrug resistance protein
Reads supporting (aligned to +/- strand):  ref base G (43/45);  new base C (2/3);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

AACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTG  >  NC_002947/5610501‑5610790
                                                                                                                                                   |                                                                                                                                              
aaCCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATg                                                                                                                                                >  5:145558/1‑148 (MQ=255)
aaCCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGa                                                                                                                                               <  8:244279/149‑1 (MQ=255)
 aCCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATg                                                                                                                                                >  3:195512/1‑147 (MQ=255)
        tcCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAg                                                                                                                                             <  7:64461/143‑1 (MQ=255)
        tcCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAg                                                                                                                                             >  8:64461/1‑143 (MQ=255)
        tcCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGACGAAGTCATCGCCACAGCCTCGATGAAGGCgat                                                                                                                                        >  6:369786/1‑146 (MQ=255)
        tcCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAAAGGAGGATCCTGCCATGAATGCCTAAACCTACCTAGCCATTGCGATCTGCGCCGAACTCATCGCCACCGCCTCGATGAAGGCGGt                                                                                                                                        >  6:392189/1‑148 (MQ=255)
          ccATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTTCATCGAGGCGGTgg                                                                                                                                      <  3:133867/148‑2 (MQ=255)
          ccATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCgg                                                                                                                                         <  6:145558/145‑1 (MQ=255)
           cATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTTCATCGAGGCGGTggcg                                                                                                                                    <  2:67524/149‑4 (MQ=255)
           cATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGa                                                                                                                                      >  6:64883/1‑147 (MQ=255)
            atatCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAg                                                                                                                                   >  8:311972/1‑149 (MQ=255)
               tCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGcc                                                                                                                                >  7:251533/1‑149 (MQ=255)
               tCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGc                                                                                                                                 <  2:214374/148‑1 (MQ=255)
                 gATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCt                                                                                                                               <  8:141989/148‑1 (MQ=255)
                   tAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGa                                                                                                                             <  2:52354/148‑1 (MQ=255)
                   tAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAg                                                                                                                            <  2:341765/149‑1 (MQ=255)
                       tAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACg                                                                                                                        >  3:215258/1‑149 (MQ=255)
                            aCTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGtt                                                                                                                   <  5:64883/149‑1 (MQ=255)
                             cTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTg                                                                                                                  >  8:294003/1‑149 (MQ=255)
                              tttACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGGAGGCGGTGAAAGGCCTGAg                                                                                                                            <  1:388446/138‑1 (MQ=255)
                              tttACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGGAGGCGGTGAAAGGCCTGAg                                                                                                                            >  2:388446/1‑138 (MQ=255)
                                  ccccATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCAc                                                                                                              <  7:356680/148‑1 (MQ=255)
                                   cccATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAg                                                                                                                                   >  1:372037/1‑126 (MQ=255)
                                   cccATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAg                                                                                                                                   <  2:372037/126‑1 (MQ=255)
                                   cccATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCActg                                                                                                            >  2:357945/1‑149 (MQ=255)
                                     cATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCActgca                                                                                                          >  1:348431/1‑148 (MQ=255)
                                      aTCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCActgctg                                                                                                         >  8:65615/1‑149 (MQ=255)
                                          agtcAACGTCTATAATCCCCGCTAAGTACCGGTGCTTCCTGCAATTAATGCCTATACATACCTCGCCATTGCAATCTGGGCCTAAGTCATCGCCACCGCCTCGATCACGGGGGTTAAAGGCCTTAGCACCCCGTTGCCACTGGTGGTGa                                                                                                     <  2:389240/146‑1 (MQ=255)
                                           ggCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAt                                                                                                    <  3:124990/149‑1 (MQ=255)
                                                  tCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTTTACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGtt                                                                                             >  1:288761/1‑149 (MQ=255)
                                                  tCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGtt                                                                                             <  5:369786/149‑1 (MQ=255)
                                                   cTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTg                                                                                            >  2:84531/1‑149 (MQ=255)
                                                    tAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTgg                                                                                           >  6:130153/1‑149 (MQ=255)
                                                    tAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTgg                                                                                           >  6:301281/1‑149 (MQ=255)
                                                     aGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGc                                                                                          >  1:99235/1‑149 (MQ=255)
                                                            ccGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGc                                                                                                                                 <  8:10446/103‑1 (MQ=255)
                                                            ccGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGc                                                                                                                                 >  7:10446/1‑103 (MQ=255)
                                                            ccGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCa                                                                                   >  6:331435/1‑149 (MQ=255)
                                                              gCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCAtt                                                                                 <  3:311627/149‑1 (MQ=255)
                                                                      aCGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGCCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTg                                                                         <  5:349039/149‑1 (MQ=255)
                                                                       cGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTgg                                                                        <  4:215258/149‑1 (MQ=255)
                                                                        ggaggaTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGtt                                                                                                                   >  2:260311/1‑105 (MQ=255)
                                                                        ggaggaTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGtt                                                                                                                   <  1:260311/105‑1 (MQ=255)
                                                                         gaggaTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAtg                                                                                                   <  6:268900/120‑1 (MQ=255)
                                                                         gaggaTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAtg                                                                                                   >  5:268900/1‑120 (MQ=255)
                                                                         gaggaTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGAt                                                                      <  1:357945/149‑1 (MQ=255)
                                                                          aggaTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATg                                                                     <  2:191965/149‑1 (MQ=255)
                                                                           ggaTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGc                                                                    <  4:383155/149‑1 (MQ=255)
                                                                              tCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGa                                                                 <  7:311972/149‑1 (MQ=255)
                                                                                cTGCCATGAATGGCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGCTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTCATGGTGGTTGGCTATGGCATTGCGTACTGGATGCTGACg                                                               >  1:171047/1‑149 (MQ=255)
                                                                                 tGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGACGTAGTCATCGCCACCGCCTCGATGAAGGGGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCt                                                                          <  1:119151/137‑1 (MQ=255)
                                                                                 tGCCATGAATGCCTATACCTACATCGCCATTCCCAACTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGAGGAGAAAGGCATGAGCACGCCGAAGCCACTGCTGCTGATGGTGCTTGGCTATTGCATTGCGTTCt                                                                          >  2:119151/1‑137 (MQ=255)
                                                                                   ccATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTa                                                                                        >  2:180690/1‑121 (MQ=255)
                                                                                   ccATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTa                                                                                        <  1:180690/121‑1 (MQ=255)
                                                                                   ccATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtg                                                            <  6:226243/149‑1 (MQ=255)
                                                                                    cATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtgg                                                           <  2:348431/149‑1 (MQ=255)
                                                                                     aTGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtggt                                                          >  7:358772/1‑149 (MQ=255)
                                                                                      tGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCg                                                                              <  1:252004/128‑1 (MQ=255)
                                                                                      tGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCg                                                                              >  2:252004/1‑128 (MQ=255)
                                                                                      tGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtggtg                                                         <  5:130153/149‑1 (MQ=255)
                                                                                             cTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCg                                                    >  6:288263/1‑147 (MQ=255)
                                                                                               ataCCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGATGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgcag                                                 >  5:103577/1‑148 (MQ=255)
                                                                                               ataCCTACCTCGCCATTGCCATCTGCACCGAAGTCATAGCCACCGCCTCGAAAAAGGCGGTGAAAGGCCTGAGCACGCACTAAAAACAGCTGCTGAAGGAGCTTGGCTATGGCATTGCGTACTGGATGCTGACGCtggtg                                                         >  7:60738/1‑140 (MQ=255)
                                                                                                  cctacctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGACTATg                                                                                      >  1:168991/1‑108 (MQ=255)
                                                                                                  cctacctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGACTATg                                                                                      <  2:168991/108‑1 (MQ=255)
                                                                                                   ctacctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTAt                                                                                       <  3:362153/106‑1 (MQ=255)
                                                                                                   ctacctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTAt                                                                                       >  4:362153/1‑106 (MQ=255)
                                                                                                   ctacctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTc                                            <  8:358772/149‑1 (MQ=255)
                                                                                                   ctacctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTc                                            >  6:352859/1‑149 (MQ=255)
                                                                                                      cctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCg                                          <  5:301281/148‑1 (MQ=255)
                                                                                                      cctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCg                                          <  2:77716/148‑1 (MQ=255)
                                                                                                        tCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTg                                       >  4:180887/1‑149 (MQ=255)
                                                                                                         cGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTgg                                      >  5:370873/1‑149 (MQ=255)
                                                                                                          gCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTggg                                     <  5:331435/149‑1 (MQ=255)
                                                                                                            cATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGa                                   >  4:140501/1‑149 (MQ=255)
                                                                                                                 ccATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGcc                              >  8:242943/1‑149 (MQ=255)
                                                                                                                 ccATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCATTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGc                               >  5:387532/1‑148 (MQ=255)
                                                                                                                   aTCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTa                            >  6:377972/1‑149 (MQ=255)
                                                                                                                   aTCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTa                            >  1:8548/1‑149 (MQ=255)
                                                                                                                    tCTGCGCCGAAGTCGTCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTAc                           >  7:345577/1‑149 (MQ=255)
                                                                                                                        cgcCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTAc                           <  5:377972/145‑1 (MQ=255)
                                                                                                                        cgcCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGcc                        >  8:237813/1‑148 (MQ=255)
                                                                                                                        cgcCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGcc                        <  2:99235/148‑1 (MQ=255)
                                                                                                                        cgcCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGc                         <  7:294003/147‑1 (MQ=255)
                                                                                                                         gcCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCAt                      >  5:84942/1‑149 (MQ=255)
                                                                                                                           cGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCt                    <  7:65615/149‑1 (MQ=255)
                                                                                                                           cGAAGGCCTCGCCACCGCCTCGATGACGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCt                    <  5:392189/149‑1 (MQ=255)
                                                                                                                                   tCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCgggg            <  2:8548/149‑1 (MQ=255)
                                                                                                                                      ccACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTg         <  5:352859/149‑1 (MQ=255)
                                                                                                                                        aCCGCCTTCATCGAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTggg       >  5:187320/1‑149 (MQ=255)
                                                                                                                                        aCCGCCTTCATCGAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCGGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGAAGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTAGGGGCTggg       >  5:219357/1‑149 (MQ=255)
                                                                                                                                        aCCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCg                                                    <  5:121037/104‑1 (MQ=255)
                                                                                                                                        aCCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCg                                                    >  6:121037/1‑104 (MQ=255)
                                                                                                                                        aCCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTggg       <  1:84531/149‑1 (MQ=255)
                                                                                                                                             cTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCATTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATtg  <  6:387532/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
AACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTG  >  NC_002947/5610501‑5610790

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: