Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,964,909 C→G 6.8% S150T (AGT→ACT)  fleQ ← transcriptional regulator FleQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,964,9090CG 6.8% 234.2 / 6.8 74S150T (AGT→ACT) fleQtranscriptional regulator FleQ
Reads supporting (aligned to +/- strand):  ref base C (27/41);  new base G (2/3);  total (29/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

ACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAG  >  NC_002947/4964762‑4965054
                                                                                                                                                   |                                                                                                                                                  
aCGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACt                                                                                                                                                   >  8:237086/1‑149 (MQ=255)
 cGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTg                                                                                                                                                  <  7:47246/149‑1 (MQ=255)
   aaCGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGt                                                                                                                                                <  2:30025/149‑1 (MQ=255)
         gcTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCCCGCTGGCGGCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAAc                                                                                                                                          <  2:220590/149‑1 (MQ=255)
             cgcgGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACCAGTCGTGCCattca                                                                                                                                        <  2:434926/147‑5 (MQ=255)
                ggcgCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACCAGTCGTGCCattca                                                                                                                                        <  8:394378/144‑5 (MQ=255)
                ggcgCGTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACCAGTCGTGCCattca                                                                                                                                        >  7:394378/1‑140 (MQ=255)
                   gcTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTTCgg                                                                                                                               <  8:263944/149‑1 (MQ=255)
                       ggAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGaa                                                                                                                             <  7:79233/147‑1 (MQ=255)
                           tGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTTCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGAtt                                                                                                                       <  5:269166/149‑1 (MQ=255)
                             gTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTgg                                                                                                                     <  4:268728/149‑1 (MQ=255)
                             gTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTgg                                                                                                                     <  8:384739/149‑1 (MQ=255)
                             gTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTgg                                                                                                                     >  5:415769/1‑149 (MQ=255)
                             gTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTACTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTgg                                                                                                                     <  1:192608/149‑1 (MQ=255)
                              tAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTggg                                                                                                                    >  7:314420/1‑149 (MQ=255)
                                gTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGCCTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGt                                                                                                                   >  6:261646/1‑148 (MQ=255)
                                gTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGt                                                                                                                   <  5:19430/148‑1 (MQ=255)
                                              cTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGAccc                                                                                                    <  2:36925/149‑1 (MQ=255)
                                               tACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGAcccc                                                                                                   >  8:142989/1‑149 (MQ=255)
                                                  cACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGc                                                                                                                                 <  1:75667/116‑1 (MQ=255)
                                                  cACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGc                                                                                                                                 >  2:75667/1‑116 (MQ=255)
                                                  cACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCt                                                                                                >  8:327127/1‑149 (MQ=255)
                                                   aattCCTTGCCGGTGGCGGACTCCCCCCTGATCACCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTTCCGCCAGTGCCGAATGGCACGAATGGTGCCAACCAGGCTGCGGCACAGATTGGGTTCGCGCTGGCGACCCCGCTc                                                                                               <  2:404148/147‑1 (MQ=255)
                                                       ccTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTcgc                                                                                             <  6:160365/147‑1 (MQ=255)
                                                         ttGCCGGTGCCGGACTCCCCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCgcgcgc                                                                                         >  8:261544/1‑149 (MQ=255)
                                                         ttGACGGTGACGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCgcgcgc                                                                                         >  4:131112/1‑149 (MQ=255)
                                                           gccggtgccggACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCt                                                                                       >  4:28881/1‑149 (MQ=255)
                                                               gtgccggACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGCTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTc                                                                                   <  4:336316/149‑1 (MQ=255)
                                                                tgccggACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCa                                                                                  <  7:124333/149‑1 (MQ=255)
                                                                        tCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTc                                                                          >  4:140200/1‑149 (MQ=255)
                                                                          acacAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTcgc                                                                        >  3:84530/4‑149 (MQ=255)
                                                                          aCCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTcgc                                                                        >  2:45755/1‑149 (MQ=255)
                                                                           cccAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTcgcg                                                                       <  3:140200/149‑1 (MQ=255)
                                                                           cccAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGGGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTcgcg                                                                       <  8:84599/149‑1 (MQ=255)
                                                                            ccAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGTGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCgg                                                                      <  2:407982/149‑1 (MQ=255)
                                                                            ccAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTcgcg                                                                       >  1:203191/1‑148 (MQ=255)
                                                                            ccAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTcgcg                                                                       <  2:203191/148‑1 (MQ=255)
                                                                            ccAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCgg                                                                      <  8:314420/149‑1 (MQ=255)
                                                                             cAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCgg                                                                      >  2:40355/1‑148 (MQ=255)
                                                                                    aGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGt                                                                     >  8:387674/1‑142 (MQ=255)
                                                                                    aGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGt                                                                     <  7:387674/142‑1 (MQ=255)
                                                                                     gCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTgg                                                             >  4:408531/1‑149 (MQ=255)
                                                                                        cGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCa                                                          >  5:62338/1‑149 (MQ=255)
                                                                                          cTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACg                                                        <  1:182362/149‑1 (MQ=255)
                                                                                             gCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTg                                                     <  3:271627/149‑1 (MQ=255)
                                                                                             gCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTg                                                     >  8:394991/1‑149 (MQ=255)
                                                                                             gCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACATGGGCACGGt                                                      >  8:400627/1‑148 (MQ=255)
                                                                                               gTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTg                                                     >  8:43187/1‑147 (MQ=255)
                                                                                                tCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAg                                                  <  3:131112/149‑1 (MQ=255)
                                                                                                  ggTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGGGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCg                                                <  4:154755/149‑1 (MQ=255)
                                                                                                    tATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGTAGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGc                                                 >  1:349004/1‑147 (MQ=255)
                                                                                                        cGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTc                                            >  7:19995/1‑147 (MQ=255)
                                                                                                                 ctgcATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTgg                                  <  6:62338/148‑1 (MQ=255)
                                                                                                                 ctgcATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTgg                                  <  7:327127/148‑1 (MQ=255)
                                                                                                                    catcatCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGt                              <  6:19798/149‑1 (MQ=255)
                                                                                                                     atcatcTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGGACGGTGCAGCGAGTCCAGCAATTGGCTGTa                             <  1:40355/149‑1 (MQ=255)
                                                                                                                     atcatcTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTg                               >  7:342220/1‑147 (MQ=255)
                                                                                                                     atcatcTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTa                             <  3:408531/149‑1 (MQ=255)
                                                                                                                             ccGCAGGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTggcg                     <  2:140185/149‑1 (MQ=255)
                                                                                                                             ccGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTggcg                     <  1:421353/149‑1 (MQ=255)
                                                                                                                             ccGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTggc                      <  7:394991/148‑1 (MQ=255)
                                                                                                                             ccGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTggc                      <  7:400627/148‑1 (MQ=255)
                                                                                                                               gCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTggcggc                   <  8:283613/149‑1 (MQ=255)
                                                                                                                               gCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTggcgg                    >  6:248326/1‑148 (MQ=255)
                                                                                                                                 aCGTGCTGAATGGCACGACTGGTGCCAACAAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTGGCGGCAt                 >  2:286952/1‑149 (MQ=255)
                                                                                                                                        ggttGGCACCAGTCGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAAtt                                    <  5:435866/120‑1 (MQ=255)
                                                                                                                                        ggttGGCACCAGTCGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAAtt                                    >  6:435866/4‑123 (MQ=255)
                                                                                                                                             gCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTGGCGGCATTTCCAGGTTgg      <  4:84530/148‑1 (MQ=255)
                                                                                                                                             gCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTGGCGGCATTTCCAGGTTGGa     <  1:69426/149‑1 (MQ=255)
                                                                                                                                              cACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTGGCGGCATTTCCAGGTTGGaa    <  3:28881/149‑1 (MQ=255)
                                                                                                                                               aCGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTGGCGGCATTTCCAGGTTGGaaa   <  5:261646/149‑1 (MQ=255)
                                                                                                                                                cGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAg  <  6:415769/149‑1 (MQ=255)
                                                                                                                                                 gACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTGGCGGCATTTCCAGGTTGGaa    <  4:239257/146‑1 (MQ=255)
                                                                                                                                                 gACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTGGCGGCATTTCCAGGTTGGaa    >  3:239257/1‑146 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGACCTGGGCACGGTGCAGCGAGTCCAGCAATTGGCTGTAGCTTGGCGGCATTTCCAGGTTGGAAAG  >  NC_002947/4964762‑4965054

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: