Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,719,840 G→A 14.6% intergenic (+73/‑208) yiaG → / → cspA HTH_CROC1 family putative transcriptional regulator/RNA chaperone and antiterminator, cold‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,719,8400GA14.6% 184.0 / 2.7 82intergenic (+73/‑208)yiaG/cspAHTH_CROC1 family putative transcriptional regulator/RNA chaperone and antiterminator, cold‑inducible
Reads supporting (aligned to +/- strand):  ref base G (39/29);  new base A (2/11);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 6.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCGAAAGGCACACTTAATTATTAAAGGTAATACACTATGTCCGGTAAAATGACTGGT  >  NC_000913/3719613‑3720068
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                     
gATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATAAAAAAAATGa                                                                                                                                                                                                                                      >  2:303514/1‑229 (MQ=255)
    aaaaaTCAGTGTTTTCGCCCGGGTGTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGaaa                                                                                                                                                                                                                                    <  2:346814/227‑1 (MQ=255)
     aaaaTCGATGATTTCGCCCGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAaat                                                                                                                                                                                                                                 >  1:163554/1‑229 (MQ=255)
          cGATGATTTCGCCCGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACGGTTTTCCAACCGCTTCATCATAAATATGACAAATAAtt                                                                                                                                                                                                                             >  1:332986/1‑228 (MQ=255)
          cGATGATTTCGCCCGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATCAATATGAAAAATAAtt                                                                                                                                                                                                                             >  1:591131/1‑228 (MQ=255)
          cGATGATTTCGCCCGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACAGTTTTCTAACCGATTAATCTTAAATATGAAAAAAAAtt                                                                                                                                                                                                                             >  2:108154/1‑228 (MQ=255)
          cGATGATTTCGCCCGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGAAAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAAAAAtt                                                                                                                                                                                                                             >  2:389978/1‑228 (MQ=255)
          cGATGATTTCGCCCGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCAATAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGATGTTTACGGTCCTGATGACAGGACCTTTTTCCAACCGATTAACCTTAAATAGAAAAAAAAAtt                                                                                                                                                                                                                             >  2:597408/1‑228 (MQ=255)
          cGATGATTTCGCCCGGGTTTTGGGCGAATCAGTCGCCATGGTAAAGGAATGGGAATCAAGACGCGTGAAGCCTTCAAGTGCCGCACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTACGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCAAACCGCTTAATCATAAATATGAAAAAAAAtt                                                                                                                                                                                                                             >  2:76965/1‑228 (MQ=255)
          cGATGATTTCGCACGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGACGTTTACGACCTGGCTGACAGGACCGTTTTCCAACCGTTTATTCAAAACTAGGAACaataaat                                                                                                                                                                                                                             >  2:695047/1‑226 (MQ=255)
            atgatTTCGCCCGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGCATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAACTATGAAAAATAAttgt                                                                                                                                                                                                                           >  2:554597/1‑228 (MQ=255)
             tgatTCCGCCGGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAAGTGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAttgttg                                                                                                                                                                                                                         <  2:373185/229‑1 (MQ=255)
                     cccGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCAAGACGCGTGAAGCCTTCCAGTGCCGAACTAAAATTGCTGCGTTTGATTCAAGCCAACCCGCCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCCATTTATTACTGTTCACGGTCCTGATGAAAGGACCGTTTTCCAGCCGATTAATAATAACTATACAaattaatttcttcttaccca                                                                                                                                                                                                                 >  2:184358/1‑217 (MQ=255)
                                    tATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCtt                                                                                                                                                                                                 >  1:593163/1‑229 (MQ=255)
                                    tATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGGTTTCCACCCGATTAATCATAAATTTGAAAAAAAATTGTTGAATTACCCCCCAATGCGTGGGtt                                                                                                                                                                                                 >  2:689967/1‑229 (MQ=255)
                                        agccGCCATGGTAAAGGAATTGGAATCCAGTCGCTGGAAGCCTTAAAGCGCCGAACAAAACTTGATGCGTTTGATCTAAGCCAACCCGGCATTAAGTAAGCAGTGGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGTTCCTGATGACAGGACGTTTTTCAAACCGATTATTCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATg                                                                                                                                                                                             <  1:468227/226‑1 (MQ=255)
                                          tgaataagagaaaggaaaggGGAATCCAGCGGCGTGAAGCCTTAAAGTGCGAAACTAAAATTGAGGCGTTTGTTTCAAGCCAACCCGGCATTAAGTAAGCATTTGATGGAACAGACTTTTTTCCTCTTTATTGTTGTTTACGGTCCTGATGACAGGCCCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATg                                                                                                                                                                                             <  1:259565/218‑1 (MQ=255)
                                                      agtaaTGGGAATCCAGCCGGGTTAACCTTTCAGTTGCAAAATAAAAGTTGATGCGTTTAATTTAAGCAACCCGGGCATTAAGTAAGGAGTTGATGGAATAGCTTTTTATCCACTTTATTGATGTTAACGGTCCTGATGACGGGCCCGTTTTCCAACCACTTAATCATAAATAGGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGttt                                                                                                                                                                               <  2:324434/226‑1 (MQ=255)
                                                           tGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCAAATAAACGGTTTGACGt                                                                                                                                                                          >  1:415449/1‑229 (MQ=255)
                                                             gggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATCCGTGGCTTAATGCACATCAACGGTTTGACGTa                                                                                                                                                                         >  1:1571/1‑229 (MQ=255)
                                                              gaaaTCCAGACGCGTGAACCCTTCAAGGGCCGAACTAAAATTGATCCGTTTAATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCATTTTATTGCTTTTTACGGTCCTGATGACAGGACCGTTTTCCACCCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTa                                                                                                                                                                         <  1:43967/226‑1 (MQ=255)
                                                                aaaccagggggggaaacccTTAAAGGTCCAAACTAAAAGTGATCGTTTTGATCACAGCACAACCGGGATTAAGTAAGAAGTTGATGGAATAGACTTTTACCCACTTTATCGCATTTAAGGGTCCTGATGAAAGGACCGTTTTCAAACCGATTAATCATAAATATGAAAAAAAATTGTTGCAGCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAg                                                                                                                                                                      <  2:622731/212‑1 (MQ=255)
                                                                aaTCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAg                                                                                                                                                                      >  1:657686/1‑229 (MQ=255)
                                                                aaTCCAGACGCGTAAATCCTTCAGTTCCCGAACTAAAATTGATTCGTTGGATTCAAGCCACCCGTGCATTAAGTAAGCAGTTGATGGAATAGCCTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATAGGAAAAATAATTGTTGCATCCCCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAg                                                                                                                                                                      <  1:108154/229‑1 (MQ=255)
                                                                   ccAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTATTCTTAAATATGAAAATAAATTGTTGCATAACCCCCCAATCCTTGGCTTAATGCACATaaaagccgttatctaaaaca                                                                                                                                                                    >  2:259565/1‑211 (MQ=255)
                                                                       cggcgGGAAGCCTTCAATTGCAGAACTAAAATTGATGCGTTTGATTCAACCCAACCCGGCATTAAGTAAGTAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGAAAGGACCTTTTTCCACCCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTa                                                                                                                                                               <  2:562954/227‑1 (MQ=255)
                                                                        cgcgTAAACCCTTCAAGTGCCGAACTAAATTTGATGCGTTTAATTCAAGCCAACCCGTCATAAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCACCCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTa                                                                                                                                                               <  1:101855/228‑1 (MQ=255)
                                                                           gtgaatcctgcaatgtctgAACTAAAATTTATGGTTTTGATTAAAGCCAACCTGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAg                                                                                                                                                            <  1:307230/210‑1 (MQ=255)
                                                                           gagAAGCCTTTAGGTGTCGAATCAAAATTGATGCGTTTAATGCAAGACAACCCGGATTAAAGTAAGAAGTGGATGGAATAGGCTTTAATCCACTTTAGTTATTTTGACGGTCCTGAACGCAGTACCGTTTGCCAACAGATCAAACAAAAAAATGAAAAATAATTGTTTCTCCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAg                                                                                                                                                            <  1:695048/226‑1 (MQ=255)
                                                                           aggAAGGCTTAAAGTTACGGAGTAAAGTTGAGGGTTTTAATTCAAGCCACCCCGGGTTTAAGTAAGCCGTTGATGGAATGGATTTTTACCCGCTTAATTGCTTTTTACGTTCATGATGACAGGACGGTTTTCAAAACGATTATGCATAAATATGAAAAATAATTTTTGCAGCAGCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAg                                                                                                                                                            <  1:630193/226‑1 (MQ=255)
                                                                            tGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATACTTGTTGCATCACCCGCCAATGCGTGGCTTAATGAACATAAACGGTTTAACGTACAACCCCTTAAAGCa                                                                                                                                                          >  2:755623/1‑229 (MQ=255)
                                                                            tGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGAATTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGAACCGTTTTCCAACCGCTTAATCTTAAATATGAAAAATAATTGTTGCATTACCCGCCAATACGTGGACTAATGCTAATCAACCGGTTATCGTACAGACAAttcaccaa                                                                                                                                                          >  1:352850/1‑223 (MQ=255)
                                                                            gaaaccctctaatgtccaaattAAATTTGGTGCTTTTATTTCAACCCACCCCGGCTTTAAGTAAGCGTTGAAGGGAATGGGTTTTTATCCCTTTTATGTCTGTTTACGGTACTAAGTACAGGACCGTTTTCCAACCGATTATTCAAAAAAATAAAAAAAATTTGTGCCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCa                                                                                                                                                          <  2:285424/208‑1 (MQ=255)
                                                                                   ttCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCATTgtag                                                                                                                                                    >  1:673924/1‑228 (MQ=255)
                                                                                          gCCAAACTAAAATTGCTGTTGTTACTTCAAGCCAACCCGGCATTAAGTAAGCCGTTGAGGGAATAGATTTTTTGCCATTTTATTTCTGTTTACGGTCCGGATTACAGGTCCGTTTTCCACCCGATTAATCATAAATATGAAAAAAAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCaa                                                                                                                                            <  1:433019/229‑1 (MQ=255)
                                                                                            cgtaCTAAAATTGTTGCTTTTATTTCAACGCAACCCTGACTTAAGTAAGGAGTTAAGTGATATGCTTTTTTCCCATGTTATTGCTGTTTAGGTTGCTAATGACAGGCCCTTTTCCCACCCAATTAATAATAAATTTAAAAAAAAATTGTGTCTTCACCTGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGt                                                                                                                                          <  2:332298/226‑1 (MQ=255)
                                                                                            caaatgaaaaTGGATGCGTTTGATTCAACCCAACCCGGATTTAATTAAGCAGTTGATGGATTAGATGTTTATCCACTTTATTGCGTTTTACGGTCCTGATAAAAGGACGGTTTTCCAACCGATTAATCATAAATATAAAAAATAATTGTTGCAACACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGt                                                                                                                                          <  1:115775/223‑1 (MQ=255)
                                                                                                         aTGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATAACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCTTTAAAGCCGTGAAGTAAGGCAAAACCCTTCAAAAGta                                                                                                                             >  2:460653/1‑228 (MQ=255)
                                                                                                               tttttttacccccacctgtctttAAGTAAGTGGTTGATGGAATAGATTTTTACCCCTTTTATTTCTTTTTACGGTCCTGAGGAGGGGCCCGTTTCCCAACGTATCAATAATAAATATGAAAAAGAATTGTTGCAGCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGtt                                                                                                                       <  2:332986/208‑1 (MQ=255)
                                                                                                                  tttCAAGCCAACCCGGCATTAAGTAAGCTGTGGATGGAATGACTTTTTACCCCTTTTATTGCTTTTTACGGTCATGATGACAGGCCCGTTTTCCAACCGATTAATCAAAAAATGGAAAAATAAGTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGAt                                                                                                                    <  2:1571/228‑1 (MQ=255)
                                                                                                                  atcaaaGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCATTTTATTGCTGTTTACGGTCCTGATGACAGGCCCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGAt                                                                                                                    <  1:460653/225‑1 (MQ=255)
                                                                                                                  aTTCAAGCCAACCCGGCTTTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTTTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATAGGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGAt                                                                                                                    <  1:653725/229‑1 (MQ=255)
                                                                                                                    taaatccAACCCGGCTGTAAGGAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCACCCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATAc                                                                                                                  <  1:597409/224‑1 (MQ=255)
                                                                                                                     taatcacaCCCGGCATTAAGTACGCAGTTGAGGGAATAGACTTTTATCCTATTTATTGCTGTTTACGTGCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATAcc                                                                                                                 <  1:303514/222‑1 (MQ=255)
                                                                                                                        gCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAATTCCCTTCAAAAGTTATCGTTAAAACCCCt                                                                                                              >  2:626145/1‑229 (MQ=255)
                                                                                                                        gCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGACCATCAACGGTTTGACGttaag                                                                                                                                                                      >  2:501082/1‑169 (MQ=255)
                                                                                                                               cGGAATTAAGTATGAGGTTGGTGGAAAAGGTTTTTATCCACTTTTATTCTTTTTACGGTCGTTATTACAGGACCGTTTTCCAATCGATTAATCATAAAAATGAAAAATAATTGTTGCATCCCCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAATAGTTATCGTTGATACCCCTCGTAGTg                                                                                                       <  1:689968/229‑1 (MQ=255)
                                                                                                                                gtgctttaagtaaggcttgaaTGAAATAGACTTTTACCACCTTTATTGCGTGTTAGTTCCCTGATGACGGGACCTTTGTCCACCCGATTAATCATAAATATAAAAAATATTTTTTGCAACCCCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCCGTGTAGTAATGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTg                                                                                                       <  2:176183/209‑1 (MQ=255)
                                                                                                                                 gCATTAAGTAACCAGTTGATGGAATGGACTGTTATCCATTTTATTTCCTTTTACGGTCCAATAGACAGGACCGTTTTCCAACCGAGTAATCATAAATATGAAAAATAATTGTTGCACCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGca                                                                                                     <  1:755624/229‑1 (MQ=255)
                                                                                                                                  cATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGca                                                                                                     >  1:161538/1‑228 (MQ=255)
                                                                                                                                  cATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGca                                                                                                     >  1:5294/1‑228 (MQ=255)
                                                                                                                                    atataataaagatgtggtGGAATAGACTTTTTTCCATTTTTTTGTTTTTAACGGGCCGTAATAAGGGACCTTTTCCCCACCGATTAATCAAAAATAAGAAAAATATTTGTTCCAACACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGcaca                                                                                                   <  2:161538/212‑1 (MQ=255)
                                                                                                                                            gCAGTTGAGGGAATAGACTGTTAACCATTTAATTGATGTTTACGGCCCTGATGACAGGACCTGTTTCCAACCGATTAACAAAAAATTTGAAAAATAATTGTTGCATCATCCGCCAATGCGTGGGTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTa                                                                                           <  1:95080/228‑1 (MQ=255)
                                                                                                                                                ttGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGACCATTCCTTTaa                                                                                          >  1:595456/1‑225 (MQ=255)
                                                                                                                                                           aGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCt                                                                             >  1:176534/1‑228 (MQ=255)
                                                                                                                                                                     tCCTCTTTATTGCTTTTTAGGGTCATGTAGACCGGACCGTTTGCCAGCCGATTAATCATAAATAAGAAAAATAATTGTGTCATCACCCGCCAATGCTTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACg                                                                   <  1:76965/229‑1 (MQ=255)
                                                                                                                                                                     accctcTTTATTGCTTTTTAGGGTCCGAATGACAGGACCGTTTTCCAACCAATTAATAAAAAAATATAAAAATAATTGTTGCATCCCCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCa                                                                     <  1:284222/223‑1 (MQ=255)
                                                                                                                                                                     accccttttATTGTTTTTTACTGTCCGGTGGCAAGGGCCTTTTCCCAACCGATTAAACAAAAATAGAAAAAATAATTGTTGAATCCCCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCa                                                                     <  1:20460/222‑1 (MQ=255)
                                                                                                                                                                                 tgttttttCGGTTCTATTGACGGGCCCGTTTTCCACCCAATAAATAAAAAATGTAAAAAAAAATTTTGGCCTCACCCGCCAATGCGTGGCTTAATGAACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCg                                                        <  2:657685/225‑1 (MQ=255)
                                                                                                                                                                                  ttGTTTAGGGTCCTAATACAAGACCCTGTTCCCACCTGATTAGTAAAAAAAATGAAAAATAATTGTACACTCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACAATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCGa                                                       <  2:515519/228‑1 (MQ=255)
                                                                                                                                                                                     tttACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGAATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTACTACCCCTCGTGGGCCACATTCCGTTAACGCTTCAAAATCCGTAAAGCACTCCATATCGCCGAAAg                                                    >  1:309333/1‑229 (MQ=255)
                                                                                                                                                                                      aataacgccTGATGACTGGACCTTTTTCCAACCGATTAATAAAAAATATAAAAAATAATTTGTGCACCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCGAAAgg                                                   <  2:342946/222‑1 (MQ=255)
                                                                                                                                                                                       ttcGGTTCTGGTAACGGGCCCTTTTCCCAACCAATTAAAAAAAAATTGAAAAAATAATTGTGGCATCCCCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCGAAAGGc                                                  <  2:673923/227‑1 (MQ=255)
                                                                                                                                                                                            tCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAATAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTATTGCACATTCCTTTAACGCTTCAAAATCTGAAAAGCACGCCATATCGCCGAACGGCACACt                                             >  1:103804/1‑229 (MQ=255)
                                                                                                                                                                                               tgatgaCAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGAAACCCATCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGAACGCAATATCGCCGAAAGGCACACTTaa                                          >  1:191378/1‑229 (MQ=255)
                                                                                                                                                                                               tgaataCAGGACCGTTTTCCAACCGATTAATCATAAATATAAAAAATAATTGTTGCATCCCCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCGAAAGGCACACTTaa                                          <  1:389978/224‑1 (MQ=255)
                                                                                                                                                                                                 atgaCAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCg                                                                                                            >  1:148639/1‑161 (MQ=255)
                                                                                                                                                                                                 atgaCAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAACCACCCCAAATCCCCGAAAGGCCAACTTAatt                                        >  2:103716/1‑229 (MQ=255)
                                                                                                                                                                                                    cgcggcccttttcCCAACTTATTAATAAAAAAAATAAAAAATAAGTGTGGCTCCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCGAAAGGCACACTTAattatt                                     <  2:591130/220‑1 (MQ=255)
                                                                                                                                                                                                    aaaGGACTGTTTTCCCACCAATTAATCAAAAATAAAAAAAATAATTGTTGCATCCCCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCGAAAGGCACACTTAattatt                                     <  2:176534/227‑1 (MQ=255)
                                                                                                                                                                                                      cacaCCGTTTTCCAACCAAATAAAAAAAATTAGAAAAAATATTGTTTTACCCCACCGCCAAAGCGTGGCTTAATGCACATCAACGGTTTGACGTCCAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCGAAAGGCACACTTAATTATTaa                                   <  2:635807/226‑1 (MQ=255)
                                                                                                                                                                                                      aGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAACCTGTAAAACACGCCATATCCCCGAAAGGCACACTTAATAATTaa                                   >  2:432568/1‑229 (MQ=255)
                                                                                                                                                                                                       ggccattttacaaaCCGTTTAATCAAAAAGTAGAAAAATAAGTGTTGCATCCCCCGCCAATGCGTGGCTTAATGCAAATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCGAAAGGCACACTTAATTATTaa                                   <  1:626146/216‑1 (MQ=255)
                                                                                                                                                                                                        gACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCGAAAGGCCCCCTTAATTATTAAAg                                 >  1:733079/1‑229 (MQ=255)
                                                                                                                                                                                                        gACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCGAAAGGCACACTTAATTATTAAAg                                 >  1:94929/1‑229 (MQ=255)
                                                                                                                                                                                                        gACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGAAGAGCACATTCCTTTAACCCTTCAAAATCTTTAAAGTCCGCCATTTAGCCGACAGGCACCATTAattctgtct                                  >  1:342946/1‑222 (MQ=255)
                                                                                                                                                                                                        gACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCACTCTTATGGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCAATATCGCCGAAAGGCACACTTAATTATTAAAg                                 >  1:286037/1‑229 (MQ=255)
                                                                                                                                                                                                        gACCGTTTTCCAACCGACTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGTTACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCGCAAGGCACACTTAATTATTaagg                                 >  1:237272/1‑227 (MQ=255)
                                                                                                                                                                                                                    aCCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTAAAGAGTTATCGTTGATACCCCTCGTAGTGAACATTCCTTTAACGCTTCAACATCTGTAAAGCACGCCATTTCGCCGAAAGGCACACTTAATTCTTaa                                   >  1:635808/1‑215 (MQ=255)
                                                                                                                                                                                                                         aTTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAAAGCAAATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCGCAAGGCACCCTTAATTATTAAAGGAAATACACTTTGta                  >  1:329324/1‑227 (MQ=255)
                                                                                                                                                                                                                            taTCATAATTATGAAAATGAATTGTTGCCTCCCCCGCAAATGGGTGGCTTAATGCACAACAACGGTTTGACTTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCGATAGGCACACTTAATTATTAAAGGTAATACACTATGTCCGGt              <  2:286037/228‑1 (MQ=255)
                                                                                                                                                                                                                                        tGAATAATAATAGTTGCATCACCCGCCAATGCGAGGCTTAATGCACAGCAACGGTTTGACGGACAGACCATTATAGCAGTGTAGTAAAGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCTACGTGCACATTCCTTTAACACTTCAAAATCTGTCAAGCCCGCCACATCCCCGAAAGGCAAAATTACTTATTAAAGGAAATACACTATTTCCGGTaaactaacagg   >  2:542673/1‑221 (MQ=255)
                                                                                                                                                                                                                                         gAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCGAAAGGCACACTTAATTATTAAAGGTAATAAACAATGTCCGGAAAAAGAACTgga  >  2:467355/1‑228 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                     
GATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGATACCCCTCGTAGTGCACATTCCTTTAACGCTTCAAAATCTGTAAAGCACGCCATATCGCCGAAAGGCACACTTAATTATTAAAGGTAATACACTATGTCCGGTAAAATGACTGGT  >  NC_000913/3719613‑3720068

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