Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,081,473 T→C 9.7% intergenic (‑311/+73) sdaC ← / ← amaC serine:H+ symport permease/D‑lysine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,081,4730TC 9.7% 248.2 / 3.2 92intergenic (‑311/+73)sdaC/amaCserine:H+ symport permease/D‑lysine aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (52/27);  new base C (2/7);  total (55/37)
Fisher's exact test for biased strand distribution p-value = 2.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATGCACACAC  >  NC_002947/4081324‑4081619
                                                                                                                                                     |                                                                                                                                                  
ggCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCGTTAGGGCCGCAAAGCGGCCCCAGCAAATTCAGCCGCCAAACGGCAATAGTAGGGGCGCTACGCCCCCgttt                                                                                                                                                    >  1:4951805/1‑150 (MQ=37)
  cATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTTCCTAGGCAAATGGCCGCCAAACGGCCCCAGCAAATTAAGAGGCAAAGCTGCAAAAGTAGGGGCCCTCCGCACCCctttag                                                                                                                                                  >  2:1996635/1‑148 (MQ=255)
  agtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTcg                                                                                                                                                  >  1:4851968‑M1/55‑150 (MQ=255)
   gtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTcgc                                                                                                                                                 >  2:480695‑M1/54‑150 (MQ=255)
   aTTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCACCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCACCGCACCCATTTcgc                                                                                                                                                 >  2:4407793/1‑150 (MQ=255)
      gCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCATGCGTCGCGATAGGGCCGCAAACCGGCCCCACCAAATTAATCGGCAAAGCTGCAAAAGTAGGGGCGCACCGCACCCTTTTAGCGAc                                                                                                                                              >  2:2172539/1‑150 (MQ=37)
      actcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGAc                                                                                                                                              >  1:2424032‑M1/51‑150 (MQ=255)
          gccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCaa                                                                                                                                          >  1:5187495‑M1/47‑150 (MQ=255)
          gccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCaa                                                                                                                                          >  2:1225111‑M1/47‑150 (MQ=255)
            gtCGATCGGTGGAATGGAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGTTAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAgg                                                                                                                                        <  2:1398532/150‑1 (MQ=255)
            caacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAgg                                                                                                                                        >  1:2095323‑M1/45‑150 (MQ=255)
                        aaTGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCCAGAGGGCCCCAAACCCGCCCCACCCAATTCAGCGGCAAAGCGGCAAACGGAGGGGCGCGCCCCAACCGGTTCGCGCCGCCACGCCGCGCCCACaa                                                                                                                            >  2:982200/1‑150 (MQ=255)
                           gTAGCGGCCTGGGCGCTGGTGGAGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCCGCGGCAAAGCTGCAATCGTTGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGt                                                                                                                         <  2:4850937/150‑1 (MQ=25)
                               cggcgggggCGGCTGGGGGGGGGCGGGGGGGCGCGGCGTGGCGTCGCGATCGGGCCGCAACGCGGCCCCCGCATTTTCCGCGGCAAAGCTGCTTTGGTTGGGGCGCTCGGCCCCCTTTTCGCGACGCACGGCCGCTCCCACAAAGTCCaa                                                                                                                     <  2:4642097/144‑1 (MQ=255)
                                           cTGTGGGGCGGCCGGTGGGCGCGGCCTTGCGTCGCGATCGGGCCGCAACGCGGCCCCAGCAATTTCCGCGGAAACGCTGCAATCGTCGGGGCGCTCCGCACCCTTTCCGCGCCGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCcaca                                                                                                         <  2:407205/150‑1 (MQ=255)
                                                 ggcggcCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGAAACGCGGCCCCAGCAATTTCCGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAggg                                                                                                   <  2:4394320/150‑1 (MQ=255)
                                                    ggcCTGGGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAACGGGGCCCCAGCATTTTCCGCGGCAAAGCGGCAATCGTAGGGGCGCGCCGCCCCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGcc                                                                                                <  2:793216/150‑1 (MQ=37)
                                                         GTGGGAGCGGCCTTGCGCCGCGATAGGGCCGCAACGCGGCCCCAGCAATTTCAGCGGAAAAGCTGAAATCGTAGGGGCGCTCCGCACCATTTCCGCGCCGCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGt                                                                                           >  2:4165976‑M1/1‑150 (MQ=255)
                                                          tgGGCGCGGCCTTGCGGCGCGAACGGGCCGCAACGCGGCCCCAGCAATTTCCGCGGCAAAGCGGAAAAGGTCGGGGCGCGCCGCACCCTTTGCGCCCCGCAAGCCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGtt                                                                                          <  2:54129/150‑1 (MQ=255)
                                                           GGGCGCGGCCTTGCGTCGCGATAGGGCCGCAACGCGGCCCCCGCATTGTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTGCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTg                                                                                         >  2:1584353‑M1/1‑150 (MQ=255)
                                                               GCGGCCTTGCGTCGCGATAGGGCCGCAACGCGGCCCCAGCAATTTCAGCGGCAATGCTGCAATCGTAGGGGCGCACCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGccc                                                                                     >  2:3466401‑M1/1‑150 (MQ=255)
                                                                 GGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGt                                                                                   <  1:3109105‑M1/150‑1 (MQ=255)
                                                                   CCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTg                                                                                 <  2:3628246‑M1/150‑1 (MQ=255)
                                                                   CCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAGGTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCCCAGGAGGGGCCGCAGTTGGCCCGTTg                                                                                 <  2:3628267‑M1/150‑1 (MQ=255)
                                                                    CTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTgg                                                                                <  1:38877‑M1/150‑1 (MQ=255)
                                                                     TTGCGTCGTGTTTGGGCCGAAAAGCGGCCCCAGCATTTTCCGCGGCAATGGTGAATTCTTTGGGGCGCTCCGCACCCTTTTCGCGACGAATGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGc                                                                               >  2:983017‑M1/1‑150 (MQ=255)
                                                                     TTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTCCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGc                                                                               >  2:4064833‑M1/1‑150 (MQ=255)
                                                                         GTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGc                                                                                                 >  2:1334156‑M1/1‑128 (MQ=255)
                                                                         GTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGc                                                                                                 <  1:1334156‑M1/128‑1 (MQ=255)
                                                                         GTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGt                                                                           <  2:3327401‑M1/150‑1 (MQ=255)
                                                                         GTCGCGATAGGGCCGCAAAGCGGACCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGt                                                                           <  1:3769698‑M1/150‑1 (MQ=255)
                                                                          TCGCGATAGGGCCGAAAAGCGGCCCCAGAAATTTCAGCGGAAAACCTGCAATCTTAGGGGCGCTCCGAACCCTTTTCGCGACGAAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGtt                                                                          >  2:1674390‑M1/1‑150 (MQ=255)
                                                                             CGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACt                                                                       <  1:5381264‑M1/150‑1 (MQ=255)
                                                                                tggggCCGCTTTGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCt                                                                    >  2:4915916/3‑150 (MQ=255)
                                                                                   ggCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGaa                                                                 >  1:1930194/1‑150 (MQ=255)
                                                                                   ggCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGaa                                                                 >  1:664758/1‑150 (MQ=255)
                                                                                   ggCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGaa                                                                 >  1:1930139/1‑150 (MQ=255)
                                                                                   ggCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGaa                                                                 >  1:1930124/1‑150 (MQ=255)
                                                                                    gCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAAc                                                                >  2:1657004/1‑150 (MQ=255)
                                                                                    gCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAAc                                                                >  2:3294473/1‑150 (MQ=255)
                                                                                    gCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAAc                                                                >  1:1585635/1‑150 (MQ=255)
                                                                                     ccGCAAAGCGGCCCCAGCAATTTCCGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACg                                                               <  2:1585635/150‑1 (MQ=255)
                                                                                     ccGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCATCCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACg                                                               >  2:735422/1‑150 (MQ=255)
                                                                                     ccGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACg                                                               >  1:1888495/1‑150 (MQ=255)
                                                                                        cAACGCGGCCCCCGCAATTTCCGCGGCAACGCTGCATTGGTCGGGGGGCGCCGCACCCTTTGCGCGACGCATGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCa                                                            <  2:1054022/150‑1 (MQ=255)
                                                                                            gCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGa                                                        >  1:2809314/1‑150 (MQ=255)
                                                                                             cGGCCCCTGCTATTTCCGCGGCAAAGCTGCAATCGTAGGGGCCCCCCGCACCCTTTCCGCGACGCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGaa                                                       <  1:2712324/150‑1 (MQ=255)
                                                                                             cGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGaa                                                       >  1:596633/1‑150 (MQ=255)
                                                                                               gCCCCAGCAATTTCCGCGGCACAGCTGCAAGCGTCGGGGCGCCCCCCACCCTTTCCGCGACGCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACg                                                     <  1:4971197/150‑1 (MQ=255)
                                                                                                ccccAGCAATTGCAGCGGCAACGCTGCAATCGCCGGGGCGCTCCCCCCCCTTTCCGCGACGCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGc                                                    <  2:1357383/150‑1 (MQ=255)
                                                                                                  ccAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCtt                                                  >  1:4254093/1‑150 (MQ=255)
                                                                                                  ccAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCtt                                                  >  2:4272535/1‑150 (MQ=255)
                                                                                                  ccAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCtt                                                  >  1:4688129/1‑150 (MQ=255)
                                                                                                  ccAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCtt                                                  >  1:3547683/1‑150 (MQ=255)
                                                                                                  ccAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCtt                                                  >  1:4254063/1‑150 (MQ=255)
                                                                                                  ccAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCtt                                                  >  2:5198991/1‑150 (MQ=255)
                                                                                                  ccAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCtt                                                  >  2:2659805/1‑150 (MQ=255)
                                                                                                   cAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTc                                                 >  2:668632/1‑150 (MQ=255)
                                                                                                   cAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTc                                                 >  2:2698545/1‑150 (MQ=255)
                                                                                                   cAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTc                                                 >  1:1103278/1‑150 (MQ=255)
                                                                                                    aGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCg                                                >  2:85898/1‑150 (MQ=255)
                                                                                                    aGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCg                                                <  2:5381264/150‑1 (MQ=255)
                                                                                                    aGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCg                                                <  1:5198991/150‑1 (MQ=255)
                                                                                                    aGCAATTCCCGCGCCACTGCTGCATTCGACGGGGCGCACCGCACCATTTCCGCCACCCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCg                                                <  1:3569791/150‑1 (MQ=255)
                                                                                                     gCATTTTCCGCGGCAATGGTGCATTGGTTGGGGCGCTCCGCACCCTTTTGGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCgg                                               <  2:5250244/150‑1 (MQ=255)
                                                                                                      cactttCCGCGGCACCGCTGCAATCGTCGGGGCGCGCCGCACCCTTTCCGCGCCGCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGc                                              <  2:2326407/147‑1 (MQ=255)
                                                                                                         tttCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAAc                                           >  2:2395038/1‑150 (MQ=255)
                                                                                                          ttCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAAcc                                          >  1:1273028/1‑150 (MQ=255)
                                                                                                           tCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAAccc                                         >  2:1679137/1‑150 (MQ=255)
                                                                                                           acAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTACGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAAccc                                         <  2:3622505/149‑1 (MQ=255)
                                                                                                              gCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGtt                                      >  1:2363904/1‑150 (MQ=255)
                                                                                                                 gCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCa                                   >  2:4618456/1‑150 (MQ=255)
                                                                                                                 gCAAAGCTGCAACCGTAGGGGCGCTCCCCACCCTTGTCGCGCCGCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCa                                   <  2:3698702/150‑1 (MQ=255)
                                                                                                                  cAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAg                                  <  2:3532844/150‑1 (MQ=255)
                                                                                                                  cAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAg                                  <  2:4360446/150‑1 (MQ=255)
                                                                                                                  cAAAGCTGAAACCGTCGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTAACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAg                                  <  2:791470/150‑1 (MQ=255)
                                                                                                                    aaGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACACAGTCCAACCGCGCCCCACCGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGTACGCTTCGGCAACCCGTTGCAGGt                                <  1:3628267/150‑1 (MQ=255)
                                                                                                                    aaGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGt                                <  1:3628246/150‑1 (MQ=255)
                                                                                                                       cTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTgg                             >  2:5354142/1‑150 (MQ=255)
                                                                                                                        tGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGc                            >  2:1826009/1‑150 (MQ=255)
                                                                                                                           aaTCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCGGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCgg                         >  2:4562789/1‑150 (MQ=255)
                                                                                                                           aaTCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCgg                         >  1:766985/1‑150 (MQ=255)
                                                                                                                             tCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGgc                       >  1:3106195/1‑150 (MQ=255)
                                                                                                                             tCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGgc                       >  1:3106004/1‑150 (MQ=255)
                                                                                                                              cGTCGGGGCGCTCCGCACCCTTTCCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGgcg                      <  2:319110/150‑1 (MQ=255)
                                                                                                                              cGTCGGGGCGCTCCGCACCCTTGTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGgcg                      <  2:4693690/150‑1 (MQ=255)
                                                                                                                                tAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCa                    >  2:4729412/1‑150 (MQ=255)
                                                                                                                                tAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCa                    >  2:3078257/1‑150 (MQ=255)
                                                                                                                                   gggCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGc                 >  1:2226724/1‑150 (MQ=255)
                                                                                                                                     gcgcTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCg               <  2:4189229/150‑1 (MQ=255)
                                                                                                                                          ccGCACCCTTTTCGCGACGCAAGGCCGCTCCCACACAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACAt          <  2:2226724/150‑1 (MQ=255)
                                                                                                                                             cACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGaa                                                       <  2:1809087/102‑1 (MQ=255)
                                                                                                                                             cACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGaa                                                       >  1:1809087/1‑102 (MQ=255)
                                                                                                                                             cACCCATTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCAGCAGATGGCCCGATGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATGca       >  1:2020494/1‑150 (MQ=255)
                                                                                                                                                  ttttCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGTACGGCAGCGAACGCTTCGGCCACCCGTGGCAGGTTGGCCGGGCGCAGGCCGGACAGGCacacac  >  1:3383611/1‑150 (MQ=255)
                                                                                                                                                  ttttCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATGCacacac  >  1:3383602/1‑150 (MQ=255)
                                                                                                                                                  ttttCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATGCacacac  >  1:3383600/1‑150 (MQ=255)
                                                                                                                                                  ttttCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATGCacacac  >  2:1318477/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATGCACACAC  >  NC_002947/4081324‑4081619

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