Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,654,891 T→C 9.0% intergenic (+43/‑215) aceA → / → PP_4118 isocitrate lyase/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,654,8910TC 9.0% 215.6 / 9.7 78intergenic (+43/‑215)aceA/PP_4118isocitrate lyase/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (34/37);  new base C (5/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

CTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGCAG  >  NC_002947/4654761‑4655036
                                                                                                                                     |                                                                                                                                                 
cTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTa                                                                                                                                  <  1:40867/149‑1 (MQ=255)
cTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGttttt                                                                                                                                   <  2:6582/147‑1 (MQ=255)
 ttCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAGGAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTa                                                                                                                                  >  1:256584/1‑149 (MQ=255)
 ttCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTAcc                                                                                                                                >  7:285129/1‑149 (MQ=255)
        acATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCCGGCCTTTGTTTTTACCCgaga                                                                                                                           >  2:305763/1‑147 (MQ=255)
          aTGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgg                                                                                                                       >  1:305711/1‑149 (MQ=255)
          aTGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgg                                                                                                                       >  4:199985/1‑149 (MQ=255)
            gACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggc                                                                                                                     >  6:56727/1‑149 (MQ=255)
              ccaccGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGACCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggccg                                                                                                                   >  7:53933/1‑149 (MQ=255)
               caccGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggccgg                                                                                                                  >  8:264214/1‑149 (MQ=255)
               caccGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggccgg                                                                                                                  >  4:275803/1‑149 (MQ=255)
                      aTCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCctgct                                                                                                           <  5:211986/149‑1 (MQ=255)
                       tCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCctgctg                                                                                                          <  5:123637/149‑1 (MQ=255)
                        ccAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTgg                                                                                                         >  7:165200/1‑149 (MQ=255)
                                   cTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTAc                                                                                              <  2:149196/149‑1 (MQ=255)
                                    ttCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTAc                                                                                              >  1:257635/1‑148 (MQ=255)
                                       gtcgGTGACTGCACTGACCTGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAg                                                                                           <  1:305763/148‑1 (MQ=255)
                                       gtcgGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGc                                                                                          <  2:305711/149‑1 (MQ=255)
                                        tcgGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGcc                                                                                         >  1:248036/1‑149 (MQ=255)
                                         cgGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCg                                                                                        >  1:270105/1‑149 (MQ=255)
                                              aCTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCa                                                                                    <  4:184937/148‑1 (MQ=255)
                                              aCTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCAc                                                                                   >  1:5162/1‑149 (MQ=255)
                                               cTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCAcc                                                                                  <  2:248036/149‑1 (MQ=255)
                                                  cACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCAcc                                                                                  <  7:220247/146‑1 (MQ=255)
                                                  cACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCAcc                                                                                  >  8:220247/1‑146 (MQ=255)
                                                  cACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGtt                                                                               <  3:275803/149‑1 (MQ=255)
                                                   aCTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGttt                                                                              >  4:299004/1‑149 (MQ=255)
                                                   aCTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGttt                                                                              >  6:222537/1‑149 (MQ=255)
                                                   aCTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGtt                                                                               >  4:132327/1‑148 (MQ=255)
                                                              cGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTAc                                                                    <  3:132327/148‑1 (MQ=255)
                                                               gACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTAcc                                                                   <  2:256584/148‑1 (MQ=255)
                                                                 ccGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTa                                                                     <  5:26301/144‑1 (MQ=255)
                                                                 ccGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTa                                                                     >  6:26301/1‑144 (MQ=255)
                                                                      gaagaGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGTGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACAGTATTcgc                                                           <  2:257635/149‑1 (MQ=255)
                                                                        agagCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTc                                                                                     >  8:76779/1‑121 (MQ=255)
                                                                        agagCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTc                                                                                     <  7:76779/121‑1 (MQ=255)
                                                                        agagCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCgg                                                         <  5:56727/149‑1 (MQ=255)
                                                                        agagCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCgg                                                         <  3:199985/149‑1 (MQ=255)
                                                                           gCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGt                                                       >  1:236101/1‑148 (MQ=255)
                                                                              gTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTcgc                                                           >  3:278871/1‑141 (MQ=255)
                                                                              gTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTcgc                                                           <  4:278871/141‑1 (MQ=255)
                                                                                  cACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGc                                               <  3:299004/149‑1 (MQ=255)
                                                                                     tGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAg                                                                                           >  7:182291/1‑102 (MQ=255)
                                                                                     tGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAg                                                                                           <  8:182291/102‑1 (MQ=255)
                                                                                           tGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACt                                      <  6:312120/149‑1 (MQ=255)
                                                                                               ttGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATc                                  <  8:165200/149‑1 (MQ=255)
                                                                                               ttGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATc                                  <  5:52947/149‑1 (MQ=255)
                                                                                                  ccGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATc                                  >  4:221263/1‑146 (MQ=255)
                                                                                                  ccGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATc                                  <  3:221263/146‑1 (MQ=255)
                                                                                                   cGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGgga                              <  7:264214/149‑1 (MQ=255)
                                                                                                   cGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGgga                              <  2:5162/149‑1 (MQ=255)
                                                                                                   cGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGgga                              <  8:285129/149‑1 (MQ=255)
                                                                                                    gTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGTTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggag                             >  7:62064/1‑149 (MQ=255)
                                                                                                    gTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCCCCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggag                             >  8:181362/1‑149 (MQ=255)
                                                                                                    gTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggag                             >  1:266138/1‑149 (MQ=255)
                                                                                                       gAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggagga                           <  2:236101/148‑1 (MQ=255)
                                                                                                       aaaCA‑‑AAGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggag                             >  1:206656/2‑146 (MQ=21)
                                                                                                       aaaCA‑‑AAGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggag                             >  4:272079/2‑146 (MQ=21)
                                                                                                        aaCA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGgaga                         <  5:222537/149‑1 (MQ=255)
                                                                                                        aaCA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGgaga                         <  2:266138/149‑1 (MQ=255)
                                                                                                        aaCA‑‑AAGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGc                                                                                          >  3:293582/1‑84 (MQ=21)
                                                                                                        aaCA‑‑AAGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGc                                                                                          <  4:293582/84‑1 (MQ=21)
                                                                                                        aaCA‑‑AAGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGgaga                         >  1:291443/1‑149 (MQ=21)
                                                                                                         aCA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGac                        <  8:53933/149‑1 (MQ=255)
                                                                                                          cA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTa                                          <  4:207454/130‑1 (MQ=255)
                                                                                                          cA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTa                                          >  3:207454/1‑130 (MQ=255)
                                                                                                          cA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGaca                       <  2:51051/149‑1 (MQ=255)
                                                                                                          cA‑‑AAGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGaca                       >  4:8366/1‑149 (MQ=21)
                                                                                                               aGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCgg                                <  2:181367/137‑1 (MQ=255)
                                                                                                               aGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCgg                                >  1:181367/1‑137 (MQ=255)
                                                                                                                 gCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACAt                  >  7:157210/1‑149 (MQ=255)
                                                                                                                   cccGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCCGCCGTTCACCGTTTCCGCCGCTACAGTATTCGCGGGTAAACCCGCACCTACAATGATCGGGAGGAGACACACATc                 >  4:228580/1‑149 (MQ=255)
                                                                                                                       gggCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGc              >  6:23644/1‑148 (MQ=255)
                                                                                                                           ttGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTg         >  7:109440/1‑149 (MQ=255)
                                                                                                                              tgAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGa            >  2:293953/1‑143 (MQ=255)
                                                                                                                              tgAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGa            <  1:293953/143‑1 (MQ=255)
                                                                                                                              tgAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTggg       <  8:109440/148‑1 (MQ=255)
                                                                                                                                  aGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGCAg  <  3:8366/149‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                 
CTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACA‑‑AAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGCAG  >  NC_002947/4654761‑4655036

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