Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,364,099 A→T 6.2% intergenic (‑223/+110) PP_2076 ← / ← lysA‑I hypothetical protein/diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,364,0990AT 6.2% 240.3 / 4.7 80intergenic (‑223/+110)PP_2076/lysA‑Ihypothetical protein/diaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (36/38);  new base T (3/2);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

TCTGATTAAGCGGTATCGTTGACGTTGGCAGGCTTTTCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTAGAAATCCGCCTCCAGGCTCAGCAGGTCCTGCAC  >  NC_002947/2363953‑2364244
                                                                                                                                                   |                                                                                                                                                 
tCTGATTAAGCGGTATCGTTGACGTTGGCAGGCTTTTCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCaaa                                                                                                                                                  <  5:328483/148‑1 (MQ=255)
tCTGATTAAGCGGTATCGTTGACGTTGGCAGGCTTTTCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCaaa                                                                                                                                                  <  5:322788/148‑1 (MQ=255)
 cTGATTAAGCGGTATCGTTGACGTTGGCAGGCTTTTCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGc                                                                                                                                                >  6:336752/1‑149 (MQ=255)
     gtAAGCGGTATCGTTGACGTTGGCAGGCTTTTCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGTGGcc                                                                                                                                            <  2:345464/148‑1 (MQ=255)
        ggCGGTATCGTTGACGTTGGCAGGCTTTTCGTCGCCTTTCACGTCACCCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATAGAGACAGGGCCGCACAGCGGcccc                                                                                                                                          >  5:244433/2‑147 (MQ=255)
          cGGTATCGTTGACGTTGGCAGGCTTTTCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGAtc                                                                                                                                       <  6:61222/149‑1 (MQ=255)
           ggTATCGTTGACGTTGGCAGGCTTTTCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACGCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGAGAGGGCCGCAAAGCGGCCCGAtct                                                                                                                                      <  2:267661/149‑1 (MQ=255)
              aTCGTTGACGTTGGCAGGCTTTTCGTCGCCTTTCACGTCAACCCAGCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCAc                                                                                                                                   >  1:262989/1‑149 (MQ=255)
                 gTTGACGTTGGCAGGCTTTTCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCa                                                                                                                                <  4:352034/149‑1 (MQ=255)
                 gTTGACGTTGGCAGGCTTTTAGTCGCCTTTCACGTCAACCCACGTGCCTCTCGTTGCGCCAGCGCTATCGGGTGCAGAGCCAACCAGGTCGCTTTTTATAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCAcc                                                                                                                                  >  5:67125/1‑147 (MQ=255)
                  ttGACGTTGGCAGGCTTTTCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGATCGGGCCGCTTTGCGGccctgtc                                                                                                                                       >  8:116318/1‑137 (MQ=255)
                  ttGACGTTGGCAGGCTTTTCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGATCGGGCCGCTTTGCGGccctgtc                                                                                                                                       <  7:116318/141‑5 (MQ=255)
                          ggCAGGCTTTTCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTGTTTAATATGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGc                                                                                                                        <  4:363975/148‑1 (MQ=255)
                            cAGGCTTTTCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTg                                                                                                                     >  4:314025/1‑149 (MQ=255)
                             aGGCTTTTCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCTTTGCGGccctgt                                                                                                                                        >  5:13239/1‑126 (MQ=255)
                             aGGCTTTTCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCTTTGCGGccctgt                                                                                                                                        <  6:13239/129‑4 (MQ=255)
                             aGGCTTTTCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCCCAAAGCCGCCCGAAGTCACCCCCCGTGAGCGGGc                                                                                                                    >  2:161190/1‑149 (MQ=255)
                                 ttttCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCggg                                                                                                                >  3:178589/1‑149 (MQ=255)
                                 ttttCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCCGAAAGCGGCCCGATCTAACACAGCTTGCGCGTGCCggg                                                                                                                >  6:21349/1‑149 (MQ=255)
                                   ttCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCgggg                                                                                                               <  1:361109/148‑1 (MQ=255)
                                   ttCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCggg                                                                                                                <  1:279766/147‑1 (MQ=255)
                                      gtcgCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGAt                                                                                                            <  6:364999/148‑1 (MQ=255)
                                            tttCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACGCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCt                                                                                                      >  1:102039/1‑149 (MQ=255)
                                            tttCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTg                                                                                                     >  1:335131/1‑149 (MQ=255)
                                            tttCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTg                                                                                                     >  1:169373/1‑149 (MQ=255)
                                            tttCACGTCAACCCACCTGCCTCTCGGTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTg                                                                                                     >  1:5303/1‑149 (MQ=255)
                                              tCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTgg                                                                                                    <  2:262989/148‑1 (MQ=255)
                                                  gTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCg                                                                                                                                               <  6:269960/101‑1 (MQ=255)
                                                  gTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCg                                                                                                                                               >  5:269960/1‑101 (MQ=255)
                                                  gTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTg                                                                                               <  1:161190/149‑1 (MQ=255)
                                                  gTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCCAACTCACCCAGCGTGAGCGGGGCGGGGGATGAAGCTGGAGGTg                                                                                               >  4:366409/1‑149 (MQ=255)
                                                    cAACCCACCTGGCTCCGCTTTCGGCACGGGCATCGGGTGCATTTCCAACCTGGTCGCTTTTTTTAATACGGTTGCCCATCGTGCAATGGCCGCAAACCGGCCCGTTCTCACCCATCTTTTGCGGTGCGGGGGATGAAGCTGGAGGTGaa                                                                                             <  4:129875/149‑1 (MQ=255)
                                                     aaCCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGaa                                                                                             <  2:5303/148‑1 (MQ=255)
                                                     aaCCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAg                                                                                            <  6:143197/149‑1 (MQ=255)
                                                            cTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGtt                                                                                     <  5:94188/149‑1 (MQ=255)
                                                                            cAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTAACCCAGCTTGAGCGTGCAGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGccc                                                                     >  6:204188/1‑149 (MQ=255)
                                                                                   aTCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAag                                                              >  3:23438/1‑149 (MQ=255)
                                                                                    tCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCaaaa                                                               <  2:335131/147‑1 (MQ=255)
                                                                                    tCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAaga                                                             <  5:21349/149‑1 (MQ=255)
                                                                                      gggTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAgt                                                           <  3:208922/149‑1 (MQ=255)
                                                                                      gggTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAgt                                                           <  3:314025/149‑1 (MQ=255)
                                                                                       ggTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAgtg                                                          >  1:21979/1‑149 (MQ=255)
                                                                                        gTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGTCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAgtgt                                                         >  3:248587/1‑149 (MQ=255)
                                                                                        gTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAgtgt                                                         <  3:366409/149‑1 (MQ=255)
                                                                                           cAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGc                                                       <  2:169373/148‑1 (MQ=255)
                                                                                            agagCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAg                                                                                            >  7:26483/1‑110 (MQ=255)
                                                                                            agagCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAg                                                                                            <  8:26483/110‑1 (MQ=255)
                                                                                                ccACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGcc                                                  >  7:145778/1‑148 (MQ=255)
                                                                                                 cACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACCGGGCCGCAAAGCGGCCCGATCTCACCCAGCGTGAGCGTGCCGGGGGATGAAGCGGGGGGGGGAGGGCAGTTGGAAGCGGGCGGGGCCGCAAAAGAGTGTgcccgcccc                                                >  4:42319/1‑148 (MQ=255)
                                                                                                     ccTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAggg                                             >  4:326684/1‑149 (MQ=255)
                                                                                                       tGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCt                                                                                                      >  1:337718/1‑90 (MQ=255)
                                                                                                       tGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCt                                                                                                      <  2:337718/90‑1 (MQ=255)
                                                                                                       tGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGc                                           >  3:309198/1‑149 (MQ=255)
                                                                                                       ggggcgcgTTTTTTAATTCCGCGCCCCAGCGTGACAGGGCCGCAAAGCGGCCCGATCTCACCCTGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGATCGCCAGTTGGAACCGGGCGTGCCCGCAAAACTGTGTGCACGCCCAGGGGc                                           <  6:67125/145‑1 (MQ=255)
                                                                                                         gTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGa                                         >  3:132414/1‑149 (MQ=255)
                                                                                                          tCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGAt                                        >  3:181691/1‑149 (MQ=255)
                                                                                                          tCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGAt                                        <  2:21979/149‑1 (MQ=255)
                                                                                                            gCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCt                                      >  4:51485/1‑149 (MQ=255)
                                                                                                             cTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCtt                                     <  3:157292/149‑1 (MQ=255)
                                                                                                              tttttttAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAgt                                                           <  8:18839/126‑1 (MQ=255)
                                                                                                              tttttttAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAgt                                                           >  7:18839/1‑126 (MQ=255)
                                                                                                              tttttttAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCttt                                    >  2:169590/1‑149 (MQ=255)
                                                                                                              tttttttAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCtt                                     >  5:226448/1‑148 (MQ=255)
                                                                                                                tttttAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCttt                                    <  4:181691/147‑1 (MQ=255)
                                                                                                                tttttAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCtt                                     <  6:226448/146‑1 (MQ=255)
                                                                                                                tttttAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTa                                   <  4:202616/148‑1 (MQ=255)
                                                                                                                tttttAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTa                                   <  4:309198/148‑1 (MQ=255)
                                                                                                                 ttttAATACGGTGCCGCATCGAGACAGGGTCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGTGGGATGAAGCTGGAGGTGAATGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTAg                                  <  4:248587/148‑1 (MQ=255)
                                                                                                                     aaTACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCGAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTAGAAAt                              <  2:102039/148‑1 (MQ=255)
                                                                                                                         cGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTAGAAATCCGc                          <  4:178589/148‑1 (MQ=255)
                                                                                                                            tGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTAGAAATCCGCCTc                       >  8:306868/1‑148 (MQ=255)
                                                                                                                              ccGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTAGAAATCCGCCTCCa                     <  5:336752/148‑1 (MQ=255)
                                                                                                                                  aTCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTAGAAATCCGCCTCCAGGCt                 >  4:341654/1‑148 (MQ=255)
                                                                                                                                     gagaCAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTa                                   >  7:12302/1‑127 (MQ=255)
                                                                                                                                     gagaCAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTa                                   <  8:12302/127‑1 (MQ=255)
                                                                                                                                     gagaCAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTAGAAATCCGCCTCCAGGCTcagc             <  1:169590/149‑1 (MQ=255)
                                                                                                                                      agatcgggCCGCTTTGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTAGAAATCCGCCTCCAGGCTcag              >  7:111597/6‑147 (MQ=255)
                                                                                                                                      agaCAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTAGAAATCCGCCTCCAGGCTcag              >  5:342348/1‑147 (MQ=255)
                                                                                                                                        aCAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTAGAAATCCGCCTCCAGGCTCAGCAgg          <  4:23438/149‑1 (MQ=255)
                                                                                                                                           gggCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTAGAAATCCGCCTCCAGGCTCAGCAGGTcc       >  5:155120/1‑149 (MQ=255)
                                                                                                                                               cGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCTGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTAGAAATCCGCCTCCAGGCTCAGCAGGTCCTGCa   >  2:108554/1‑149 (MQ=255)
                                                                                                                                                gCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTAGAAATCCGCCTCCAGGCTCAGCAGGTCCTGCAc  >  2:77584/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCTGATTAAGCGGTATCGTTGACGTTGGCAGGCTTTTCGTCGCCTTTCACGTCAACCCACCTGCCTCTCGTTGCGCCAGGGCTATCGGGTGCAGAGCCACCCTGGTCGCTTTTTTTAATACGGTGCCGCATCGAGACAGGGCCGCAAAGCGGCCCGATCTCACCCAGCTTGAGCGTGCCGGGGGATGAAGCTGGAGGTGAAGGCCAGTTGGAAGCGGGCGTGCCCGCAAAAGAGTGTGCACGCCCAGGGGCGATCTTTAGAAATCCGCCTCCAGGCTCAGCAGGTCCTGCAC  >  NC_002947/2363953‑2364244

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