Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,507,063 C→A 5.3% P92P (CCC→CCA cycA → D‑alanine/beta‑alanine/D‑serine/glycine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,507,0630CA 5.3% 364.2 / 5.1 113P92P (CCC→CCAcycAD‑alanine/beta‑alanine/D‑serine/glycine permease
Reads supporting (aligned to +/- strand):  ref base C (57/50);  new base A (2/4);  total (59/54)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 118 total aligned reads displayed.

ACCATCAGCCTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAGGGCTGCTGCT  >  NC_002947/5506914‑5507212
                                                                                                                                                        |                                                                                                                                                     
acgagcAGCCGGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGcccc                                                                                                                                                       <  2:1450912/145‑1 (MQ=255)
   aTCAGCCTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTgg                                                                                                                                                    <  2:742090/150‑1 (MQ=255)
        ccTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGATTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCgg                                                                                                                                               >  1:1362771/1‑150 (MQ=255)
         cTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCTCCTGGGCCGGc                                                                                                                                              >  2:707731/1‑150 (MQ=255)
         cTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGc                                                                                                                                              >  1:1929968/1‑150 (MQ=255)
                  ccATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTAct                                                                                                                                     <  1:707731/150‑1 (MQ=255)
                   cATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTActg                                                                                                                                    <  2:704391/150‑1 (MQ=255)
                    aTCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTActgg                                                                                                                                   <  2:1929968/150‑1 (MQ=255)
                    aTCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTActgg                                                                                                                                   >  1:275008/1‑150 (MQ=255)
                     tCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTActggc                                                                                                                                  >  1:1422383/1‑150 (MQ=255)
                      cGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTActggct                                                                                                                                 <  1:968977/150‑1 (MQ=255)
                         tcatgTTCGTCTTCATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTCGGCCGGCTACTTTACTGGCTGGa                                                                                                                              <  2:1194191/150‑1 (MQ=255)
                          catcTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTgg                                                                                                                                                    <  1:150972/127‑1 (MQ=255)
                          catcTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGAc                                                                                                                             <  1:815650/150‑1 (MQ=255)
                            tcTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTAc                                                                                                                                           <  2:1565529/134‑1 (MQ=255)
                            tcTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTAc                                                                                                                                           >  1:1565529/1‑134 (MQ=255)
                            tcTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCt                                                                                                                           >  2:319192/1‑150 (MQ=255)
                             cTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTa                                                                                                                          >  1:1794579/1‑150 (MQ=255)
                              ttCGTCTACATGATTATAGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTAc                                                                                                                         <  2:2220121/150‑1 (MQ=255)
                                 gTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTgg                                                                                                                      >  2:715775/1‑150 (MQ=255)
                                      cATGATTATCGGCTTCATGCTGTTCTTGGTCATGCGAGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTATCAAACCGCCCTGCT‑‑‑CCCCTGGGCCGGCCACTTTACTGGCTGGACCTACTGGTTctg                                                                                                                 >  2:664713/1‑150 (MQ=255)
                                           ttATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTggg                                                                                                            >  2:488026/1‑150 (MQ=255)
                                           ttATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTggg                                                                                                            <  1:484769/150‑1 (MQ=255)
                                            tATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGggt                                                                                                           >  2:1862236/1‑150 (MQ=255)
                                             agcGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTAATTTATGGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGggtt                                                                                                          <  1:157899/148‑2 (MQ=255)
                                             aTCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGggtg                                                                                                          <  1:2084038/150‑1 (MQ=255)
                                             aTCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGggtg                                                                                                          <  1:994744/150‑1 (MQ=255)
                                             aTCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCAGGGGCCgagca                                                                                                                                            <  1:1717913/116‑5 (MQ=255)
                                              tCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGggtgg                                                                                                         >  2:275210/1‑150 (MQ=255)
                                                  cTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCAc                                                                                                     >  1:821363/1‑150 (MQ=255)
                                                  cTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCAc                                                                                                     >  2:168544/1‑150 (MQ=255)
                                                    tCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTTGTTCTGCTGGGTGGTCACCg                                                                                                   <  2:1549632/150‑1 (MQ=255)
                                                       tGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCa                                                                                                <  1:715775/150‑1 (MQ=255)
                                                       tGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCa                                                                                                >  1:1085987/1‑150 (MQ=255)
                                                        gCTNTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCAt                                                                                               >  2:1517991/1‑150 (MQ=255)
                                                         cTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCagg                                                                                              <  2:1413182/150‑3 (MQ=255)
                                                          tGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACACGGCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGACCCCTGAACCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCAtcg                                                                                             >  2:1125806/1‑148 (MQ=255)
                                                          tGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTg                                                                                             >  2:838417/1‑150 (MQ=255)
                                                           gTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACATGAGCTAATGCTGGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGc                                                                                            >  1:1995507/1‑150 (MQ=255)
                                                            ttcttcGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGcc                                                                                           >  1:56247/1‑150 (MQ=255)
                                                            ttcttcGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGcc                                                                                           >  1:486918/1‑150 (MQ=255)
                                                               ttcGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGAc                                                                                        <  2:1794579/150‑1 (MQ=255)
                                                               ttcGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGAc                                                                                        <  1:838417/150‑1 (MQ=255)
                                                                  gTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTa                                                                                                                                            <  1:1654001/95‑1 (MQ=255)
                                                                  gTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTg                                                                                     <  1:463179/150‑1 (MQ=255)
                                                                  gTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTGTACTTTCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTg                                                                                     >  1:2380862/1‑150 (MQ=255)
                                                                   tCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTgg                                                                                    <  2:529883/150‑1 (MQ=255)
                                                                   tCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTgg                                                                                    <  1:1088271/150‑1 (MQ=255)
                                                                     aTGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCTCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCAt                                                                                               <  1:758926/137‑1 (MQ=255)
                                                                       gcgcgcCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTcgc                                                                                <  2:1524157/150‑1 (MQ=255)
                                                                       gcgcgcCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTcgc                                                                                <  2:1524183/150‑1 (MQ=255)
                                                                           gcCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATc                                                                            <  1:451074/150‑1 (MQ=255)
                                                                             cATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCAGGGGCCgagca                                                                                                                                            >  2:345236/1‑80 (MQ=255)
                                                                             cATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCAGGGGCCgagca                                                                                                                                            <  1:345236/84‑5 (MQ=255)
                                                                              aTGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATcgcc                                                                         >  2:454309/1‑150 (MQ=255)
                                                                              aTGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATcgcc                                                                         >  2:1694459/1‑150 (MQ=255)
                                                                                   cGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTTGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTa                                                                    <  2:2157199/150‑1 (MQ=255)
                                                                                       ctgctgctgTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTgg                                                                                                                                                    >  1:1886230/1‑66 (MQ=255)
                                                                                       ctgctgctgTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTgg                                                                                                                                                    <  2:1886230/66‑1 (MQ=255)
                                                                                       ctgctgctgTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGc                                                                                                                                              >  2:2293697/1‑72 (MQ=255)
                                                                                       ctgctgctgTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGc                                                                                                                                              <  1:2293697/72‑1 (MQ=255)
                                                                                         gctgctgTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACAgccg                                                              <  2:221609/150‑2 (MQ=255)
                                                                                          ctgctgTCGATCCTGAACGACAAGGCGTTGATCGACGTT‑‑TCCGCCGACCTGCGCGGCCGCTGGGCCGGGTAGTTTACTGGCTGGACCTACTGGTTCTTCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAg                                                             >  2:821248/1‑150 (MQ=255)
                                                                                           tgctgTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACATGCGCGGCCCCTGGGGCGGCTACTTTACTCGCTGGACATACGGGTGCTGCTGGGTGGTCACCGGCAGTGCCGACGTGGTCGCCCTCGCCGCCTACACCCAGt                                                            >  2:253572/1‑150 (MQ=255)
                                                                                               gTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTg                                                        <  1:1517991/150‑1 (MQ=255)
                                                                                                tCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCGCGGGATGTGGCGTTCTGGCCCAACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTgg                                                       <  1:168544/150‑1 (MQ=255)
                                                                                                  gAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTgg                                                                                    >  1:491258/1‑119 (MQ=255)
                                                                                                  gAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTgg                                                                                    <  2:491258/119‑1 (MQ=255)
                                                                                                  gAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGtt                                                     >  1:1157639/1‑150 (MQ=255)
                                                                                                    aCCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTcc                                                   <  2:275008/150‑1 (MQ=255)
                                                                                                     ccTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGggtg                                                                                                          <  2:1371738/94‑1 (MQ=255)
                                                                                                     ccTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTccc                                                  >  1:507472/1‑150 (MQ=255)
                                                                                                      cTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCAGGGGCCga                                                                                                                                               <  1:764817/56‑2 (MQ=255)
                                                                                                      cTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCAGGGGCCga                                                                                                                                               >  2:764817/1‑55 (MQ=255)
                                                                                                       tGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCg                                                >  1:1404405/1‑150 (MQ=255)
                                                                                                       tGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCg                                                >  1:115302/1‑150 (MQ=255)
                                                                                                           cTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGAc                                                                                        <  2:389979/106‑1 (MQ=255)
                                                                                                           cTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCt                                            >  1:592424/1‑150 (MQ=255)
                                                                                                            tACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTg                                           >  2:2243439/1‑150 (MQ=255)
                                                                                                                   cGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGt                                    <  2:1362771/150‑1 (MQ=255)
                                                                                                                         tCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGAc                                                                                                                             <  2:2059649/55‑1 (MQ=255)
                                                                                                                         tCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGAc                                                                                                                             >  1:2059649/1‑55 (MQ=255)
                                                                                                                         tCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATAc                              >  2:811785/1‑150 (MQ=255)
                                                                                                                         tCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATAc                              >  1:1685302/1‑150 (MQ=255)
                                                                                                                         tCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGAGGCCGCGATCACCGCCTACACCCAGTTCTGGTTGCCCCACCTGCCGCAGTGGATAc                              >  2:744649/1‑150 (MQ=255)
                                                                                                                              ttt‑‑tCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGgcg                         >  2:2056637/1‑150 (MQ=255)
                                                                                                                                        ccgACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCt                  >  1:809022/1‑150 (MQ=255)
                                                                                                                                        ccgACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCt                  >  1:695709/1‑150 (MQ=255)
                                                                                                                                          gACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTCCTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGc                <  1:1299521/150‑1 (MQ=255)
                                                                                                                                          gACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGc                <  2:2380862/150‑1 (MQ=255)
                                                                                                                                             cTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTa             >  2:185268/1‑150 (MQ=255)
                                                                                                                                             cTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTa             >  2:1133624/1‑150 (MQ=255)
                                                                                                                                             cTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTa             <  2:809022/150‑1 (MQ=255)
                                                                                                                                             cTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTa             >  1:95948/1‑150 (MQ=255)
                                                                                                                                             cTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTa             >  1:96106/1‑150 (MQ=255)
                                                                                                                                               gCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAgg           <  1:275210/150‑1 (MQ=255)
                                                                                                                                                 tCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAGGgc         >  1:1870011/1‑150 (MQ=255)
                                                                                                                                                     ccccTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAGGgctgct     <  2:1571018/150‑1 (MQ=255)
                                                                                                                                                     ccccTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAGGgctgct     >  2:1988744/1‑150 (MQ=255)
                                                                                                                                                        cTGGGCCGGCTACTTTACTGGCTGGACCGACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAGGgctgctgct  <  2:455377/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
ACCATCAGCCTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTT‑‑TCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAGGGCTGCTGCT  >  NC_002947/5506914‑5507212

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: