Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,217,808 G→T 5.6% E185D (GAG→GAT oprJ → outer membrane protein OprJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,217,8080GT 5.6% 283.2 / 2.1 90E185D (GAG→GAToprJouter membrane protein OprJ
Reads supporting (aligned to +/- strand):  ref base G (38/46);  new base T (3/3);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GAATACGAAGTGGACCTGTTCGGCCGCGTCAAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAA  >  NC_002947/3217659‑3217955
                                                                                                                                                     |                                                                                                                                                   
gAATACGAAGTGGACCTGTTCGGCCGCGTCAAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAg                                                                                                                                                     >  2:554695/1‑150 (MQ=255)
     cGAAGTGGACCTGTTCGGCCGCGTCAAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAgcggcg                                                                                                                                                >  1:202047/1‑150 (MQ=255)
          tGGACCTGTTCGGCCGCGTCAAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAgg                                                                                                                                           >  2:1421311/1‑150 (MQ=255)
           ggACCTGTTCGGCCGCGTCAAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGgc                                                                                                                                          <  2:1759065/150‑1 (MQ=255)
             aCCTGTTCGGCCGCGTCAAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGCCGGCGCCGAGCGGCGCCAGgcgc                                                                                                                                        >  1:1653859/1‑150 (MQ=255)
               cTGTTCGGCCGCGTCAAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTg                                                                                                                                      >  2:296982/1‑150 (MQ=255)
                     ggCCGCGTCAAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGt                                                                                                                                <  1:747177/150‑1 (MQ=255)
                       ccGCGTCAAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCa                                                                                                                              <  1:1614191/150‑1 (MQ=255)
                         gcgTCAAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGa                                                                                                                            >  2:226473/1‑150 (MQ=255)
                          cgTCAAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGAc                                                                                                                           <  1:554695/150‑1 (MQ=255)
                           gTCAAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTCATCGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGAcc                                                                                                                          >  1:1534198/1‑150 (MQ=255)
                             cAANAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCt                                                                                                                        >  2:1187592/1‑150 (MQ=255)
                             cAAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCt                                                                                                                        >  2:380915/1‑150 (MQ=255)
                             cAAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCt                                                                                                                        <  2:236299/150‑1 (MQ=255)
                              aaaaGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTg                                                                                                                       >  2:1271959/1‑150 (MQ=255)
                              aaaaGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTg                                                                                                                       >  1:1461518/1‑150 (MQ=255)
                                   ccTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAg                                                                                                                  >  1:384692/1‑150 (MQ=255)
                                    cTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGc                                                                                                                 <  1:820185/150‑1 (MQ=255)
                                     tGAGCGATGCGGCCCTGGAGCAGTACCTGGACAACGAGCAATTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGcc                                                                                                                <  2:1399958/150‑1 (MQ=255)
                                      gAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCGAGGCCTCCCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCcg                                                                                                               >  1:699059/1‑150 (MQ=255)
                                         cGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGa                                                                                                            <  2:244045/150‑1 (MQ=255)
                                         cGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGa                                                                                                            <  1:406490/150‑1 (MQ=255)
                                          gATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAg                                                                                                           <  1:1014208/150‑1 (MQ=255)
                                            tGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGa                                                                                                         >  2:1352075/1‑150 (MQ=255)
                                            tGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGa                                                                                                         >  1:1303470/1‑150 (MQ=255)
                                            tGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGa                                                                                                         >  2:806717/1‑150 (MQ=255)
                                                 cccTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTccc                                                                                                    <  1:806717/150‑1 (MQ=255)
                                                  ccTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACAGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCt                                                                                                   <  2:985373/150‑1 (MQ=255)
                                                   cTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTg                                                                                                  >  2:373131/1‑150 (MQ=255)
                                                   cTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTg                                                                                                  >  1:120949/1‑150 (MQ=255)
                                                   cTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTg                                                                                                  <  1:373131/150‑1 (MQ=255)
                                                     ggAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGc                                                                                                <  1:1135243/150‑1 (MQ=255)
                                                     ggAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGc                                                                                                <  1:1849637/150‑1 (MQ=255)
                                                                    cACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAgcg                                                                                 >  2:1576117/1‑150 (MQ=255)
                                                                    cACCGAGGAAGTCGGGCGTGCCGCACACATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAgcg                                                                                 >  1:1696955/1‑150 (MQ=255)
                                                                     aCCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAgcgg                                                                                >  1:284228/1‑150 (MQ=255)
                                                                     aCCGAGGAAGTCGGGCGTGCCGCACGCATCGAATCGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGATCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAgcgg                                                                                >  1:1590621/1‑150 (MQ=255)
                                                                       cGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAgcggcg                                                                              >  2:828442/1‑150 (MQ=255)
                                                                       cGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAgcggcg                                                                              >  2:357137/1‑150 (MQ=255)
                                                                        gAGGAAGTCGGGCGTGCCGCACGCATCGCCCGGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTAACCGCCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGgcgc                                                                             <  1:595336/150‑1 (MQ=255)
                                                                                   gCGTGACGCATGCATCTCCCTGATCCCCGAGGTGCGCCAGGCCTACCTGGCCTTGGACGCCGCCGAGCGGCGCCATATGATGACCGGTCTGACCGTTGTCAGCCGCGAGGATTCCCTGGCGTTGTTCGGGCAGCGGCGCGCTGCCGGcac                                                                  <  2:583224/150‑1 (MQ=255)
                                                                                    cGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGcacg                                                                 <  1:781627/150‑1 (MQ=255)
                                                                                     gTGCCGCACGCAACGCCCTGATCGCCGAGGCGCGAACGGCTTGCCGCGCCGTGGACGGAACCGAGCGGCGCCAGCCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGAGTTGGTCGGCCAGCGGCGCGCCGCCGGcacgg                                                                <  2:912760/150‑1 (MQ=255)
                                                                                         cgcacgcaTCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGc                                                                                                                 <  2:1853032/97‑1 (MQ=255)
                                                                                         cgcacgcaTCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGc                                                                                                                 >  1:1853032/1‑97 (MQ=255)
                                                                                                aTCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTg                                                     >  1:1735922/1‑150 (MQ=255)
                                                                                                  cGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGTCGCGCCGCCGGCACGGCCACGGCGCTGGa                                                   <  1:815864/150‑1 (MQ=255)
                                                                                                   gCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGac                                                  <  2:1303470/150‑1 (MQ=255)
                                                                                                     ccTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGcacg                                                                 >  1:1412157/1‑133 (MQ=255)
                                                                                                     ccTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGcacg                                                                 <  2:1412157/133‑1 (MQ=255)
                                                                                                     ccTGATCGCCGAGGGGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGacca                                                <  2:692144/150‑1 (MQ=255)
                                                                                                      cTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGaccac                                               <  1:1528125/150‑1 (MQ=255)
                                                                                                      cTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGaccac                                               <  1:1352075/150‑1 (MQ=255)
                                                                                                       tGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGaccacc                                              <  2:384692/150‑1 (MQ=255)
                                                                                                       tGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGaccacc                                              >  2:94538/1‑150 (MQ=255)
                                                                                                         aTCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGTCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCaa                                            <  1:1187592/150‑1 (MQ=255)
                                                                                                         aTCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCaa                                            <  1:380915/150‑1 (MQ=255)
                                                                                                           cGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAaga                                          <  1:296982/150‑1 (MQ=255)
                                                                                                               gAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCCGACTCAAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGgcg                                      >  1:1670087/1‑150 (MQ=255)
                                                                                                                aGGTGAGCCAGGCCTACCTGGCCTTGGACGAGCCTGGTTCGCTCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGgcgc                                     <  2:1590621/150‑1 (MQ=255)
                                                                                                                   tGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTgg                                  <  2:1217313/150‑1 (MQ=255)
                                                                                                                   ggAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGTCCCTGGCCAGCCGCGGGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTgg                                  <  1:1670490/149‑1 (MQ=255)
                                                                                                                       ccAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCt                              >  1:603464/1‑150 (MQ=255)
                                                                                                                        cAGGCCTACCTGGCCTTGGACGGGGCCGAGCGGCGCCGGGCGCTGACCCGTCCGGCGCTGGCGAGCCGCGAGGGTTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTg                             <  2:1653859/150‑1 (MQ=255)
                                                                                                                        cAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTg                             <  1:226473/150‑1 (MQ=255)
                                                                                                                         aGGCCTACCTGGCCTTGGACGGCGCCGGGCCGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCTTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTgg                            <  2:1534198/150‑1 (MQ=255)
                                                                                                                         aGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTgg                            <  1:482284/150‑1 (MQ=255)
                                                                                                                         aGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTgg                            <  1:1576117/150‑1 (MQ=255)
                                                                                                                          ggCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCAGCCCGCACTGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGt                           <  2:1809620/150‑1 (MQ=255)
                                                                                                                           gCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTTGCCAGCCGCGAGGATTCCCTg                                                                                                  <  2:1656496/78‑1 (MQ=255)
                                                                                                                           gCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTg                                                                                                  >  1:1656496/1‑78 (MQ=255)
                                                                                                                           gCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTc                          >  2:164851/1‑150 (MQ=255)
                                                                                                                           gCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTc                          >  1:972861/1‑150 (MQ=255)
                                                                                                                             cTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGa                        >  1:270133/1‑150 (MQ=255)
                                                                                                                             cTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGa                        <  2:603464/150‑1 (MQ=255)
                                                                                                                              tACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGac                                                  <  2:213609/123‑1 (MQ=255)
                                                                                                                              tACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGac                                                  >  1:213609/1‑123 (MQ=255)
                                                                                                                               aCCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGc                      <  2:202047/150‑1 (MQ=255)
                                                                                                                                ccTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGcc                     >  2:550848/1‑149 (MQ=255)
                                                                                                                                ccTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCa                     <  1:164851/150‑1 (MQ=255)
                                                                                                                                ccTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCa                     <  2:972861/150‑1 (MQ=255)
                                                                                                                                      cTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCg               >  1:1142437/1‑150 (MQ=255)
                                                                                                                                       ttGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGt              >  1:1208103/1‑150 (MQ=255)
                                                                                                                                       ttGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGt              >  1:1266123/1‑150 (MQ=255)
                                                                                                                                           ctggCGCCGCTCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGcacggccac                                                            <  2:674083/98‑1 (MQ=255)
                                                                                                                                           ctggCGCCGCTCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGcacggccac                                                            <  2:675494/98‑1 (MQ=255)
                                                                                                                                           ctggCGCCGCTCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGcacggccac                                                            >  1:675494/3‑100 (MQ=255)
                                                                                                                                           ctggCGCCGCTCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGcacggccac                                                            >  1:674083/3‑100 (MQ=255)
                                                                                                                                                gcCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTc     >  2:1279932/1‑150 (MQ=255)
                                                                                                                                                 ccGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCg    >  1:245485/1‑150 (MQ=255)
                                                                                                                                                   gAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGaa  <  2:1142437/150‑1 (MQ=255)
                                                                                                                                                   gAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCACGACGACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGaa  <  2:1706691/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAATACGAAGTGGACCTGTTCGGCCGCGTCAAAAGCCTGAGCGATGCGGCCCTGGAGCAGTACCTGGCCACCGAGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAA  >  NC_002947/3217659‑3217955

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: