Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,453,480 A→G 8.8% intergenic (‑34/+79) PP_3069 ← / ← ppiC‑I outer membrane autotransporter/peptidyl‑prolyl cis‑trans isomerase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,453,4800AG 8.8% 183.0 / 7.4 68intergenic (‑34/+79)PP_3069/ppiC‑Iouter membrane autotransporter/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base A (27/35);  new base G (3/3);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

CCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACG  >  NC_002947/3453335‑3453627
                                                                                                                                                 |                                                                                                                                                   
ccgggccgggTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAATCCTGCGCACCACGCGTGGTGCACAg                                                                                                                                                  >  6:20477/1‑149 (MQ=255)
  gggccgggTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGTGCACCACGCGTGGTGCGCAggt                                                                                                                                                <  7:167055/149‑2 (MQ=255)
         gTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGTGCACCACGCGTGGTGCGCAGGTTTccctg                                                                                                                                         <  7:18971/149‑3 (MQ=255)
             gaAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGtt                                                                                                                                     >  1:176507/1‑149 (MQ=255)
             gaAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGtt                                                                                                                                     <  7:290913/149‑1 (MQ=255)
             gaAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGtt                                                                                                                                     >  5:129541/1‑149 (MQ=255)
             gaAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCAGGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTGCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGtt                                                                                                                                     >  1:268012/1‑149 (MQ=255)
              aaGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTg                                                                                                                                    >  4:208793/1‑149 (MQ=255)
                  tACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGtata                                                                                                                                <  3:133228/149‑1 (MQ=255)
                  tACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGtata                                                                                                                                <  3:138167/149‑1 (MQ=255)
                   aCAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGta                                                                                                                                  >  6:258481/1‑146 (MQ=255)
                   aCAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGAGCACAGGGTTACCGGGGTTGta                                                                                                                                  >  1:255166/1‑146 (MQ=255)
                     aaGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGt                                                                                                                                   >  2:107735/1‑143 (MQ=255)
                     aaGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGt                                                                                                                                   <  1:107735/143‑1 (MQ=255)
                     aaGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATagag                                                                                                                             <  8:109897/149‑1 (MQ=255)
                       gCCGCGTGGGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAgagga                                                                                                                           <  3:152683/149‑1 (MQ=255)
                            gtgtGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCg                                                                                                                      >  2:295387/1‑149 (MQ=255)
                                 cAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCAc                                                                                                                  <  6:261998/148‑1 (MQ=255)
                                 cAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCAc                                                                                                                  <  4:139200/148‑1 (MQ=255)
                                  aGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGa                                                                                                                <  6:129541/149‑1 (MQ=255)
                                     gggCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGAGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACca                                                                                                             >  5:65974/1‑149 (MQ=255)
                                               aTGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCa                                                                                                    >  3:214564/1‑148 (MQ=255)
                                                  cccAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCaaaaa                                                                                                >  8:39363/1‑149 (MQ=255)
                                                   ccAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAAt                                                                                               <  3:208793/149‑1 (MQ=255)
                                                     agcagcGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCCCGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGt                                                                                             >  7:9992/1‑149 (MQ=255)
                                                     agcagcGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGt                                                                                             >  1:98618/1‑149 (MQ=255)
                                                     agcagcGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGt                                                                                             <  8:218663/149‑1 (MQ=255)
                                                       cagcGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAATCCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCa                                                                                           <  5:20477/149‑1 (MQ=255)
                                                       cagcGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCa                                                                                           <  3:174684/149‑1 (MQ=255)
                                                         gcGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGAc                                                                                                               >  6:243005/1‑127 (MQ=255)
                                                         gcGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGAc                                                                                                               <  5:243005/127‑1 (MQ=255)
                                                               gCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCa                                                                                    >  3:199773/1‑148 (MQ=255)
                                                               gCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTAAGGTGGCa                                                                                    >  7:169584/1‑148 (MQ=255)
                                                                  aaCCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAgg                                                                                <  2:193285/149‑1 (MQ=255)
                                                                    cccAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGc                                                                              <  4:199773/149‑1 (MQ=255)
                                                                    cccAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGc                                                                              <  1:295387/149‑1 (MQ=255)
                                                                    cccAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGc                                                                              <  5:258481/149‑1 (MQ=255)
                                                                    cccAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGc                                                                              <  6:65974/149‑1 (MQ=255)
                                                                          gtgaGCATCAGGCGTAGATGCATCAGGAAAACCTCCATTCAGCAGGGAAAGCTGCGCAGCACGCGTGGTGCAGAGGGTTACTGGGGTTGTATATAGGAGTGCGTCTCGACCACAACTGGCAAAATTGTCAGGTGGGAAAGGGGTTgcgc                                                                        <  8:161639/146‑1 (MQ=255)
                                                                           tCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGcgcg                                                                       <  2:255166/149‑1 (MQ=255)
                                                                              gCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTa                                                                    <  6:27656/149‑1 (MQ=255)
                                                                                    ggCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTcgcg                                                              <  7:86384/149‑1 (MQ=255)
                                                                                      cGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGt                                                            <  2:98618/149‑1 (MQ=255)
                                                                                       gTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCgg                                                             >  4:24517/1‑147 (MQ=255)
                                                                                          gACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCa                                                                                           >  8:250642/1‑114 (MQ=255)
                                                                                          gACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGaa                                                        >  4:180489/1‑149 (MQ=255)
                                                                                              cATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGtt                                                    <  5:249257/149‑1 (MQ=255)
                                                                                               aTCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGttt                                                   <  2:268012/149‑1 (MQ=255)
                                                                                                tCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTg                                                  >  7:277188/1‑149 (MQ=255)
                                                                                                   ggAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCAc                                               >  6:258854/1‑149 (MQ=255)
                                                                                                   ggAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCAc                                               <  4:214564/149‑1 (MQ=255)
                                                                                                   ggAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCAc                                               <  7:39363/149‑1 (MQ=255)
                                                                                                          cTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtgg                                        >  3:305138/1‑149 (MQ=255)
                                                                                                           tCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtgg                                        >  8:253137/1‑148 (MQ=255)
                                                                                                            ccATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtgg                                        >  8:271835/1‑147 (MQ=255)
                                                                                                             cATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTa                                      <  5:120990/148‑1 (MQ=255)
                                                                                                                    cAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAAtt                              <  5:258854/149‑1 (MQ=255)
                                                                                                                       ggAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGc                           <  7:303238/149‑1 (MQ=255)
                                                                                                                          aaCCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGctct                        <  4:305138/149‑1 (MQ=255)
                                                                                                                          aaCCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGctct                        <  3:24517/149‑1 (MQ=255)
                                                                                                                            ccTGTGCACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAAtt                              >  2:216522/1‑141 (MQ=255)
                                                                                                                            ccTGTGCACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAAtt                              <  1:216522/141‑1 (MQ=255)
                                                                                                                            ccTGTGCACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTg                      >  5:28596/1‑149 (MQ=255)
                                                                                                                             cTGTGCACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGa                     >  2:276662/1‑149 (MQ=255)
                                                                                                                                     cacGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCg             >  8:188405/1‑149 (MQ=255)
                                                                                                                                       cgcgTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCgtg           >  3:148302/1‑149 (MQ=255)
                                                                                                                                         cgTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTa         <  8:9992/149‑1 (MQ=255)
                                                                                                                                          gttgtgCACAGGGTTACCGGGGTTGTATTGAGGAGTGCGTCACGACCTCAAGTGGCAAAAATGTCAGGTGGCAAACGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTTCTCTTGATCGGCCCGTGTAc        <  3:172193/146‑1 (MQ=255)
                                                                                                                                                 acaGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACg  >  6:51267/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACG  >  NC_002947/3453335‑3453627

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: ‑