Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 156,383 C→T 6.1% intergenic (‑8/+35) ygdQ ← / ← citN transport protein/citrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947156,3830CT 6.1% 242.4 / 4.2 81intergenic (‑8/+35)ygdQ/citNtransport protein/citrate transporter
Reads supporting (aligned to +/- strand):  ref base C (43/33);  new base T (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 = 7.91e-01

CATGCGGCTGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTC  >  NC_002947/156235‑156528
                                                                                                                                                    |                                                                                                                                                   
cATGCGGCTGACCAGAATCGAGATCATGATGTTGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGcccc                                                                                                                                                     >  1:127143/1‑149 (MQ=255)
 aTGCGGCTGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCg                                                                                                                                                    <  2:127143/149‑1 (MQ=255)
           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                                                                    >  5:142656/1‑147 (MQ=255)
                                                                                 aaCCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGcc                                                                    >  6:134751/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                                                               >  7:301773/1‑149 (MQ=255)
                                                                                    ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCa                                                               >  7:9510/1‑149 (MQ=255)
                                                                                    ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCa                                                               >  8:42882/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                                   >  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                                >  3:328570/1‑148 (MQ=255)
                                                                                                                    aTTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACtt                                >  5:19451/1‑148 (MQ=255)
                                                                                                                    aTTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACtt                                >  3:236961/1‑148 (MQ=255)
                                                                                                                    aTTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACtt                                >  7:177377/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                   >  5:74635/1‑149 (MQ=255)
                                                                                                                                gCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAgcgc                   <  7:42882/149‑1 (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:19884/149‑1 (MQ=255)
                                                                                                                                   gTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                <  6:74635/149‑1 (MQ=255)
                                                                                                                                     gtacgcgaCGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                <  6:165681/140‑1 (MQ=255)
                                                                                                                                      aTGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtgatga             <  8:177377/149‑1 (MQ=255)
                                                                                                                                        ggCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATc           >  5:187783/1‑149 (MQ=255)
                                                                                                                                        ggCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATc           >  6:332334/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)
                                                                                                                                                    |                                                                                                                                                   
CATGCGGCTGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAT‑‑TAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTC  >  NC_002947/156235‑156528

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: