Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 624,537 T→G 5.1% K125Q (AAG→CAG)  PP_0537 ← transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947624,5370TG 5.1% 228.1 / 2.1 80K125Q (AAG→CAG) PP_0537transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (40/35);  new base G (2/2);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

GCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGTTGGCGGGCGCTGGCCGTGGACGGCC  >  NC_002947/624396‑624684
                                                                                                                                             |                                                                                                                                                   
gcagGGCATAGGTTTCGCCCTCGACCACCTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCtt                                                                                                                                             >  1:164362/1‑150 (MQ=255)
 cagGGCATAGGTTTCGCCCTCGACCACCTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCttt                                                                                                                                            >  2:557177/1‑150 (MQ=255)
  agGGCATAGGTTTCGCCCTCGACCACCTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTACCGGCTTGCACACCCAAGCCCTTCAGGGCtttt                                                                                                                                           <  2:1130360/150‑1 (MQ=255)
     gCATAGGTTTCGCCCTCGACCACCTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGcc                                                                                                                                        <  1:190340/150‑1 (MQ=255)
      cATAGGTTTCGCCCTCGACCACCTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCg                                                                                                                                       >  2:416688/1‑150 (MQ=255)
            tttCGCCCTCGACCACCTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAgc                                                                                                                                 >  1:69852/1‑150 (MQ=255)
                gCCCTCGACCACCTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAgcg                                                                                                                                <  1:1224047/147‑1 (MQ=255)
                gCCCTCGACCACCTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAgcg                                                                                                                                >  2:1224047/1‑147 (MQ=255)
                gCCCTCGACCACCTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAg                                                                                                                             <  1:766901/150‑1 (MQ=255)
                gCCCTCGACCACCTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAg                                                                                                                             >  2:984798/1‑150 (MQ=255)
                gCCCTCGACCACCTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAg                                                                                                                             <  1:284782/150‑1 (MQ=255)
                gCCCTCGACCACCTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAg                                                                                                                             >  2:869197/1‑150 (MQ=255)
                  ccTCGACCACCTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGt                                                                                                                           <  2:676268/150‑1 (MQ=255)
                    tCGACCACCTGGGTTAGGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGtgc                                                                                                                         >  2:53994/1‑150 (MQ=255)
                    ggcgaaaccTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGcc                                                                                                                                        >  2:1149566/7‑135 (MQ=255)
                    ggcgaaaccTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGcc                                                                                                                                        <  1:1149566/129‑1 (MQ=255)
                      gACCACCTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGtgctt                                                                                                                       <  2:958227/150‑1 (MQ=255)
                         caccTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGtgcttgct                                                                                                                    >  1:26400/1‑150 (MQ=255)
                           ccTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTtgctgc                                                                                                                  <  1:122911/150‑1 (MQ=255)
                            cTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGGCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGcc                                                                                                                 >  1:793103/1‑150 (MQ=255)
                            cTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGcc                                                                                                                 <  1:650619/150‑1 (MQ=255)
                               ggTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCAc                                                                                                              <  2:1033616/150‑1 (MQ=255)
                                 tCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCg                                                                                                            >  2:1200464/1‑150 (MQ=255)
                                   aTGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGc                                                                                                          <  1:557177/150‑1 (MQ=255)
                                   aTGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGc                                                                                                          <  1:1200464/150‑1 (MQ=255)
                                     gCCCAGGTCGGTGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTCATTGCCTGCATTGCCAGGCAGATTGCGTTGCCGGCTTACACCCCGAATCCAAACAGGGCTTTTACCTGCCAGCGCAGCTGCTTGCTGCCCAccctcca                                                                                                        >  2:247966/1‑145 (MQ=255)
                                       ccAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACa                                                                                                      <  2:121405/150‑1 (MQ=255)
                                       ccAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACa                                                                                                      >  1:898425/1‑150 (MQ=255)
                                       ccAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACa                                                                                                      <  1:1078243/150‑1 (MQ=255)
                                           ggcGATGGCTAGGGTAGGTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGGAGATTGCGTTGCCGGCTTGCACACCCAAGCGCTTCTGGGCTTTTGCCGTCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCaa                                                                                                  <  1:195297/148‑1 (MQ=255)
                                           gTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCaa                                                                                                  <  1:869197/150‑1 (MQ=255)
                                           gTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCaa                                                                                                  >  2:243422/1‑150 (MQ=255)
                                           gTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCaa                                                                                                  <  2:898425/150‑1 (MQ=255)
                                                  gCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTGAAGGGCttgggt                                                                                                                                         <  2:503730/104‑5 (MQ=255)
                                                  gCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTGAAGGGCttgggt                                                                                                                                         >  1:503730/1‑100 (MQ=255)
                                                       ggTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCg                                                                                                            >  1:1138602/1‑128 (MQ=255)
                                                       ggTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGACGGCTTTCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCg                                                                                                            <  2:1138602/128‑1 (MQ=255)
                                                          aGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTc                                                                                   <  2:26400/150‑1 (MQ=255)
                                                                  gCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAgg                                                                           <  2:307821/150‑1 (MQ=255)
                                                                        tCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCACGTGCTTACTGACCGCCGGCGCCAGCAAGCCGTTAAGCACTTCGTAGTAGGGCAACt                                                                     >  2:566573/1‑150 (MQ=255)
                                                                         cAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCtt                                                                    <  1:53994/150‑1 (MQ=255)
                                                                            cGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCaa                                                                 <  1:243422/150‑1 (MQ=255)
                                                                                  cATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCg                                                           <  1:416688/150‑1 (MQ=255)
                                                                                            gCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCAcc                                                                      >  2:990081/1‑129 (MQ=255)
                                                                                            gCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCAcc                                                                      <  1:990081/129‑1 (MQ=255)
                                                                                            gCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGt                                                 >  2:911862/1‑150 (MQ=255)
                                                                                            gCCGGCATTGCCAGGCAGATTGCGTTGCCAGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGt                                                 >  1:881147/1‑150 (MQ=255)
                                                                                              cGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCa                                                                  >  2:1261621/1‑131 (MQ=255)
                                                                                              cGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCa                                                                  <  1:1261621/131‑1 (MQ=255)
                                                                                                gCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCg                                             >  1:490915/1‑150 (MQ=255)
                                                                                                  aTTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGAc                                           >  2:296582/1‑150 (MQ=255)
                                                                                                        aGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCggg                                     >  2:1219855/1‑150 (MQ=255)
                                                                                                         ggCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCgta                                                                                <  1:847143/106‑1 (MQ=255)
                                                                                                         ggCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCgta                                                                                >  2:847143/1‑106 (MQ=255)
                                                                                                         ggCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGGACCTTCAATACCCGGCCACGATGTTCCGACAGCgggg                                    >  2:302513/1‑150 (MQ=255)
                                                                                                          gCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCaa                                                                                                  <  1:1082609/87‑1 (MQ=255)
                                                                                                          gCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCaa                                                                                                  >  2:1082609/1‑87 (MQ=255)
                                                                                                            aGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTaa                                 >  2:805966/1‑150 (MQ=255)
                                                                                                             gATTGCGTTGCCGGCTTGCACACCCAAACCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAg                                >  1:236545/1‑150 (MQ=255)
                                                                                                              aTTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGc                               >  1:379366/1‑150 (MQ=255)
                                                                                                                tGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCgg                             <  2:69852/150‑1 (MQ=255)
                                                                                                                   gTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGtt                          >  2:982595/1‑150 (MQ=255)
                                                                                                                    ttGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGTTg                         <  2:379366/150‑1 (MQ=255)
                                                                                                                       ccGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGTTggcg                      >  1:789971/1‑150 (MQ=255)
                                                                                                                        cGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCGTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGTTggcgg                     >  1:353327/1‑150 (MQ=255)
                                                                                                                         ggCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCg                                                           >  1:933183/1‑111 (MQ=255)
                                                                                                                         ggCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCg                                                           <  2:933183/111‑1 (MQ=255)
                                                                                                                          gCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTCGTCAGCCACCTTCAATACCCGGGGTCGAGGTTCCGACAGCGGGGTAAGCGGGTTggcgggc                   >  2:943485/1‑150 (MQ=255)
                                                                                                                          gCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGTTggcg                      >  1:22133/1‑147 (MQ=255)
                                                                                                                          gCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGTTggcg                      <  2:22133/147‑1 (MQ=255)
                                                                                                                                  acCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGTTGGCGGGCGCTGGCCg           >  2:904850/1‑150 (MQ=255)
                                                                                                                                  acCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGTTGGCGGGCGCTGGCCg           >  1:61535/1‑150 (MQ=255)
                                                                                                                                   cccAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGTTGGCGGGCGCTGGCCGt          <  2:639836/150‑1 (MQ=255)
                                                                                                                                   cccAAGCCCTTCAGGGCTTGTGCCGGCAATCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATATTCCTCCAGCGGGGTAAGCGGGTTGGCGGGCGCTGGCCGt          <  1:247966/150‑1 (MQ=255)
                                                                                                                                    ccAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGTTGGCGGGCGCTGGCCGTg         >  1:181219/1‑150 (MQ=255)
                                                                                                                                    ccAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGTTGGCGGGCGCTGGCCGTg         >  2:39901/1‑150 (MQ=255)
                                                                                                                                    aaaaGCCCTGAAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGc                                        >  2:865659/3‑119 (MQ=255)
                                                                                                                                    aaaaGCCCTGAAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGc                                        <  1:865659/117‑1 (MQ=255)
                                                                                                                                        gCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGTTGGCGGGCGCTGGCCGTGGACg     >  2:979102/1‑150 (MQ=255)
                                                                                                                                           cTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGTTGGCGGGCGCTGGCCGTGGACGGcc  <  2:353327/150‑1 (MQ=255)
                                                                                                                                            ttCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGTTGGCGGGCGCTGGCCg           <  2:692728/140‑1 (MQ=255)
                                                                                                                                            ttCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGTTGGCGGGCGCTGGCCg           >  1:692728/1‑140 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATGCCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCAGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCCCACCGGCGACAGCAAGCCGTTGCGCACTTCGTAGTAGGGCACCTTCAATACCCGGCCACGATGTTCCGACAGCGGGGTAAGCGGGTTGGCGGGCGCTGGCCGTGGACGGCC  >  NC_002947/624396‑624684

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: