Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 156,388 A→G 5.9% intergenic (‑13/+30) ygdQ ← / ← citN transport protein/citrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947156,3880AG 5.9% 249.6 / 3.2 83intergenic (‑13/+30)ygdQ/citNtransport protein/citrate transporter
Reads supporting (aligned to +/- strand):  ref base A (43/35);  new base G (3/2);  total (46/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 = 1.00e+00

CCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTCGATACC  >  NC_002947/156246‑156534
                                                                                                                                              |                                                                                                                                                    
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                                                                                                                                       <  2:86614/139‑1 (MQ=255)
                   gatgTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTAc                                                                                                                                       >  1:86614/1‑139 (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                                                                     >  5:86279/1‑149 (MQ=255)
                                                                         ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCa                                                                     >  7:8155/1‑149 (MQ=255)
                                                                         ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCa                                                                     >  8:42882/1‑149 (MQ=255)
                                                                         ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCa                                                                     >  8:110538/1‑149 (MQ=255)
                                                                         ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCa                                                                     >  7:9510/1‑149 (MQ=255)
                                                                         ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCa                                                                     >  7:301773/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                                         >  8:245716/1‑149 (MQ=255)
                                                                                                     aTGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCa                                         <  7:168173/149‑1 (MQ=255)
                                                                                                      tGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCtggtatgcctgtgcataatgctcggcgcactgctgat                                        <  4:234861‑M1/149‑38 (MQ=255)
                                                                                                      tGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGa                                                                        <  4:36358/117‑1 (MQ=255)
                                                                                                         aTTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACtt                                      >  5:19451/1‑148 (MQ=255)
                                                                                                         aTTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACtt                                      >  3:328570/1‑148 (MQ=255)
                                                                                                         aTTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACtt                                      >  7:177377/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                      <  4:187029/140‑1 (MQ=255)
                                                                                                                              gcgaCGCGCCCTGACCGGGGCGCGTTGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                      >  3:187029/4‑143 (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)
                                                                                                                                            tCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTCGATAcc  <  4:328570/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTCGATACC  >  NC_002947/156246‑156534

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: