Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 713,098 C→T 7.2% I29I (ATC→ATT) ‡ lspA → lipoprotein signal peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947713,0980CT 7.2% 245.3 / 6.2 83I29I (ATC→ATT) ‡lspAlipoprotein signal peptidase
Reads supporting (aligned to +/- strand):  ref base C (36/41);  new base T (2/4);  total (38/45)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

CCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAGCGCTGG  >  NC_002947/712956‑713245
                                                                                                                                              |                                                                                                                                                    
ccATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAg                                                                                                                                                <  5:249091/149‑1 (MQ=255)
   tCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGc                                                                                                                                             <  3:285775/148‑1 (MQ=255)
   tCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGc                                                                                                                                             >  4:285775/1‑148 (MQ=255)
     cGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGccac                                                                                                                                          >  1:271979/1‑149 (MQ=255)
       gAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGGTGGCCTGGTCGATGACCAggaccagca                                                                                                                                        <  7:324978/149‑8 (MQ=255)
       gAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCGATGACCAggaccagc                                                                                                                                         <  8:82412/149‑7 (MQ=255)
           aCCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCt                                                                                                                                    <  6:317161/149‑1 (MQ=255)
           aCCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCt                                                                                                                                    <  8:239659/149‑1 (MQ=255)
             ccGGAAATCTGCGGCCGTTGCGTCGCCAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGCCCAGGCCCCCAAGCTg                                                                                                                                   >  7:43956/1‑148 (MQ=255)
             ccGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGGTGGCCTGGTCGATGACCAGGACCAGCAAGCtca                                                                                                                                  <  2:193836/149‑3 (MQ=38)
             ccGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGCGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTg                                                                                                                                   >  2:146847/1‑148 (MQ=255)
                tcttgatGCGGCAGTTGCGTAGACATCCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGCGCGCTTCGGGCGCCTTGCATGGCTGTGGGTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGt                                                                                                                                  <  7:38554/140‑1 (MQ=255)
                 aaaTCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACtt                                                                                                                              <  6:203976/149‑1 (MQ=255)
                  aaTCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACtt                                                                                                                              >  1:98094/1‑148 (MQ=255)
                      tGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTc                                                                                                                             <  1:26245/145‑1 (MQ=255)
                      tGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTc                                                                                                                             >  2:26245/1‑145 (MQ=255)
                        cGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTcaa                                                                                                                           <  2:86907/145‑1 (MQ=255)
                         ggCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACg                                                                                                                      <  6:142779/149‑1 (MQ=255)
                          gCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGc                                                                                                                     >  3:119541/1‑149 (MQ=255)
                            cGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAAGGCTGCAGGGGGGCTCTTCGGGCGGCTTGCTTGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGAGGAGGCCACCAAGGTGTACTTCTAGAACGcgg                                                                                                                   <  6:132400/149‑3 (MQ=255)
                                gCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGa                                                                                                                <  1:235645/148‑1 (MQ=255)
                                 cGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTcaac                                                                                                                          >  6:275563/1‑137 (MQ=255)
                                 cGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTcaac                                                                                                                          <  5:275563/137‑1 (MQ=255)
                                 cGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGAc                                                                                                               <  2:197025/148‑1 (MQ=255)
                                     gACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTg                                                                                                                                   >  5:340129/1‑124 (MQ=255)
                                     gACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTg                                                                                                                                   <  6:340129/124‑1 (MQ=255)
                                     gACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATg                                                                                                           <  5:183658/148‑1 (MQ=255)
                                     gACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATg                                                                                                           >  2:190885/1‑148 (MQ=255)
                                       cAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTAc                                                                                                        <  4:285690/149‑1 (MQ=255)
                                       cAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTAc                                                                                                        >  3:285690/1‑149 (MQ=255)
                                        aaCCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTAcc                                                                                                       <  7:298124/149‑1 (MQ=255)
                                        aaCCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTAcc                                                                                                       <  2:98094/149‑1 (MQ=255)
                                        aaCCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTAcc                                                                                                       <  3:117957/149‑1 (MQ=255)
                                        aaCCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTAc                                                                                                        >  7:88878/1‑148 (MQ=255)
                                          ccTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACcag                                                                                                     <  8:326316/149‑1 (MQ=255)
                                            tGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACcag                                                                                                     >  8:2877/1‑147 (MQ=255)
                                              aCCGGCGCGGGGGCCGTGCGCCACTCTGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCGGTGGCTGAGCTTGCTGGTCCTGGTCCTCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGa                                                                                                 <  2:49299/149‑1 (MQ=255)
                                                        gTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCAt                                                                                       >  1:115043/1‑149 (MQ=255)
                                                                 gcCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGacta                                                                              <  5:86020/149‑1 (MQ=255)
                                                                   cACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCaacaac                                                                                                                       >  4:138706/1‑106 (MQ=255)
                                                                   cACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCaacaac                                                                                                                       <  3:138706/106‑1 (MQ=255)
                                                                        tGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAg                                                                        <  7:2877/148‑1 (MQ=255)
                                                                           cTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAATCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTgg                                                                    <  5:152584/149‑1 (MQ=255)
                                                                           cTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGa                                                                                 >  3:328786/1‑136 (MQ=255)
                                                                           cTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGa                                                                                 <  4:328786/136‑1 (MQ=255)
                                                                           cTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTg                                                                     >  3:313142/1‑148 (MQ=255)
                                                                              aCCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGAc                                                                  <  8:88878/148‑1 (MQ=255)
                                                                               cccTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCAGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGAccc                                                                >  6:139499/1‑149 (MQ=255)
                                                                                  tGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAGGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGCCTACTTCAGCTGGACCCTgg                                                             >  2:93477/1‑149 (MQ=255)
                                                                                  tGCAGAGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTgg                                                             >  3:97449/1‑149 (MQ=255)
                                                                                    cAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCt                                                           >  8:108707/1‑149 (MQ=255)
                                                                                       cGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCCTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTac                                                        <  7:290983/149‑1 (MQ=255)
                                                                                        ggggCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTac                                                        >  4:271338/1‑148 (MQ=255)
                                                                                          ggCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTacaa                                                      >  4:142112/1‑148 (MQ=255)
                                                                                               ttCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCgg                                                >  5:185248/1‑149 (MQ=255)
                                                                                                   ggCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGcc                                            <  2:271979/149‑1 (MQ=255)
                                                                                                    gcgcCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGa                                                                                 >  5:78550/1‑111 (MQ=255)
                                                                                                    gcgcCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGa                                                                                 <  6:78550/111‑1 (MQ=255)
                                                                                                     cgcCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCGCTGACCATGTACCTGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCcgc                                          <  4:313142/149‑1 (MQ=255)
                                                                                                     cgcCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCcgc                                          <  4:97449/149‑1 (MQ=255)
                                                                                                       ccTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGt                                        >  1:201514/1‑149 (MQ=255)
                                                                                                         ttGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGtt                                       >  7:173832/1‑148 (MQ=255)
                                                                                                                gctgtggctgAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTcc                               <  1:190885/149‑1 (MQ=255)
                                                                                                                      gctgAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACGATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGCCCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTg                         >  7:242443/1‑149 (MQ=255)
                                                                                                                       ctgAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGa                        >  4:58909/1‑149 (MQ=255)
                                                                                                                        acaGCTTGGTGGCCTGGTCGATGACCAGGACCAGCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACtt                                                                           >  1:342752/3‑97 (MQ=21)
                                                                                                                        acaGCTTGGTGGCCTGGTCGATGACCAGGACCAGCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACtt                                                                           <  2:342752/95‑1 (MQ=21)
                                                                                                                             ttGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGa                                                                                 >  1:91327/1‑86 (MQ=255)
                                                                                                                             ttGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGa                                                                                 <  2:91327/86‑1 (MQ=255)
                                                                                                                                ctggtcctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCgg                <  3:81479/148‑1 (MQ=255)
                                                                                                                                   gtcctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGac                                                                                >  6:223164/1‑81 (MQ=255)
                                                                                                                                   gtcctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGac                                                                                <  5:223164/81‑1 (MQ=255)
                                                                                                                                   gtcctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTg            <  8:173832/149‑1 (MQ=255)
                                                                                                                                   ggcctggtcGATGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGGCCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTg            >  3:73582/3‑149 (MQ=255)
                                                                                                                                    tcctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTgg           <  1:146847/149‑1 (MQ=255)
                                                                                                                                    tcctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTgg           <  2:201514/149‑1 (MQ=255)
                                                                                                                                       tggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAg        >  1:375497/1‑149 (MQ=255)
                                                                                                                                        ggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGt              <  8:22651/142‑1 (MQ=255)
                                                                                                                                        ggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGt              >  7:22651/1‑142 (MQ=255)
                                                                                                                                        ggtcATCGACCAGGCCACCAAGCTGAACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAgc       >  6:70673/1‑149 (MQ=255)
                                                                                                                                            aTCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAGCGCTg   >  5:28887/1‑149 (MQ=255)
                                                                                                                                             tCGACCAGGCCACCAAGCTGTACTTCACCAACGCCCTGACCATGTACCAGCAGATCGTGGGCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAGCGCTgg  <  8:43956/149‑1 (MQ=255)
                                                                                                                                             tCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAGCGCTgg  >  8:279714/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAGCGCTGG  >  NC_002947/712956‑713245

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: