Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,453,506 C→T 9.8% intergenic (‑17/+96) 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,5060CT 9.8% 238.9 / 12.8 92intergenic (‑17/+96)PP_3069/ppiC‑Iouter membrane autotransporter/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base C (38/45);  new base T (4/5);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

ACAGGTGTAGGTGTCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGAT  >  NC_002947/3453358‑3453652
                                                                                                                                                    |                                                                                                                                                  
acaGGTGTAGGTGTCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTgcg                                                                                                                                                   <  1:340642/150‑1 (MQ=255)
        aggtgtCGTTACCCGGGCCGGGTGAGAATGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGTGCACCAcgc                                                                                                                                           <  1:395667/150‑1 (MQ=255)
         ggtgtCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCCGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACTAGGGAAACCTGTGCACCAcgcg                                                                                                                                          >  2:914983/1‑150 (MQ=255)
          gtgtCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCgt                                                                                                                                         >  2:346819/1‑150 (MQ=255)
          gtgtCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCgt                                                                                                                                         >  1:1210359/1‑150 (MQ=255)
          gtgtCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCgt                                                                                                                                         >  1:982500/1‑150 (MQ=255)
          gtgtCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCgt                                                                                                                                         >  1:669896/1‑150 (MQ=255)
           tgtCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCgtg                                                                                                                                        <  1:188772/150‑1 (MQ=255)
           tgtCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCgtg                                                                                                                                        <  1:84261/150‑1 (MQ=255)
           ggtCGTTACCCGGGCCGGGTGGGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCgtg                                                                                                                                        <  2:964392/149‑1 (MQ=255)
           ggtCGTTACCCGGGCCGGGTGAGCAGGTACAAGCCGCGTGCGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGCGCAGCATCAGGCGTAGACGCTTCAGGAAAACCTCCATTCTCCCGGGAAACCTGCGCACCACGCgtg                                                                                                                                        <  2:877857/149‑1 (MQ=255)
            gtCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGca                                                                                                                                                 >  1:841191/1‑140 (MQ=255)
            gtCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGca                                                                                                                                                 <  2:841191/140‑1 (MQ=255)
             tCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGGCCAGCAGCGGGAGCGAACGCAGGGTCAGCATCAGGCGTAGACGCATCAGGAACACCTCCATTCACCAGGTAACCCCGCGCACAACGCgtggt                                                                                                                                      <  2:176576/150‑1 (MQ=255)
             tCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCgtggt                                                                                                                                      <  1:111517/150‑1 (MQ=255)
             tCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCgtggt                                                                                                                                      <  2:1172285/150‑1 (MQ=255)
             tCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCgtggg                                                                                                                                      >  1:1130565/1‑149 (MQ=255)
             tCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCgt                                                                                                                                         <  1:440943/147‑1 (MQ=255)
             tCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCgt                                                                                                                                         >  2:440943/1‑147 (MQ=255)
               gTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGAGc                                                                                                                                    >  2:582699/1‑150 (MQ=255)
                  aCCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCCGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACTAGGGAAACCTGTGCACCACGCGTGGTGcgca                                                                                                                                 <  1:914983/150‑1 (MQ=255)
                   cccGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCTAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAg                                                                                                                                >  1:1642919/1‑150 (MQ=255)
                    ccgggccgggTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGGCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAgg                                                                                                                               >  2:1069980/1‑150 (MQ=255)
                     cgggccgggTGAGAAGGTACATGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCAGCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGTGCACCACGCGTGGTGCGCAggt                                                                                                                              <  1:78153/150‑2 (MQ=255)
                     cgggccgggTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAggg                                                                                                                              <  2:624490/150‑1 (MQ=255)
                      gggccgggTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGcgca                                                                                                                                 >  2:818210/1‑143 (MQ=255)
                      gggccgggTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGcgca                                                                                                                                 <  1:818210/146‑1 (MQ=255)
                          cgggTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTAcc                                                                                                                         >  2:811980/1‑150 (MQ=255)
                            ggTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCgg                                                                                                                       >  1:515873/1‑150 (MQ=255)
                              tGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCgggg                                                                                                                     <  2:982500/150‑1 (MQ=255)
                                     gTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGtata                                                                                                              <  1:346819/150‑1 (MQ=255)
                                     gTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGtata                                                                                                              <  1:614751/150‑1 (MQ=255)
                                     gTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGtata                                                                                                              >  2:692351/1‑150 (MQ=255)
                                     gTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGtata                                                                                                              >  2:614751/1‑150 (MQ=255)
                                       aCAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATaga                                                                                                            >  1:1086112/1‑150 (MQ=255)
                                        cAAGCCGCGTGTGCAGCGGGCGAGCAAATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATngag                                                                                                           <  2:1186360/150‑1 (MQ=255)
                                          aGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAgagga                                                                                                         <  1:692351/150‑1 (MQ=255)
                                          aGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAgagga                                                                                                         <  2:1536811/150‑1 (MQ=255)
                                           gCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGta                                                                                                                >  2:216990/1‑142 (MQ=255)
                                           gCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGta                                                                                                                <  1:216990/142‑1 (MQ=255)
                                           gCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAgaggag                                                                                                        >  2:1232861/1‑150 (MQ=255)
                                                     cAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGTGCACCAcgc                                                                                                                                           <  1:740814/105‑1 (MQ=255)
                                                     cAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGTGCACCAcgc                                                                                                                                           >  2:740814/1‑105 (MQ=255)
                                                      aGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGAc                                                                                             <  2:1086112/150‑1 (MQ=255)
                                                         gggCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACcac                                                                                          <  2:1210359/150‑1 (MQ=255)
                                                              aGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTg                                                                                     <  1:279788/150‑1 (MQ=255)
                                                               ggCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTgg                                                                                    <  2:515873/150‑1 (MQ=255)
                                                               ggCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTgg                                                                                    >  2:517946/1‑150 (MQ=255)
                                                                      cccAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGACACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAAt                                                                             >  1:230258/1‑150 (MQ=255)
                                                                             gcGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGt                                                                      <  2:515050/150‑1 (MQ=255)
                                                                             gcGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGt                                                                      <  2:230258/150‑1 (MQ=255)
                                                                             gcGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCACCAGGATATCCTCCAGTCACCAGGGAAACCTGCGCCACACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGt                                                                      <  1:1069980/150‑1 (MQ=255)
                                                                             gcGGGAGCGAACCAAGGGTCAGCATCAGGCGTAGGCGCATCAGGAAAGCGTTTATGCAACAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGCTGTATAGAGGAGTGCGCCACGCGCAAAACTGGCAAAAATGTCAGGt                                                                      >  2:741328/1‑150 (MQ=255)
                                                                                 gAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCa                                                                  >  1:298335/1‑150 (MQ=255)
                                                                                     gAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAgg                                                              >  1:30028/1‑150 (MQ=255)
                                                                                       aCCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAgg                                                              >  2:754522/1‑148 (MQ=255)
                                                                                       aCCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAgg                                                              <  1:754522/148‑1 (MQ=255)
                                                                                          cAGGGTCAGCATCAGGCGTAGACACATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCTACTATGTCAGGTGGCAAAGGGCTAg                                                         >  2:950808/1‑150 (MQ=255)
                                                                                           aGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAgc                                                        >  1:297906/1‑150 (MQ=255)
                                                                                               tCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGt                                                    <  1:1232861/150‑1 (MQ=255)
                                                                                                 aGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTa                                                  <  1:950808/150‑1 (MQ=255)
                                                                                                 aGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTa                                                  <  2:298335/150‑1 (MQ=255)
                                                                                                      cAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCAAGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTcgc                                             <  2:1566694/150‑1 (MQ=255)
                                                                                                       aGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTcgcg                                            >  1:1147015/1‑150 (MQ=255)
                                                                                                        ggCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCCGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCgg                                           <  2:603600/150‑1 (MQ=255)
                                                                                                           gTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTAACAAATTTGTAGCGCGTTAGTCGCGGTAg                                        >  2:441384/1‑150 (MQ=255)
                                                                                                                 gCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCgg                                           >  1:502073/1‑141 (MQ=255)
                                                                                                                 gCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCgg                                           <  2:502073/141‑1 (MQ=255)
                                                                                                                  gaTCAGGAAATCCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCGAGCGCGTTAGTCGCGGTAGAACGttt                                 <  2:1176196/149‑1 (MQ=255)
                                                                                                                  cATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGATCACAACTGGCAAAAATGTCAGGTGGCAAAGGTCTAGCGCGTTAGTCGCGGTAGAACGttt                                 <  1:663962/150‑1 (MQ=255)
                                                                                                                  cATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGttt                                 <  2:297906/150‑1 (MQ=255)
                                                                                                                  cATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGttt                                 <  2:763004/150‑1 (MQ=255)
                                                                                                                        gAAAACCTCCATTCACCAGGGAAACCTGCGCACCAcgcg                                                                                                                                          >  2:994067/1‑39 (MQ=255)
                                                                                                                        gAAAACCTCCATTCACCAGGGAAACCTGCGCACCAcgcg                                                                                                                                          <  1:994067/39‑1 (MQ=255)
                                                                                                                         aaaaCCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAg                          <  2:669896/150‑1 (MQ=255)
                                                                                                                            aCCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtg                       >  2:294978/1‑150 (MQ=255)
                                                                                                                            aCCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtg                       <  2:437217/150‑1 (MQ=255)
                                                                                                                              cTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTCGCGCGTTAGTCCCTGTAGAACGTTTGCACCAggtggt                     >  1:479781/1‑150 (MQ=255)
                                                                                                                              cTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtggt                     >  2:165847/1‑150 (MQ=255)
                                                                                                                              cTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtggt                     >  2:124598/1‑150 (MQ=255)
                                                                                                                               tCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTa                    <  2:170596/150‑1 (MQ=255)
                                                                                                                                 cATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTcgc                                             >  2:1690502/1‑123 (MQ=255)
                                                                                                                                 cATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTcgc                                             <  1:1690502/123‑1 (MQ=255)
                                                                                                                                         aGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTg          >  1:1312379/1‑150 (MQ=255)
                                                                                                                                          gggAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGc         <  1:811980/150‑1 (MQ=255)
                                                                                                                                          gggAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGc         <  1:1894853/150‑1 (MQ=255)
                                                                                                                                            gAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGctc       <  2:1147015/150‑1 (MQ=255)
                                                                                                                                            gAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGctc       <  1:1412200/150‑1 (MQ=255)
                                                                                                                                               aCCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTg    >  2:1803758/1‑150 (MQ=255)
                                                                                                                                                ccTGTGCACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGttt                                 <  2:923387/120‑1 (MQ=255)
                                                                                                                                                ccTGTGCACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGttt                                 >  1:923387/1‑120 (MQ=255)
                                                                                                                                                 cTGTGCACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGAt  >  2:91778/1‑150 (MQ=255)
                                                                                                                                                 cTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGAt  >  2:210761/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACAGGTGTAGGTGTCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGAT  >  NC_002947/3453358‑3453652

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