Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,585,219 G→A 5.6% intergenic (‑206/‑170) PP_1388 ← / → PP_1389 EmrB/QacA family drug resistance transporter/oxaloacetate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,585,2190GA 5.6% 212.8 / 2.1 71intergenic (‑206/‑170)PP_1388/PP_1389EmrB/QacA family drug resistance transporter/oxaloacetate decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base A (2/2);  total (36/35)
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.45e-01

CCTCGGATTGTAACCTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAGGGGTCGGTAGT  >  NC_002947/1585073‑1585359
                                                                                                                                                  |                                                                                                                                            
ccTCGGATTGTAACCTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGt                                                                                                                                            >  6:138950/1‑149 (MQ=255)
 cTCGGATTGTAACCTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGtt                                                                                                                                           >  6:248955/1‑149 (MQ=255)
 cTCGGATTGTAACCTGAGTTAGGTAGCTTCCTATGTACTATCTGCACCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGtt                                                                                                                                           >  1:272807/1‑149 (MQ=255)
    ggATTGTAACCTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGt                                                                                                                                        <  6:7957/149‑1 (MQ=255)
          tAACCTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATGCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTgcagc                                                                                                                                   <  5:341500/148‑1 (MQ=255)
          tAACCTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAACCTGGTgcaag                                                                                                                                   <  3:41720/148‑3 (MQ=255)
           aaCCTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCACCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTgcagca                                                                                                                                  >  2:159520/1‑148 (MQ=255)
              cTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATAtt                                                                                                                              >  7:296096/1‑149 (MQ=255)
                 aGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGcc                                                                                                                           <  5:91315/149‑1 (MQ=255)
                      ggTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACTGAGGCTGAAATCCGGGGCTGCACCAACCTGGTgcaagccc                                                                                                                                <  5:316659/138‑6 (MQ=255)
                      ggTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACTGAGGCTGAAATCCGGGGCTGCACCAACCTGGTgcaagccc                                                                                                                                >  6:316659/1‑133 (MQ=255)
                       gTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGt                                                                                                                     >  8:52293/1‑149 (MQ=255)
                           cTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGttttt                                                                                                                 <  6:138556/149‑1 (MQ=255)
                            ttCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGttttt                                                                                                                 >  8:334937/1‑148 (MQ=255)
                              ccTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCa                                                                                                               <  3:95140/148‑1 (MQ=255)
                                tATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAg                                                                                                                       >  1:309319/1‑138 (MQ=255)
                                tATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAg                                                                                                                       <  2:309319/138‑1 (MQ=255)
                                tATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTAAGTc                                                                                                            >  2:106578/1‑149 (MQ=255)
                                  tGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGCATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCg                                                                                                          >  2:5998/1‑149 (MQ=255)
                                    tACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCt                                                                                                        <  1:6387/149‑1 (MQ=255)
                                    tACTATCTGCATACCCTATGCAAGATCGTCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGc                                                                                                         <  7:329972/148‑1 (MQ=255)
                                     aCTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTg                                                                                                       >  3:106039/1‑149 (MQ=255)
                                         tCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTg                                                                                                                                       <  5:5555/113‑1 (MQ=255)
                                         tCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTg                                                                                                                                       >  6:5555/1‑113 (MQ=255)
                                         tCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGtctc                                                                                                   >  1:273418/1‑149 (MQ=255)
                                          cTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGtctc                                                                                                   <  7:261888/148‑1 (MQ=255)
                                               tCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCacac                                                                                             <  7:334937/149‑1 (MQ=255)
                                                  ccTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGttttt                                                                                                                 >  5:96724/1‑126 (MQ=255)
                                                  ccTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGttttt                                                                                                                 <  6:96724/126‑1 (MQ=255)
                                                   cTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGATGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCttt                                                                                         >  7:284090/1‑149 (MQ=255)
                                                      tGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTAt                                                                                      <  5:72686/149‑1 (MQ=255)
                                                        cAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTAtt                                                                                     <  1:150980/148‑1 (MQ=255)
                                                           gATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCctc                                                                                 <  8:296096/149‑1 (MQ=255)
                                                            aTCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCctct                                                                                >  1:36588/1‑149 (MQ=255)
                                                                gCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTctgcct                                                                            <  5:248955/149‑1 (MQ=255)
                                                                       cgcATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGt                                                                      <  7:52293/148‑1 (MQ=255)
                                                                       cgcATTCATCCAGTGCAGCATGGCTGCGTGATATCGTTTGTATGACACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCTCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTc                                                                     <  8:276777/149‑1 (MQ=255)
                                                                         cATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGttttt                                                                                                                 >  8:28823/1‑103 (MQ=255)
                                                                         cATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGttttt                                                                                                                 <  7:28823/103‑1 (MQ=255)
                                                                         cATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGa                                                                   <  7:285747/149‑1 (MQ=255)
                                                                         cATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGa                                                                   <  5:138950/149‑1 (MQ=255)
                                                                              aTCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTGTCAGTCCGCTGTCTCCCACACCTTTTATTACTCTGCCTGCATGTCGAGCATg                                                              >  2:257441/1‑149 (MQ=255)
                                                                               tCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCCGTCCGCTGGCTCCCACACGTTTTATGCATCTTCCTGCATTTCGAGAATGa                                                             >  8:321226/1‑149 (MQ=255)
                                                                                  aGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCAAACATTTTATTCCTCTGCATGCAAGTCGAGCATCACCg                                                          >  4:89773/1‑149 (MQ=255)
                                                                                      ctgcATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGccc                                                       <  1:159520/148‑1 (MQ=255)
                                                                                           tGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATgg                                                 >  2:223297/1‑149 (MQ=255)
                                                                                                 cGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACgg                                           <  2:36588/149‑1 (MQ=255)
                                                                                                  gTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGt                                          <  4:106039/149‑1 (MQ=255)
                                                                                                          gTTCGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCAGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAAcag                                  <  2:289186/149‑1 (MQ=255)
                                                                                                              gTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAAcaggtc                               <  6:4592/148‑1 (MQ=255)
                                                                                                               tATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAAcaggtcag                             >  7:142430/1‑149 (MQ=255)
                                                                                                                aTGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGa                                                             <  1:184789/116‑1 (MQ=255)
                                                                                                                aTGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGa                                                             >  2:184789/1‑116 (MQ=255)
                                                                                                                 tGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAAcaggtcagg                            >  4:11250/1‑148 (MQ=255)
                                                                                                                   ccACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAAcaggtcagg                            >  4:165118/1‑146 (MQ=255)
                                                                                                                   ccACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAAcaggtcagg                            <  3:165118/146‑1 (MQ=255)
                                                                                                                     acaGAGTCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGAc                       <  6:342380/149‑1 (MQ=255)
                                                                                                                     acaGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGAc                       >  2:2365/1‑149 (MQ=255)
                                                                                                                     acaGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTACTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGAc                       >  3:56653/1‑149 (MQ=255)
                                                                                                                               aaaTCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAg             >  5:141633/1‑149 (MQ=255)
                                                                                                                                aaTCCGTGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAg             <  3:89773/148‑1 (MQ=255)
                                                                                                                                aaTCCGGGGCTTTCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATATCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAg             <  5:366142/148‑1 (MQ=255)
                                                                                                                                aaTCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATATCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAg             >  6:366142/1‑148 (MQ=255)
                                                                                                                                   ccGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAGGGGt         >  5:201933/1‑149 (MQ=255)
                                                                                                                                    cGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCtttt                                                                                        <  3:279436/69‑1 (MQ=255)
                                                                                                                                    cGGCGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCtttt                                                                                        >  4:279436/1‑69 (MQ=255)
                                                                                                                                    aatatgctGCACCAACCTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAGGGGTc        >  3:82855/8‑149 (MQ=255)
                                                                                                                                     ggggCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAGGGGTCg       >  1:26601/1‑149 (MQ=255)
                                                                                                                                      gggCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAGGGGTCgg      <  3:11250/149‑1 (MQ=255)
                                                                                                                                          ttGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAGGGGTCGGTAGt  >  3:344268/1‑149 (MQ=255)
                                                                                                                                          ttGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAGGGGTCGGTAGt  >  1:206638/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
CCTCGGATTGTAACCTGAGTTAGGTAGCTTCCTATGTACTATCTGCATCCCCTATGCAAGATCGGCATTGGCGCATTCATCCAGTGCTGCATGGCTGCGTGATATCGTTGGTATGCCACAGAGGCTGAAATCCGGGGCTTGCACCAGGTTGGTGCAGCATATTGCCGCAGGTTTTTCAGTCCGCTGTCTCCCACACCTTTTATTCCTCTGCCTGCATGTCGAGCATGACCGCCCAGATGGCGACGGTTGGAACAGGTCAGGTAGACCCGATTACAGGGGTCGGTAGT  >  NC_002947/1585073‑1585359

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: