Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 788,581 C→A 5.2% intergenic (+31/‑182) ettA → / → gdhA ADP/ATP ratio sensor/translation inhibitor/NADP‑dependent glutamate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947788,5810CA 5.2% 295.8 / 4.6 96intergenic (+31/‑182)ettA/gdhAADP/ATP ratio sensor/translation inhibitor/NADP‑dependent glutamate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (51/40);  new base A (3/2);  total (54/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-01

CTCGAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCC  >  NC_002947/788433‑788728
                                                                                                                                                    |                                                                                                                                                   
cTCGAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGccc                                                                                                                                                     >  3:223409/1‑149 (MQ=255)
       gtggtgTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTgg                                                                                                                                              <  3:390417/149‑1 (MQ=255)
        tggtgTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAACTGGc                                                                                                                                             >  8:76270/1‑149 (MQ=255)
          gtgtTTTTAGATTGTAACTACACCGAGTACGAAGGCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGccc                                                                                                                                           <  6:280778/149‑1 (MQ=255)
            gtTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCa                                                                                                                                         >  3:190355/1‑149 (MQ=255)
             ttcttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGGGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAt                                                                                                                                        <  7:16530/149‑1 (MQ=255)
               cttcGAAGGTAACTACACCCACTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATTGGGCCCCAAGTGGCCCCAttt                                                                                                                                      >  2:322640/1‑149 (MQ=255)
                ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAtttt                                                                                                                                     >  3:230238/1‑149 (MQ=255)
                  cGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAATCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttttt                                                                                                                                    <  8:321541/148‑1 (MQ=255)
                  cGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttttt                                                                                                                                    <  8:86006/148‑1 (MQ=255)
                  cGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAtttt                                                                                                                                     <  8:398433/147‑1 (MQ=255)
                  cGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTtg                                                                                                                                   >  8:82998/1‑149 (MQ=255)
                  cGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAAAATGGGGCCACTTGGGGCCCCAttc                                                                                                                                      <  5:359231/146‑2 (MQ=255)
                       gTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAAAATGGGGCCACTTGGGGCCCCATTCTTGTGcat                                                                                                                              <  5:342523/149‑3 (MQ=38)
                          aCTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGc                                                                                                                           <  7:18916/149‑1 (MQ=255)
                           cTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCg                                                                                                                          <  6:403267/149‑1 (MQ=255)
                             acacCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttt                                                                                                                                      >  1:367664/1‑135 (MQ=255)
                             acacCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttt                                                                                                                                      <  2:367664/135‑1 (MQ=255)
                                 cGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGc                                                                                                                    >  5:127468/1‑149 (MQ=255)
                                       cGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATaaa                                                                                                              >  8:254131/1‑149 (MQ=255)
                                       cGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATaa                                                                                                               <  4:188222/148‑1 (MQ=255)
                                                 cGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcgaggc                                                                                                     >  4:241064/1‑148 (MQ=255)
                                                 cGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGAGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGAGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcgaggc                                                                                                     >  7:71032/1‑148 (MQ=255)
                                                   cAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGaa                                                                                                  <  6:127468/149‑1 (MQ=255)
                                                    aagaagCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAg                                                                                                 <  1:171852/149‑1 (MQ=255)
                                                    aagaagCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAg                                                                                                 <  3:101706/149‑1 (MQ=255)
                                                    aagaagCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAg                                                                                                 <  7:82998/149‑1 (MQ=255)
                                                    aagaagCGGATGGGTGCTGCTTCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAg                                                                                                 <  8:180674/149‑1 (MQ=255)
                                                             cTGTGTGATGCTGCTGGCCAGCCGGACGGTGTACGGAACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACATGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACGGCTATAAATGGCGAGGCGAAGGTGTAtttt                                                                                        <  2:341043/149‑1 (MQ=255)
                                                                ggTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTAtttt                                                                                        <  3:241064/146‑1 (MQ=255)
                                                                ggTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTAtttt                                                                                        <  8:110467/146‑1 (MQ=255)
                                                                ggTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGta                                                                                     <  7:254131/149‑1 (MQ=255)
                                                                    aTGCTGCTGCCCGGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATac                                                                                  >  8:338382/1‑148 (MQ=255)
                                                                        tgctgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCtat                                                                                                                 >  5:195060/1‑113 (MQ=255)
                                                                        tgctgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCtat                                                                                                                 <  6:195060/113‑1 (MQ=255)
                                                                                aGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAAc                                                                     >  1:118004/1‑149 (MQ=255)
                                                                                aGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAAc                                                                     >  1:162530/1‑149 (MQ=255)
                                                                                  ccGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATaaaa                                                                      >  2:222080/1‑146 (MQ=255)
                                                                                  ccGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATaaaa                                                                      <  1:222080/146‑1 (MQ=255)
                                                                                  ccGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAAcac                                                                   >  5:7805/1‑147 (MQ=255)
                                                                                    gCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCAc                                                                 <  4:223409/149‑1 (MQ=255)
                                                                                     cACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCAc                                                                 >  3:417050/1‑148 (MQ=255)
                                                                                     cACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCAc                                                                 >  4:142765/1‑148 (MQ=255)
                                                                                         gtgtACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATacac                                                                                >  6:107929/1‑129 (MQ=255)
                                                                                         gtgtACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATacac                                                                                <  5:107929/129‑1 (MQ=255)
                                                                                            tACGGCACAAGAAACTGGCCCAGTAAGCAGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAaataat                                                         >  6:128873/1‑149 (MQ=255)
                                                                                               ggCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAAttt                                                       <  7:35487/148‑1 (MQ=255)
                                                                                               ggCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAAttt                                                       >  8:35487/1‑148 (MQ=255)
                                                                                                 cacaAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCa                                                     >  8:66786/1‑148 (MQ=255)
                                                                                                  tcaAGAAACTGGCACAGTAAGCCGGTTTTCTGATGCAAAAGAATGGGGCCCCAAGTGGACCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaa                                                    <  4:372969/147‑1 (MQ=255)
                                                                                                  acaAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaa                                                    >  4:377980/1‑148 (MQ=255)
                                                                                                  acaAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaa                                                    >  1:340201/1‑148 (MQ=255)
                                                                                                    aaGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcg                                                                                                         <  8:205813/93‑1 (MQ=255)
                                                                                                    aaGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcg                                                                                                         >  7:205813/1‑93 (MQ=255)
                                                                                                    aaGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaaaa                                                  >  3:291882/1‑148 (MQ=255)
                                                                                                    aaGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCACCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATTGCTAGGCGAAGGTGTATTTTGTATACACTTACAAAAAACGCACCAAATAATTTCaa                                                    >  3:372969/1‑146 (MQ=255)
                                                                                                     aGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATaaaa                                                                      >  1:75956/1‑127 (MQ=255)
                                                                                                     aGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATaaaa                                                                      <  2:75956/127‑1 (MQ=255)
                                                                                                      gAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCGCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACCATTACATAAAACGCACGAAATAATTTCaaaaaaa                                               >  8:349652/1‑149 (MQ=255)
                                                                                                      gAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCTAGGCGAAGGTGTATTTTGGATACACTGACATAAAAAGCACCAAAGAATTTAaaaaaa                                                >  8:355760/1‑148 (MQ=255)
                                                                                                       aaaCTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAAc                                              >  1:29745/1‑149 (MQ=255)
                                                                                                           tGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGTCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACa                                          >  1:81185/1‑149 (MQ=255)
                                                                                                               ccAGTAAGCCGGTTTTCTGATGCACAAGATTGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGta                                                                                     <  6:292580/102‑1 (MQ=255)
                                                                                                               ccAGTAAGCCGGTTTTCTGATGCACAAGAGTGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGta                                                                                     >  5:292580/1‑102 (MQ=255)
                                                                                                                cAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCg                                     >  2:148925/1‑149 (MQ=255)
                                                                                                                 aGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGc                                    >  5:322056/1‑149 (MQ=255)
                                                                                                                 aGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAAAGCACCAAATAATTTCAAAAAAACGACATTTAGc                                    >  4:41608/1‑149 (MQ=255)
                                                                                                                  gTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGcc                                   <  7:76270/149‑1 (MQ=255)
                                                                                                                  ctAAGCCGGTTTTCTGATCCCCAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGATGCGAAGGTGTATTATGTATACACTTACATAAAACTCACCAAATAATTTCAAAAAAACGACATTTCGcc                                   <  4:302218/148‑1 (MQ=255)
                                                                                                                      gCCGTTTTTCTGTTGCACATGAATGGGGGCCCAAGTGGGACAATTTTTGTGCCTGGCGGACCGCTATAAATGGAGATGCGAATGTGTTTTTTGTATAAAGTTACATAAAACGCACCAAATAATTTCAAAAAAACGTCATTTCGCCCt                                 <  4:238831/147‑1 (MQ=255)
                                                                                                                      gCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCt                                 >  3:238831/1‑147 (MQ=255)
                                                                                                                      gCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATATCAAAAAAAAGACATTTCGCCCTAt                               >  2:209183/1‑149 (MQ=255)
                                                                                                                           ttttCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCgg                          >  5:208865/1‑149 (MQ=255)
                                                                                                                                tGATGCACAAGAATGGGGCCCCAAGTGGCCCTATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGa                     >  5:155181/1‑149 (MQ=255)
                                                                                                                                tGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGa                     >  1:244670/1‑149 (MQ=255)
                                                                                                                                tGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGa                     >  4:206174/1‑149 (MQ=255)
                                                                                                                                  aTGCACAAGAATGGGGCCCCAATTGGCCACATTTTTGTGCATGGCGGACCACTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGAt                    <  1:144543/148‑1 (MQ=255)
                                                                                                                                  aTGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGAt                    <  4:291882/148‑1 (MQ=255)
                                                                                                                                   ggcgcaAAAATGGGGCCACTTGGGGCCCCATTCTTGTGACTGGCGGACCGCTATAAATTGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAAAGCACCAAATAATTTCAAAAAAACGACATTGCGCACTATTCCGGTGAGATTg                  >  6:262393/5‑149 (MQ=37)
                                                                                                                                    gCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGc                 >  1:17467/1‑149 (MQ=255)
                                                                                                                                    gCACAAAAATGGGGCCACTTGGGGCCCCATTCTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGc                 >  7:371112/1‑149 (MQ=39)
                                                                                                                                    gCACAAAAATGGGGCCACTTGGGGCCCCATTCTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGc                 >  7:396448/1‑149 (MQ=39)
                                                                                                                                     cacaATAATGTGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATATCGCCCTATTCCGGTGCGATTGc                 <  7:66786/148‑1 (MQ=255)
                                                                                                                                     cacaAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCAc                                                                 <  2:28881/100‑1 (MQ=255)
                                                                                                                                     cacaAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCAc                                                                 >  1:28881/1‑100 (MQ=255)
                                                                                                                                     cacaAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGct                >  4:358757/1‑149 (MQ=255)
                                                                                                                                     cacaAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGc                 <  3:377980/148‑1 (MQ=255)
                                                                                                                                      acaAGAATTGGGCCCCAAGTGGCCCCAATTTTGTGACTAGCGCACCGCTATAAATGGCGAGGCGAAGGTCTATTTTGTATACAATTACATAAAACGCACAAAAAAATTTCAAAAAAACCACATTTCGCCCTATTCAGGTCCGATTGct                >  4:180366/1‑148 (MQ=255)
                                                                                                                                      acaAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAATTAATTTCAAAAAAACGTCATTTCGCCGTATTCAGGTGCGATTGctt               >  3:290808/1‑149 (MQ=255)
                                                                                                                                      acaAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGctt               <  2:118004/149‑1 (MQ=255)
                                                                                                                                        aaGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCGCCAAATAATTTCAAAAAAACGACGTTTCGCCGTAGTCGGGTGCGATGGgttcg             >  4:47583/1‑148 (MQ=255)
                                                                                                                                        aaGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttct             <  3:47583/149‑1 (MQ=255)
                                                                                                                                        aaGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttct             >  3:429635/1‑149 (MQ=255)
                                                                                                                                        aaGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttct             <  6:208865/149‑1 (MQ=255)
                                                                                                                                          gAATGTGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCTGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTg           <  1:33531/149‑1 (MQ=255)
                                                                                                                                             tGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGt         >  5:111900/1‑148 (MQ=255)
                                                                                                                                                ggCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCTCTATAAATGGCGAGGCGAAGGTGTATTTTTTATACACTTACATAAAACGAACCAAATAATTTCAAAAAAACGACATTTAGCCCTATTCCTGTGCAATTGCTTCTTGGTACAt     >  6:323539/1‑149 (MQ=255)
                                                                                                                                                  ccccAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTc   >  3:321616/1‑149 (MQ=255)
                                                                                                                                                   cccAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTcc  <  4:190355/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTCGAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCC  >  NC_002947/788433‑788728

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