Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,591,709 A→T 7.5% L62H (CTT→CAT) ‡ PP_1395 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,591,7090AT 7.5% 186.3 / 3.6 63L62H (CTT→CAT) ‡PP_1395AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (34/24);  new base T (3/2);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTG  >  NC_002947/1591561‑1591854
                                                                                                                                                    |                                                                                                                                                 
gtccaTGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCaa                                                                                                                                                   <  7:357355/149‑1 (MQ=255)
gtccaTGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCaa                                                                                                                                                   <  5:178793/149‑1 (MQ=255)
   caTGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAgg                                                                                                                                                 <  5:6266/148‑1 (MQ=255)
   caTGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAgg                                                                                                                                                 <  8:34762/148‑1 (MQ=255)
   caTGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGa                                                                                                                                                <  7:236339/149‑1 (MQ=255)
       cgcAGGCCCAGGCTTGGGGTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGca                                                                                                                                            <  8:169864/149‑1 (MQ=255)
              ccAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAggg                                                                                                                                       >  2:26723/1‑147 (MQ=255)
              ccAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAggg                                                                                                                                       >  6:203034/1‑147 (MQ=255)
                     tGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAg                                                                                                                                                  >  1:190257/1‑129 (MQ=255)
                     tGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAg                                                                                                                                                  <  2:190257/129‑1 (MQ=255)
                        gTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTgcgc                                                                                                                             >  3:189990/1‑147 (MQ=255)
                            tcGGGGCTGGCATCTAGCCTCAAGGCGGTGAGTGGCAATGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgggt                                                                                                                       <  8:340289/149‑1 (MQ=255)
                            tcGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgggt                                                                                                                       >  6:281251/1‑149 (MQ=255)
                              ggggCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgggt                                                                                                                       <  3:232783/147‑1 (MQ=255)
                               gggCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgggtggg                                                                                                                    <  1:26723/149‑1 (MQ=255)
                                  cTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTa                                                                                                                  >  8:63667/1‑148 (MQ=255)
                                     gCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACa                                                                                                                <  5:203034/147‑1 (MQ=255)
                                              cTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCaa                                                                                                                                                   >  4:131528/1‑103 (MQ=255)
                                              cTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCaa                                                                                                                                                   <  3:131528/103‑1 (MQ=255)
                                                  cGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGAt                                                                                                                                               <  1:82256/103‑1 (MQ=255)
                                                  cGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGAt                                                                                                                                               >  2:82256/1‑103 (MQ=255)
                                                  cGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCg                                                                                                 <  6:402526/149‑1 (MQ=255)
                                                    gcgcTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGgt                                                                                               <  7:145248/149‑1 (MQ=255)
                                                    gcgcTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGgt                                                                                               <  7:63667/149‑1 (MQ=255)
                                                        tGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTaa                                                                                           >  7:356218/1‑149 (MQ=255)
                                                              gCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCAc                                                                                      >  1:402237/1‑148 (MQ=255)
                                                               cAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCa                                                                                    <  4:24119/149‑1 (MQ=255)
                                                                      aCCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAg                                                                             >  7:308672/1‑149 (MQ=255)
                                                                        cGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAgg                                                                                  >  3:224822/1‑142 (MQ=255)
                                                                        cGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAgg                                                                                  <  4:224822/142‑1 (MQ=255)
                                                                        cGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAgg                                                                            >  5:404030/1‑148 (MQ=255)
                                                                        cGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGt                                                                           >  4:198308/1‑149 (MQ=255)
                                                                                ccATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATg                                                                   <  5:281251/149‑1 (MQ=255)
                                                                                ccATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATg                                                                   >  1:168515/1‑149 (MQ=255)
                                                                                             tGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGAc                                                      <  8:356218/149‑1 (MQ=255)
                                                                                             tGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCTTGGCAAAGAAGCACAAGGCAGGCAGCGCCAGGGTTGGTACAGAAAGTACATTTTAGGTGTTACGCACCAGGGACAAGTCATAGAATGCCCAAGCTTGAc                                                      >  2:136490/1‑148 (MQ=255)
                                                                                                 tCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCa                                                    >  6:277304/1‑147 (MQ=255)
                                                                                                 tCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCa                                                    >  2:93873/1‑147 (MQ=255)
                                                                                                     aGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGc                                              >  8:327011/1‑149 (MQ=255)
                                                                                                       gCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTc                                            >  2:229806/1‑149 (MQ=255)
                                                                                                            cTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgg                                                                                                                         >  4:265967/1‑67 (MQ=255)
                                                                                                            cTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAgg                                                                                                                         <  3:265967/67‑1 (MQ=255)
                                                                                                                tCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATg                                                                   <  2:278900/117‑1 (MQ=255)
                                                                                                                tCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATg                                                                   >  1:278900/1‑117 (MQ=255)
                                                                                                                tCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATg                                   >  3:95082/1‑149 (MQ=255)
                                                                                                                tCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATg                                   >  2:32615/1‑149 (MQ=255)
                                                                                                                 cGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATgg                                  >  1:358339/1‑149 (MQ=255)
                                                                                                                  gTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTAc                                                                                                                 >  2:328632/1‑69 (MQ=255)
                                                                                                                  gTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTAc                                                                                                                 <  1:328632/69‑1 (MQ=255)
                                                                                                                  gTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGCCCTAAGGGCTCGCCCGCAAGGc                                 >  8:384333/1‑149 (MQ=255)
                                                                                                                  gTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGc                                 <  5:213334/149‑1 (MQ=255)
                                                                                                                     aCCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCgg                               <  4:189990/148‑1 (MQ=255)
                                                                                                                       ccAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCACCCGCATGGCGGtt                             >  1:236927/1‑148 (MQ=255)
                                                                                                                           gAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCa                        <  2:168515/149‑1 (MQ=255)
                                                                                                                               aGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTc                    >  6:204831/1‑149 (MQ=255)
                                                                                                                                 gcgcagcctgcCCTGTGCATCCTTGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCgg                  >  6:335432/8‑149 (MQ=38)
                                                                                                                                 gCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCgg                  >  4:1846/1‑149 (MQ=255)
                                                                                                                                  cgcagcctgcCCTGTGCATCCTTGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGc                 <  5:335432/143‑1 (MQ=255)
                                                                                                                                   ttttGCTGCCCTGGGCAAGGATGCACAGGGCAGGCAGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGc                 >  1:268745/1‑148 (MQ=255)
                                                                                                                                     agcctgcCCTGTGCATCCTTGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAg              >  6:172323/4‑149 (MQ=255)
                                                                                                                                       gctgcCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCt            >  7:29436/1‑149 (MQ=255)
                                                                                                                                       cctgcCCTGTGCATCCTTGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTg                                                         <  8:200816/103‑1 (MQ=38)
                                                                                                                                       cctgcCCTGTGCATCCTTGCCCAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTg                                                         >  7:200816/2‑104 (MQ=38)
                                                                                                                                         tgcCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTc           >  6:75856/1‑148 (MQ=255)
                                                                                                                                                  cAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGCTCAGCTCGGCCAGCTCCAGGCGCTg  >  8:92260/1‑148 (MQ=255)
                                                                                                                                                  cAAGGATGCACAGGGCAGGCTGCACCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCaaa                                                  >  5:357521/1‑100 (MQ=255)
                                                                                                                                                  cAAGGATGCACAGGGCAGGCTGCACCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCaaa                                                  <  6:357521/100‑1 (MQ=255)
                                                                                                                                                   aaGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTg  >  3:196091/1‑147 (MQ=255)
                                                                                                                                                   aaGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTg  <  4:196091/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGCTGAGTGGCAAGGTCACCGAAACGACCATGTAATGCAGTGGGTCGTAGGCGTACTGCTCGTCACCCAGGAACAGGCTTTTGCTGCCCTGGGCAAGGATGCACAGGGCAGGCTGCGCCAGGGTGGGTACAGAACGTACATTTTCGGTGTAACGCACCAGGTACAGGTCATCGATGGCCGATGCTGGACCAAAGGGCTCGCCCGCATGGCGGTTGATCAGCTCGGCCAGCTCCAGGCGCTG  >  NC_002947/1591561‑1591854

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: