Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 156,386 T→C 5.1% intergenic (‑11/+32) ygdQ ← / ← citN transport protein/citrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947156,3860TC 5.1% 560.2 / 11.5 175intergenic (‑11/+32)ygdQ/citNtransport protein/citrate transporter
Reads supporting (aligned to +/- strand):  ref base T (93/73);  new base C (4/5);  total (97/78)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

Only 100 of 178 total aligned reads displayed.

GCGGCTGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTCGATACCG  >  NC_002947/156238‑156535
                                                                                                                                                    |                                                                                                                                                     
 cGGCTGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCa                                                                                                                                                     <  1:1363489/150‑1 (MQ=255)
     tGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGgc                                                                                                                                                 <  1:133192/150‑1 (MQ=255)
       aCCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGgcgc                                                                                                                                               >  2:1016126/1‑150 (MQ=255)
         cAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGcgcgt                                                                                                                                             >  2:862076/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)
         cAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGcgcgt                                                                                                                                             >  2:5276386/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)
             aTCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGcgcgtcgcg                                                                                                                                         >  1:672120/1‑150 (MQ=255)
               cGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTa                                                                                                                                       >  2:2793608/1‑150 (MQ=255)
                 agaTCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGt                                                                                                                                                       <  2:1597088/132‑1 (MQ=255)
                 agaTCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACg                                                                                                                                     >  1:862324/1‑150 (MQ=255)
                 agaTCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACg                                                                                                                                     >  1:3819858/1‑150 (MQ=255)
                   aTCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTg                                                                                                                                   >  2:3202050/1‑150 (MQ=255)
                    tCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGa                                                                                                                                  <  1:3718447/150‑1 (MQ=255)
                    tCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGa                                                                                                                                  <  1:3604318/150‑1 (MQ=255)
                    tCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGa                                                                                                                                  >  2:4862199/1‑150 (MQ=255)
                    tCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGa                                                                                                                                  <  2:2672119/150‑1 (MQ=255)
                          tgatgTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGa                                                                                                                            <  2:2151622/150‑1 (MQ=255)
                          tgatgTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGa                                                                                                                            >  1:2412633/1‑150 (MQ=255)
                          tgatgTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGa                                                                                                                            <  1:4247159/150‑1 (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)
                                  tCGATACCGAGGACGATCTCAAGCGCTGTAAGGGTTAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTaa                                                                                                                    <  1:560139/150‑1 (MQ=255)
                                  tCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTaa                                                                                                                    >  1:4178294/1‑150 (MQ=255)
                                   cGATACCGAGGACGATCTCCAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTaaa                                                                                                                   <  1:1988422/150‑1 (MQ=255)
                                   cGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTaaa                                                                                                                   <  1:1016126/150‑1 (MQ=255)
                                     aTACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGa                                                                                                                 >  2:4002608/1‑150 (MQ=255)
                                         cGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCTAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGactac                                                                                                             >  1:5190265/1‑150 (MQ=255)
                                         cGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGactac                                                                                                             >  2:4956532/1‑150 (MQ=255)
                                         cGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGactac                                                                                                             >  2:3637338/1‑150 (MQ=255)
                                         cGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGactac                                                                                                             >  1:428389/1‑150 (MQ=255)
                                                  tctcAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCg                                                                                                    <  1:975860/150‑1 (MQ=255)
                                                     cAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGa                                                                                                 >  2:1070366/1‑150 (MQ=255)
                                                          gcgGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATc                                                                                            <  2:1749814/150‑1 (MQ=255)
                                                                 gggTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATca                                                                                     <  2:2326744/150‑1 (MQ=255)
                                                                 gggTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATca                                                                                     <  2:558765/150‑1 (MQ=255)
                                                                 gggTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATca                                                                                     <  2:4792486/150‑1 (MQ=255)
                                                                    tGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATcagca                                                                                  <  2:392273/150‑1 (MQ=255)
                                                                    tGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATcagca                                                                                  >  2:4599231/1‑150 (MQ=255)
                                                                        aaaGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTgc                                                                              <  1:4326076/150‑1 (MQ=255)
                                                                        aaaGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTgc                                                                              <  1:1886112/150‑1 (MQ=255)
                                                                             cAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGa                                                                         <  1:3303043/150‑1 (MQ=255)
                                                                             cAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGa                                                                         <  2:2939471/150‑1 (MQ=255)
                                                                              aaCCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGaaa                                                                                               <  2:4645713/127‑1 (MQ=255)
                                                                              aaCCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGaaa                                                                                               >  1:4645713/1‑127 (MQ=255)
                                                                                 ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCAt                                                                     >  1:1968788/1‑150 (MQ=255)
                                                                                 ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCAt                                                                     >  1:2889374/1‑150 (MQ=255)
                                                                                 ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCAt                                                                     >  1:4750101/1‑150 (MQ=255)
                                                                                 ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCAt                                                                     >  2:2279979/1‑150 (MQ=255)
                                                                                  cAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCAtt                                                                    >  2:1515646/1‑150 (MQ=255)
                                                                                  cAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCAtt                                                                    >  1:2798993/1‑150 (MQ=255)
                                                                                   aGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTa                                                                   >  1:1215915/1‑150 (MQ=255)
                                                                                          cGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAg                                                            <  1:3920432/150‑1 (MQ=255)
                                                                                               tGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGa                                                                                                                                  >  2:1226126/1‑75 (MQ=255)
                                                                                                 gTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATAcc                                                     >  1:2737467/1‑150 (MQ=255)
                                                                                                 gTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGAGCGCCGAGCATTATGCACAGGCATAcc                                                     >  2:4561187/1‑150 (MQ=255)
                                                                                                  tCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCa                                                    >  2:1101754/1‑150 (MQ=255)
                                                                                                         ttCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATca                                                                                     >  1:1047969/1‑110 (MQ=255)
                                                                                                         ttCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGc                                             <  1:1631976/150‑1 (MQ=255)
                                                                                                         ttCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGc                                             >  2:2889306/1‑150 (MQ=255)
                                                                                                               gTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACtt                                       <  2:2798993/150‑1 (MQ=255)
                                                                                                                tATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTg                                      >  2:701566/1‑150 (MQ=255)
                                                                                                                tATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTg                                      >  1:212056/1‑150 (MQ=255)
                                                                                                                 aTTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGa                                     <  2:4682203/150‑1 (MQ=255)
                                                                                                                 aTTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGa                                     <  2:105247/150‑1 (MQ=255)
                                                                                                                 aTTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGa                                     >  2:3580619/1‑150 (MQ=255)
                                                                                                                  ttttCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAg                                    >  2:2033109/1‑150 (MQ=255)
                                                                                                                   tttCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAgg                                                           <  1:2099327/126‑1 (MQ=255)
                                                                                                                   tttCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGc                                   >  2:3103967/1‑150 (MQ=255)
                                                                                                                    ttCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACtt                                       <  1:1828540/145‑1 (MQ=255)
                                                                                                                    ttCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACtt                                       >  2:1828540/1‑145 (MQ=255)
                                                                                                                       cTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGa                               <  1:3643915/150‑1 (MQ=255)
                                                                                                                        ttCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGaa                              >  2:4074638/1‑150 (MQ=255)
                                                                                                                        ttCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGaa                              >  1:3083905/1‑150 (MQ=255)
                                                                                                                              cGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtg                        <  1:2310012/150‑1 (MQ=255)
                                                                                                                               ggTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                       <  1:3041148/150‑1 (MQ=255)
                                                                                                                               ggTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                       <  2:212056/150‑1 (MQ=255)
                                                                                                                                  gtacgCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGaa                              <  1:3547827/136‑1 (MQ=255)
                                                                                                                                    tGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATc                  <  1:28877/150‑1 (MQ=255)
                                                                                                                                     ggCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCg                 >  2:3485458/1‑150 (MQ=255)
                                                                                                                                      gcgaCGCGCCCTGACCGGGGCGCGTTGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATcagca                                                                                  >  2:2612016/4‑84 (MQ=37)
                                                                                                                                      gCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGc                <  2:2308562/150‑1 (MQ=255)
                                                                                                                                      gCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGc                >  2:517276/1‑150 (MQ=255)
                                                                                                                                       cAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGcc               >  2:2810024/1‑150 (MQ=255)
                                                                                                                                       cAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGcc               >  2:2810258/1‑150 (MQ=255)
                                                                                                                                        gaCGCGCCCTGACCGGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCg              >  1:1432144/2‑150 (MQ=255)
                                                                                                                                        aaCGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCg              >  2:561625/1‑150 (MQ=255)
                                                                                                                                         aCGCGCCCTGACCGGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCa                                                    <  1:1517202/111‑1 (MQ=255)
                                                                                                                                           gcgcCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAAc           <  1:833215/150‑1 (MQ=255)
                                                                                                                                             gcCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTc         <  2:1641936/150‑1 (MQ=255)
                                                                                                                                              ccccGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTCg        >  1:4865198/1‑150 (MQ=255)
                                                                                                                                               cccGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTCGa       >  1:3746288/1‑150 (MQ=255)
                                                                                                                                                    tCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTCGATACCg  <  1:3555600/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGGCTGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTCGATACCG  >  NC_002947/156238‑156535

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: