Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,455,445 A→G 6.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,455,4450AG 6.8% 677.7 / 24.6 219intergenic (‑34/+79)PP_3069/ppiC‑Iouter membrane autotransporter/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base A (104/100);  new base G (10/5);  total (114/105)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 235 total aligned reads displayed.

ACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCA  >  NC_002947/3455298‑3455594
                                                                                                                                                   |                                                                                                                                                     
aCCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGTGCACCACGCGTGGTGcgca                                                                                                                                                     <  2:2306726/150‑1 (MQ=255)
aCCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGcaca                                                                                                                                                     <  2:1161378/150‑1 (MQ=255)
 cccGGGCCGGGTGAGAAGGGACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAg                                                                                                                                                    <  2:5011566/150‑1 (MQ=255)
  ccgggccgggTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAgg                                                                                                                                                   >  2:4004557/1‑150 (MQ=255)
    gggccgggTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGt                                                                                                                                                 >  1:442154/1‑150 (MQ=255)
        cgggTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGTAACCCTGTGCACCACGCGTGGTGCGCAGGTTTccc                                                                                                                                             <  2:2470226/150‑1 (MQ=39)
            tGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCgggg                                                                                                                                         <  2:9049248/150‑1 (MQ=255)
            tGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCgggg                                                                                                                                         <  2:2180979/150‑1 (MQ=255)
            tGAGAAGGAACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGCGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCgggg                                                                                                                                         <  2:788556/150‑1 (MQ=255)
               gaAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTg                                                                                                                                      <  1:2296998/150‑1 (MQ=255)
                aaGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGt                                                                                                                                     >  2:3091590/1‑150 (MQ=255)
                   gTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGtata                                                                                                                                  <  1:2749768/150‑1 (MQ=255)
                   gTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGtata                                                                                                                                  <  1:6471239/150‑1 (MQ=255)
                      cAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATagag                                                                                                                               <  1:872937/150‑1 (MQ=255)
                      cAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATagag                                                                                                                               <  1:2061097/150‑1 (MQ=255)
                         gCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAgaggag                                                                                                                            <  1:3926804/150‑1 (MQ=255)
                                   cAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGTGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGa                                                                                                                  <  1:4337613/150‑1 (MQ=255)
                                      cGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACca                                                                                                               >  1:8232673/1‑150 (MQ=255)
                                         gCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACaa                                                                                                            <  2:4297962/150‑1 (MQ=255)
                                           gAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACt                                                                                                          >  1:3054953/1‑150 (MQ=255)
                                           gAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACt                                                                                                          >  2:6129369/1‑150 (MQ=255)
                                            aGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTg                                                                                                         <  2:8027098/150‑1 (MQ=255)
                                                          agcGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAgg                                                                                           <  2:6693537/150‑1 (MQ=255)
                                                           gcGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGt                                                                                          <  1:5453879/150‑1 (MQ=255)
                                                                  cGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAg                                                                                   <  2:2147266/150‑1 (MQ=255)
                                                                   gAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAgg                                                                                  <  1:4594688/150‑1 (MQ=255)
                                                                       ccAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCtc                                                                              >  1:1170151/1‑149 (MQ=255)
                                                                       ccAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTa                                                                              >  1:8320064/1‑150 (MQ=255)
                                                                       ccAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTa                                                                              >  1:1156876/1‑150 (MQ=255)
                                                                       ccAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTa                                                                              >  1:7662289/1‑150 (MQ=255)
                                                                         aGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAgc                                                                            >  1:324723/1‑150 (MQ=255)
                                                                             tCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGt                                                                        <  2:6578500/150‑1 (MQ=255)
                                                                             tCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGt                                                                        >  2:845702/1‑150 (MQ=255)
                                                                               aGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTa                                                                      <  2:3821243/150‑1 (MQ=255)
                                                                                  aTCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTc                                                                   >  2:6263170/1‑150 (MQ=255)
                                                                                  aTCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGGGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGGTAGTc                                                                   >  1:798117/1‑150 (MQ=255)
                                                                                     aGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTcgcg                                                                >  2:3391932/1‑150 (MQ=255)
                                                                                      ggCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCgg                                                               <  2:4541848/150‑1 (MQ=255)
                                                                                       gCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGt                                                              >  2:6328276/1‑150 (MQ=255)
                                                                                        cGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTa                                                             <  1:2920926/150‑1 (MQ=255)
                                                                                         gTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAg                                                            >  1:850797/1‑150 (MQ=255)
                                                                                               gCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGtt                                                      >  1:2055359/1‑150 (MQ=255)
                                                                                               gCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGtt                                                      <  1:4765875/150‑1 (MQ=255)
                                                                                               gCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGtt                                                      >  1:6697327/1‑150 (MQ=255)
                                                                                                cATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGttt                                                     <  1:336943/150‑1 (MQ=255)
                                                                                                 aTCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTg                                                    <  2:4497616/150‑1 (MQ=255)
                                                                                                     ggAAAACCTCCATTCACCAGTGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCAcc                                                >  1:7117117/1‑150 (MQ=255)
                                                                                                        aaaCCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAgg                                             <  1:6347064/150‑1 (MQ=255)
                                                                                                            cTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAATTGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtggt                                         >  2:5295182/1‑150 (MQ=255)
                                                                                                            cTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtggt                                         >  1:6397302/1‑150 (MQ=255)
                                                                                                            cTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtggt                                         >  1:3508496/1‑150 (MQ=255)
                                                                                                            cTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtggt                                         >  2:1348996/1‑150 (MQ=255)
                                                                                                            cTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGGTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtggt                                         >  2:3456369/1‑150 (MQ=255)
                                                                                                            cTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTAACCGGGGTTGTATAGAGGAGAGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCAAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAggtgga                                         >  2:2935275/1‑149 (MQ=255)
                                                                                                              ccATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAg                                       <  1:5511174/150‑1 (MQ=255)
                                                                                                              ccATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAg                                       <  1:5510750/150‑1 (MQ=255)
                                                                                                              ccATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAg                                       <  2:1856172/150‑1 (MQ=255)
                                                                                                                aTTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGcc                                     >  1:5047753/1‑150 (MQ=255)
                                                                                                                    aCCAGGGAAACCTGCGCACCACGCGTGGTGCACCGGGTTACCGGGGTTGGATAGAGGAGTGCGTCACGACCACAACTGGCACAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAAt                                 >  2:1554917/1‑150 (MQ=255)
                                                                                                                    aCCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTcgcg                                                                <  1:7642568/119‑1 (MQ=255)
                                                                                                                    aCCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAAt                                 >  1:2661992/1‑150 (MQ=255)
                                                                                                                    aCCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAAt                                 <  1:4004557/150‑1 (MQ=255)
                                                                                                                    aCCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAAt                                 <  1:2278921/150‑1 (MQ=255)
                                                                                                                     ccAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGACACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAAtt                                >  1:8839658/1‑150 (MQ=255)
                                                                                                                     ccAGGGAAACCTCCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGCCCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAAtt                                >  1:3835809/1‑150 (MQ=255)
                                                                                                                      cAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGGTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAAttt                               >  2:5181361/1‑150 (MQ=255)
                                                                                                                       aGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTg                              <  1:6205734/150‑1 (MQ=255)
                                                                                                                        gggAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGc                             >  1:9115723/1‑150 (MQ=255)
                                                                                                                         ggAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGct                            <  2:8001648/150‑1 (MQ=255)
                                                                                                                         ggAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGct                            <  1:5713998/150‑1 (MQ=255)
                                                                                                                              ccTGTGCACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGtt                                                      <  1:1544557/119‑1 (MQ=255)
                                                                                                                              ccTGTGCACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGa                       >  1:6643130/1‑150 (MQ=255)
                                                                                                                              ccTGTGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGa                       >  2:1544377/1‑150 (MQ=255)
                                                                                                                              ccTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGa                       <  1:8892078/150‑1 (MQ=255)
                                                                                                                              ccTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGa                       <  1:6129369/150‑1 (MQ=255)
                                                                                                                               cTGTGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGAt                      >  2:913800/1‑150 (MQ=255)
                                                                                                                               cTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGAt                      <  1:7483378/150‑1 (MQ=255)
                                                                                                                                tgtgcACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATc                     >  1:6025508/4‑150 (MQ=255)
                                                                                                                                tgtgcACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATc                     >  1:7792873/4‑150 (MQ=255)
                                                                                                                                tGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATc                     >  2:8770197/1‑150 (MQ=255)
                                                                                                                                tGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATc                     <  1:3815038/150‑1 (MQ=255)
                                                                                                                                tGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATc                     >  1:105881/1‑150 (MQ=255)
                                                                                                                                tGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATc                     <  2:388663/150‑1 (MQ=255)
                                                                                                                                 gtgcACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCg                    >  1:6883046/3‑150 (MQ=255)
                                                                                                                                   gcACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGc                  <  2:442154/150‑1 (MQ=255)
                                                                                                                                     accacGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTa                                                                              <  1:2336312/88‑1 (MQ=255)
                                                                                                                                     accacGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTCGATCGGccc                >  2:7159355/1‑150 (MQ=255)
                                                                                                                                       cacGCGTGGTGTACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCgt              >  2:7717259/1‑150 (MQ=255)
                                                                                                                                       cacGCGTGGTGGACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTCGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCgt              >  2:7717453/1‑150 (MQ=255)
                                                                                                                                       cacGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTCGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCgt              >  2:8986475/1‑150 (MQ=255)
                                                                                                                                       cacGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCgt              >  2:6781895/1‑150 (MQ=255)
                                                                                                                                           cgTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTAc          <  2:7337442/150‑1 (MQ=255)
                                                                                                                                             tggtgCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACca        <  2:324723/150‑1 (MQ=255)
                                                                                                                                                  cacaGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCgt              <  2:8912616/139‑1 (MQ=255)
                                                                                                                                                  cacaGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGc   >  1:4908880/1‑150 (MQ=255)
                                                                                                                                                  cacaGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGc   >  1:8583525/1‑150 (MQ=255)
                                                                                                                                                  cacaGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGc   >  2:4793130/1‑150 (MQ=255)
                                                                                                                                                  cacaGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGATTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGc   >  2:7599693/1‑150 (MQ=255)
                                                                                                                                                   acaGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCa  >  2:1638777/1‑150 (MQ=255)
                                                                                                                                                   acaGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCa  >  2:3313894/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCA  >  NC_002947/3455298‑3455594

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: