Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 156,385 G→A 5.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 5.1% 563.2 / 11.9 177intergenic (‑10/+33)ygdQ/citNtransport protein/citrate transporter
Reads supporting (aligned to +/- strand):  ref base G (94/74);  new base A (4/5);  total (98/79)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 177 total aligned reads displayed.

GCGGCTGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTCGAT  >  NC_002947/156238‑156531
                                                                                                                                                   |                                                                                                                                                  
gCGGCTGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTc                                                                                                                                                  <  1:4304436/150‑1 (MQ=255)
 cGGCTGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCa                                                                                                                                                 <  1:1363489/150‑1 (MQ=255)
     tGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGgc                                                                                                                                             <  1:133192/150‑1 (MQ=255)
         cAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGcgcgt                                                                                                                                         >  2:862076/1‑150 (MQ=255)
         cAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGcgcgt                                                                                                                                         >  2:5276386/1‑150 (MQ=255)
         cAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGcgcgt                                                                                                                                         <  1:4738291/150‑1 (MQ=255)
         cAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGcgcgt                                                                                                                                         >  1:5449595/1‑150 (MQ=255)
         cAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGcgcgt                                                                                                                                         >  1:3250923/1‑150 (MQ=255)
          agcaTCGAGCTCATGATGATGTTGTCGATACCGAGGACGCTCTCACGCGCGGTAAGGGTGAAAAAGGCCACCCAGAGTTCCGGGCTGGTCAGCCATTCCAGGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCTGACCGGGGcgcgtt                                                                                                                                        <  2:1715153/147‑2 (MQ=255)
          aGACTCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCTGACCGGGGcgcgtt                                                                                                                                        <  2:3123514/150‑2 (MQ=255)
            aaTCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGcgcgtcgc                                                                                                                                      <  1:5276386/150‑1 (MQ=255)
               cGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTa                                                                                                                                   >  2:2793608/1‑150 (MQ=255)
                 agaTCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGt                                                                                                                                                   <  2:1597088/132‑1 (MQ=255)
                 agaTCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACg                                                                                                                                 >  1:1676211/1‑150 (MQ=255)
                 agaTCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACg                                                                                                                                 >  1:862324/1‑150 (MQ=255)
                    tCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGa                                                                                                                              <  2:2672409/150‑1 (MQ=255)
                    tCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGa                                                                                                                              <  1:3718447/150‑1 (MQ=255)
                    tCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGa                                                                                                                              >  2:4862199/1‑150 (MQ=255)
                       tgatgatgTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATAc                                                                                                                           >  2:4818366/1‑150 (MQ=255)
                          tgatgTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGa                                                                                                                        <  1:4247159/150‑1 (MQ=255)
                          tgatgTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGa                                                                                                                        >  1:2412633/1‑150 (MQ=255)
                          tgatgTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGa                                                                                                                        >  1:3756072/1‑150 (MQ=255)
                             tgttgtCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGt                                                                                                                     <  2:4530558/150‑1 (MQ=255)
                               ttgtCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGtct                                                                                                                   >  1:2272529/1‑148 (MQ=255)
                               ttgtCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAt                                                                                                                   >  2:4193403/1‑150 (MQ=255)
                                tgtCGATACCGAGGACGATCTCCAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTAtt                                                                                                                  <  1:3115463/150‑1 (MQ=255)
                                 gtCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTa                                                                                                                 >  1:138095/1‑150 (MQ=255)
                                  tCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTaa                                                                                                                >  1:4178294/1‑150 (MQ=255)
                                  tCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGAAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGGCGCGTACGTGATACGGAAGTATTaa                                                                                                                >  2:2281956/1‑150 (MQ=255)
                                   cGATACCGAGGACGATCTCCAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTaaa                                                                                                               <  1:1988422/150‑1 (MQ=255)
                                   cGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTaaa                                                                                                               <  1:3528340/150‑1 (MQ=255)
                                   cGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTaaa                                                                                                               <  1:1016126/150‑1 (MQ=255)
                                     aTACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGa                                                                                                             >  2:4002608/1‑150 (MQ=255)
                                         cGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGactac                                                                                                         >  2:4956532/1‑150 (MQ=255)
                                         cGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGactac                                                                                                         >  2:3637338/1‑150 (MQ=255)
                                                  tctcAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCg                                                                                                <  1:975860/150‑1 (MQ=255)
                                                     cAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGa                                                                                             >  2:1070366/1‑150 (MQ=255)
                                                            ggTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTAATTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATccc                                                                                      <  2:2224620/150‑1 (MQ=255)
                                                             gTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATcccc                                                                                     <  1:3612088/150‑1 (MQ=255)
                                                                aGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATc                                                                                  >  1:62195/1‑150 (MQ=255)
                                                                 gggTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATca                                                                                 <  2:558765/150‑1 (MQ=255)
                                                                 gggTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATca                                                                                 <  2:2326744/150‑1 (MQ=255)
                                                                   gTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATcagc                                                                               >  2:1053121/1‑150 (MQ=255)
                                                                    tGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATcagca                                                                              >  2:4599231/1‑150 (MQ=255)
                                                                        aaaGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTgc                                                                          <  1:1886112/150‑1 (MQ=255)
                                                                             cAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGa                                                                     <  2:2939471/150‑1 (MQ=255)
                                                                             cAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGa                                                                     <  1:3303043/150‑1 (MQ=255)
                                                                              aaCCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGaaa                                                                                           <  2:4645713/127‑1 (MQ=255)
                                                                              aaCCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAg                                                                    <  1:3712231/150‑1 (MQ=255)
                                                                               aCCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGc                                                                   >  2:28877/1‑150 (MQ=255)
                                                                                cccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCa                                                                  <  2:138095/150‑1 (MQ=255)
                                                                                 ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCAt                                                                 >  1:2889374/1‑150 (MQ=255)
                                                                                 ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCAt                                                                 >  1:1968788/1‑150 (MQ=255)
                                                                                  cAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCAtt                                                                <  1:754624/150‑1 (MQ=255)
                                                                                  cAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCAtt                                                                >  1:2798993/1‑150 (MQ=255)
                                                                                   aGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTa                                                               >  1:1215915/1‑150 (MQ=255)
                                                                                         ccGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGCGCGCCGAGCATTATGcaca                                                         >  2:181914/1‑150 (MQ=255)
                                                                                                 gTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGAGCGCCGAGCATTATGCACAGGCATAcc                                                 >  2:4561187/1‑150 (MQ=255)
                                                                                                  tCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCa                                                >  2:1101754/1‑150 (MQ=255)
                                                                                                     gCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCa                                             <  1:4956532/150‑1 (MQ=255)
                                                                                                         ttCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATca                                                                                 <  2:1047969/110‑1 (MQ=255)
                                                                                                         ttCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATca                                                                                 >  1:1047969/1‑110 (MQ=255)
                                                                                                         ttCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGc                                         >  2:2889306/1‑150 (MQ=255)
                                                                                                         ttCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGc                                         <  1:1631976/150‑1 (MQ=255)
                                                                                                             aTGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCAc                                     >  1:3000747/1‑150 (MQ=255)
                                                                                                               gTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACtt                                   <  2:2798993/150‑1 (MQ=255)
                                                                                                                tATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTg                                  >  1:212056/1‑150 (MQ=255)
                                                                                                                tATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTg                                  >  2:701566/1‑150 (MQ=255)
                                                                                                                 aTTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGa                                 >  2:3580619/1‑150 (MQ=255)
                                                                                                                 aTTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGa                                 <  2:4682203/150‑1 (MQ=255)
                                                                                                                  ttttCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAg                                >  2:2033109/1‑150 (MQ=255)
                                                                                                                  ttttCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAg                                <  2:338111/150‑1 (MQ=255)
                                                                                                                   tttCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAgg                                                       <  1:2099327/126‑1 (MQ=255)
                                                                                                                   tttCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGc                               >  2:3103967/1‑150 (MQ=255)
                                                                                                                    ttCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACtt                                   >  2:1828540/1‑145 (MQ=255)
                                                                                                                        ttCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGaa                          >  2:4074638/1‑150 (MQ=255)
                                                                                                                        ttCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGaa                          >  1:3083905/1‑150 (MQ=255)
                                                                                                                        ttCGAGCGGTAATGGCAACGCGCCCCGGTCACGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCAACAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGaa                          >  1:3833327/1‑150 (MQ=255)
                                                                                                                              cGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtg                    <  1:2310012/150‑1 (MQ=255)
                                                                                                                               ggTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                   <  1:4091120/150‑1 (MQ=255)
                                                                                                                               ggTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                   <  1:4818366/150‑1 (MQ=255)
                                                                                                                               ggTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                   <  2:212056/150‑1 (MQ=255)
                                                                                                                                 tAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtgatt                 >  1:447300/1‑149 (MQ=255)
                                                                                                                                    tGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATc              >  1:4738190/1‑150 (MQ=255)
                                                                                                                                    tGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATc              >  1:5252943/1‑150 (MQ=255)
                                                                                                                                    tGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATc              <  1:28877/150‑1 (MQ=255)
                                                                                                                                     ggCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCg             >  2:3485458/1‑150 (MQ=255)
                                                                                                                                      gcgaCGCGCCCTGACCGGGGCGCGTTGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATcagca                                                                              >  2:2612016/4‑84 (MQ=37)
                                                                                                                                      gcgaCGCGCCCTGACCGGGGCGCGTTGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATcagca                                                                              <  1:2612016/81‑1 (MQ=37)
                                                                                                                                      gCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGc            <  2:2308562/150‑1 (MQ=255)
                                                                                                                                      gCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGc            >  2:517276/1‑150 (MQ=255)
                                                                                                                                       cAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGcc           <  2:3970803/150‑1 (MQ=255)
                                                                                                                                       cAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGcc           >  2:2810024/1‑150 (MQ=255)
                                                                                                                                        gaCGCGCCCTGACCGGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCg          >  1:1432144/2‑150 (MQ=255)
                                                                                                                                        aaCGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCg          >  2:561625/1‑150 (MQ=255)
                                                                                                                                         aCGCGCCCTGACCGGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCa                                                >  2:1517202/1‑111 (MQ=255)
                                                                                                                                         aCGCGCCCTGACCGGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCa                                                <  1:1517202/111‑1 (MQ=255)
                                                                                                                                           gcgcCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAAc       <  1:833215/150‑1 (MQ=255)
                                                                                                                                             gcCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTc     <  2:1641936/150‑1 (MQ=255)
                                                                                                                                             gcCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTc     >  1:3565725/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCGGCTGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTCGAT  >  NC_002947/156238‑156531

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: