Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 156,385 G→A 6.1% intergenic (‑10/+33) ygdQ ← / ← citN transport protein/citrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947156,3850GA 6.1% 240.8 / 2.6 81intergenic (‑10/+33)ygdQ/citNtransport protein/citrate transporter
Reads supporting (aligned to +/- strand):  ref base G (43/33);  new base A (3/2);  total (46/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTCGATA  >  NC_002947/156246‑156532
                                                                                                                                           |                                                                                                                                                     
ccAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGgcgc                                                                                                                                              >  1:339416/1‑149 (MQ=255)
 cAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCTGACCGGGGcgcg                                                                                                                                             >  8:102282/1‑149 (MQ=255)
   gAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCTGACCGGGGcgcgtt                                                                                                                                           <  7:102282/149‑2 (MQ=255)
              atgatgatgTTGTCGATACCGAGGACGATCTCAAGCGAGGTAACGGTGAAAAAGGCAACCAAGATTTCCGGGCTGGTCAGACATTCCATGAATTTTCCTTCGAGCGGTAATGGCAACTCGCCCCGGTAACGGAGCGTCGCTTACGtaa                                                                                                                                 >  7:247872/1‑146 (MQ=255)
                gatgatgTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATAc                                                                                                                              >  5:13667/1‑149 (MQ=255)
                  tgatgTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACgg                                                                                                                            <  6:13667/149‑1 (MQ=255)
                   gatgTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTAc                                                                                                                                     >  1:86614/1‑139 (MQ=255)
                   gatgTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTAc                                                                                                                                     <  2:86614/139‑1 (MQ=255)
                                 cGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGacta                                                                                                           <  2:339416/149‑1 (MQ=255)
                                    ggACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGactact                                                                                                         >  1:221196/1‑148 (MQ=255)
                                       cGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGaa                                                                                                      >  1:172607/1‑148 (MQ=255)
                                       cGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAAc                                                                                                     <  2:143667/149‑1 (MQ=255)
                                        gATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGaa                                                                                                      >  4:243730/1‑147 (MQ=255)
                                        gATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACa                                                                                                    >  8:184734/1‑149 (MQ=255)
                                          tctcAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGTATGTAT‑‑Taaa                                                                                                                <  8:214237/135‑1 (MQ=255)
                                          tctcAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑Taaa                                                                                                                >  7:214237/1‑135 (MQ=255)
                                          tctcAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGc                                                                                                  >  7:214628/1‑149 (MQ=255)
                                              aaGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCggg                                                                                               >  8:34354/1‑148 (MQ=255)
                                                    ggTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATcc                                                                                        <  7:184734/149‑1 (MQ=255)
                                                     gTATGGGTGAAAACGCAAACCCTGATTTCCGGGCTGGTAAGCCATTCCTTGTATTTTCCTTCTAGCGGTATTGGAAACGCGCCCCGGTCAGGGCGCGTCGCGTAGGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATccc                                                                                       <  1:237947/149‑1 (MQ=255)
                                                      tAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATccc                                                                                       >  1:157677/1‑148 (MQ=255)
                                                        aGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCa                                                                                     <  2:172607/148‑1 (MQ=255)
                                                             gAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAG                                                                                 >  2:100393‑M1/1‑147 (MQ=255)
                                                               aaaaGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACgg                                                                                                                            >  2:52466/1‑104 (MQ=255)
                                                               aaaaGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCACATCAGCAGT                                                                             >  8:139748‑M1/1‑149 (MQ=255)
                                                                      aaCCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGcc                                                                        >  6:134751/1‑147 (MQ=255)
                                                                      aaCCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGcc                                                                        >  5:142656/1‑147 (MQ=255)
                                                                      aaCCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGa                                                                      <  3:243730/149‑1 (MQ=255)
                                                                         ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCa                                                                   >  8:42882/1‑149 (MQ=255)
                                                                         ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCa                                                                   >  7:301773/1‑149 (MQ=255)
                                                                         ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCa                                                                   >  7:8155/1‑149 (MQ=255)
                                                                         ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCa                                                                   >  5:86279/1‑149 (MQ=255)
                                                                         ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCa                                                                   >  7:9510/1‑149 (MQ=255)
                                                                         ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCa                                                                   >  8:110538/1‑149 (MQ=255)
                                                                         ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCACCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCa                                                                   >  1:16788/1‑149 (MQ=255)
                                                                              tttCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATg                                                              <  6:86279/149‑1 (MQ=255)
                                                                                         gTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTCTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGc                                                       <  1:100393/147‑1 (MQ=255)
                                                                                          tCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATAc                                                   >  5:19884/1‑148 (MQ=255)
                                                                                            aGCAATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAg                                                <  8:168173‑M2/149‑2 (MQ=255)
                                                                                              ccATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAgca                                              >  7:34354‑M2/1‑146 (MQ=255)
                                                                                               cATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGAACAGGCATACCAgcac                                             <  7:189467‑M2/149‑5 (MQ=255)
                                                                                                aTTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAgcact                                            <  1:184183‑M2/149‑6 (MQ=255)
                                                                                                  tCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCAt                                                     <  6:324148/138‑1 (MQ=255)
                                                                                                  tCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCAt                                                     >  5:324148/1‑138 (MQ=255)
                                                                                                    cATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCtggtatgcctgtgcataatgctcggcgcactgct                                         >  3:234861‑M1/1‑114 (MQ=255)
                                                                                                     aTGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCtggtatgcctgtgcataatgctcggcgcactgctga                                       <  2:157677‑M1/149‑37 (MQ=255)
                                                                                                     aTGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAgg                                                        >  4:39201/1‑132 (MQ=255)
                                                                                                     aTGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCa                                       <  7:168173/149‑1 (MQ=255)
                                                                                                     aTGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCa                                       >  8:245716/1‑149 (MQ=255)
                                                                                                      tGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCtggtatgcctgtgcataatgctcggcgcactgctgat                                      <  4:234861‑M1/149‑38 (MQ=255)
                                                                                                      tGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGa                                                                      <  4:36358/117‑1 (MQ=255)
                                                                                                         aTTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACtt                                    >  7:177377/1‑148 (MQ=255)
                                                                                                         aTTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACtt                                    >  5:19451/1‑148 (MQ=255)
                                                                                                         aTTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACtt                                    >  3:328570/1‑148 (MQ=255)
                                                                                                         aTTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACtt                                    >  3:236961/1‑148 (MQ=255)
                                                                                                                  cGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGaa                           <  8:214628/148‑1 (MQ=255)
                                                                                                                  cGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGa                            <  6:19451/147‑1 (MQ=255)
                                                                                                                   gAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCCGCACTGCCCACTTGAGCGTGAAgc                         >  3:39082/1‑149 (MQ=255)
                                                                                                                     gCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAgcgc                       <  7:42882/149‑1 (MQ=255)
                                                                                                                     gCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAgcgc                       >  5:74635/1‑149 (MQ=255)
                                                                                                                       ggTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtg                     <  4:1141/149‑1 (MQ=255)
                                                                                                                        gTAATTGCAACGGGCCACGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                    <  8:189467/149‑1 (MQ=255)
                                                                                                                        gTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAgcgc                       <  4:181870/146‑1 (MQ=255)
                                                                                                                        gTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAgcgc                       >  3:181870/1‑146 (MQ=255)
                                                                                                                        gTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                    <  8:9510/149‑1 (MQ=255)
                                                                                                                        gTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                    <  8:8155/149‑1 (MQ=255)
                                                                                                                        gTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                    <  6:74635/149‑1 (MQ=255)
                                                                                                                        gTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                    <  6:19884/149‑1 (MQ=255)
                                                                                                                          gtacgcgaCGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                    <  6:165681/140‑1 (MQ=255)
                                                                                                                           aTGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtgatga                 <  8:177377/149‑1 (MQ=255)
                                                                                                                             ggCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATc               >  6:332334/1‑149 (MQ=255)
                                                                                                                             ggCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATc               >  5:187783/1‑149 (MQ=255)
                                                                                                                              gcgaCGCGCCCTGACCGGGGCGCGTTGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                    >  3:187029/4‑143 (MQ=255)
                                                                                                                              gcgaCGCGCCCTGACCGGGGCGCGTTGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                    <  4:187029/140‑1 (MQ=255)
                                                                                                                              gcgaCGCGCCCTGACCGGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCg              >  2:209048/4‑149 (MQ=255)
                                                                                                                                aaCGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTg                             <  2:179037/132‑1 (MQ=255)
                                                                                                                                aaCGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTg                             >  1:179037/1‑132 (MQ=255)
                                                                                                                                aaCGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGcc            >  7:240255/1‑149 (MQ=255)
                                                                                                                                  cgcgCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCg           >  2:300144/1‑148 (MQ=255)
                                                                                                                                   gcgcCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGaa                           >  6:98732/1‑131 (MQ=255)
                                                                                                                                      ccccGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTc      <  7:245716/149‑1 (MQ=255)
                                                                                                                                          ggTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGAGTGAAGCGCTGATGATCGCCGAACTCGATa  <  2:184183/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
CCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTCGATA  >  NC_002947/156246‑156532

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: