Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,623,448 T→C N224D (AAC→GAC)  lplA ← lipoate‑‑protein ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,623,4480TC97.8% 138.5 / NA 46N224D (AAC→GAC) lplAlipoate‑‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (21/24);  minor base A (1/0);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TGCGCGGTACAGGCAGCCTTGCAGTCGTCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGC  >  NC_000913/4623203‑4623681
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          
cgcgcGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCgggg                                                                                                                                                                                                                                       <  1:471349/250‑1 (MQ=255)
                   tGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGGTTTGTCCGGGGAGATGATTTCCGCTTct                                                                                                                                                                                                                      <  1:396820/249‑2 (MQ=255)
                         cgcccGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGTGCGCGGGTGATATGGCCTTTTACAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCTCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTc                                                                                                                                                                                                              <  2:262607/248‑1 (MQ=255)
                           cccGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGc                                                                                                                                                                                                            >  1:225178/2‑251 (MQ=255)
                           cccGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGc                                                                                                                                                                                                            <  1:224276/250‑1 (MQ=255)
                                   gggCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATgg                                                                                                                                                                                                    >  2:56664/1‑251 (MQ=255)
                                                cgcgGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCAt                                                                                                                                                                                                         >  1:256753/1‑233 (MQ=255)
                                                cgcgGGGTTGAGGCTGGCGGTAAACCCCTGGGCGCGGGTGATATGGCCGTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCAt                                                                                                                                                                                                         <  2:256753/233‑1 (MQ=255)
                                                              agTCGGCAAACACCTGGGCGCGGGTGATATGGCATTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGACGGAGCCTGACCGAAGATCCATTCCCAGCTACTATGACGGGCAAAGGTTTCGGCGAAGTTTGGAACGTCTGGCGGTTTGTACGGGGGGGTGATTTCCGCTGCAAACCGCTGGCAATAATGGGGGAAAAAGGGCTTGGTTATGGCCttc                                                                                                                                                                        >  2:437968/2‑249 (MQ=255)
                                                               gTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGc                                                                                                                                                                       >  1:13864/1‑251 (MQ=255)
                                                                      aaCACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTg                                                                                                                                                                >  1:36253/1‑251 (MQ=255)
                                                                          ccTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTAGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTCCACGCTCGCCATAATGGGCGAAACAGGCCTTGGTTATGGCCTCGCAAACCTGCTc                                                                                                                                                             >  2:49457/1‑250 (MQ=255)
                                                                                  cgGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCAAGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGAt                                                                                                                                                    >  1:403962/1‑251 (MQ=255)
                                                                                       gATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGCTCCCCg                                                                                                                                               >  1:247402/1‑251 (MQ=255)
                                                                                                 ttttCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCt                                                                                                                                      <  1:40496/250‑1 (MQ=255)
                                                                                                  tttCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCg                                                                                                                                    >  1:457897/1‑251 (MQ=255)
                                                                                                        cGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAgg                                                                                                                              <  1:237245/251‑1 (MQ=255)
                                                                                                            gAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTgg                                                                                                                          <  2:276436/251‑1 (MQ=255)
                                                                                                               aTGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCa                                                                                                                       <  1:56664/251‑1 (MQ=255)
                                                                                                                   aGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCAcgcg                                                                                                                   >  1:138294/1‑251 (MQ=255)
                                                                                                                       ccACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGACTGCAGGAGCCTTACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGa                                                                                                                <  2:123353/250‑1 (MQ=255)
                                                                                                                        cACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAAGGCCGGAGCCTTACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTGCGGCGAAGTTTGGCAAGTCTGGCGGTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAAc                                                                                                              <  2:467799/251‑1 (MQ=255)
                                                                                                                              gccCCAGCTAAAGCTTTCAGCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCGGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACgcactg                                                                                                        <  1:274640/251‑6 (MQ=255)
                                                                                                                              gccCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGACGTTCCATTCCCAGATACTCTGACGGGCGAAGGTTTAGGCGAAGTTTGGCAAGTCTGGCGATTTGTCCGGGGAGATGATTTCGGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTGGCAAACCTGATCATGGGTGATCCACGCCAGCAAGTCGGTGAGGTTGGTCACGCGGGAACgcacgg                                                                                                        >  2:199888/1‑246 (MQ=255)
                                                                                                                                   agctAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACgcactgaggta                                                                                                   <  2:13864/248‑11 (MQ=255)
                                                                                                                                   agctAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACgcactgaggt                                                                                                    >  1:35495/4‑241 (MQ=255)
                                                                                                                                                    agcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAg                                                                                  <  2:35495/251‑1 (MQ=255)
                                                                                                                                                    agcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGGCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAg                                                                                  <  2:102134/251‑1 (MQ=255)
                                                                                                                                                      cagaggcaggaATGCCGGAGCCGGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGTTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTGGGCCGGGGAGATGATGTCCGCTTCGTCACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATTCCCGTCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCATTTCGCCGCCAGtt                                                                                <  2:433274/242‑1 (MQ=255)
                                                                                                                                                        gATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTATGGCGTTTTGTCCGGGGAGATGATTTACGCTTCTACACGCTCGACATAATGGGCGCAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGCGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGACTTTCGCCGCCAGtttt                                                                              >  1:278045/1‑250 (MQ=255)
                                                                                                                                                                      gtacccTGACCCAATTTCCACTCCCTGCGTCTCTTGCGGGCAAGGGTTTCGGCGAATTTTGGCAAGTCTGTCGTGTTGTCCGGTGAGATGATTTCCGCTTCTACACGCTCGCCATAATTGGCAAAAAGTGCCTTGGTTATTGCCTCGCAAACCTGCTCATTGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGTGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGAtt                                                                <  2:160321/247‑1 (MQ=255)
                                                                                                                                                                       gAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCtt                                                                                            <  2:101310/222‑1 (MQ=255)
                                                                                                                                                                       gAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCtt                                                                                            >  1:101310/1‑222 (MQ=255)
                                                                                                                                                                        aGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTg                                                               >  1:398472/1‑250 (MQ=255)
                                                                                                                                                                         gccgtaCCGAAGTTCCAGTCCCGGCGACTCTGACGGGCAAAGGTGTCGGCGAAGTTTGGCAAGTCTGGCGTTTGGTCCGGGGAGAGGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCTGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTgag                                                             <  2:247402/246‑1 (MQ=255)
                                                                                                                                                                           cTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTg                                                               >  1:382593/1‑247 (MQ=255)
                                                                                                                                                                                 gAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTg                                                     >  1:119055/1‑251 (MQ=255)
                                                                                                                                                                                              cggCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGTCAATCCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGGTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATTCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCtaag                                        <  2:417266/249‑4 (MQ=255)
                                                                                                                                                                                                gCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGt                                      <  1:359050/251‑1 (MQ=255)
                                                                                                                                                                                                           gggCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAgc                           <  1:23498/251‑2 (MQ=255)
                                                                                                                                                                                                               aaaGGTTTCGGCGAAGTTGGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAgcagca                       <  2:382593/251‑1 (MQ=255)
                                                                                                                                                                                                                            aaGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCgtggtg          >  1:94594/1‑251 (MQ=255)
                                                                                                                                                                                                                             aGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGtggtgg         >  1:257418/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tGGCAAGTCTGGCGTTTTTTCCGGGGAGATGATTTCCGCTTCTCCACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAATCCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAg      <  1:207973/250‑1 (MQ=255)
                                                                                                                                                                                                                                  ggCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGcc    >  2:454594/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGCGCAGCAAGGTGCCGTGGTGGAAGccac  >  2:231025/1‑249 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          
TGCGCGGTACAGGCAGCCTTGCAGTCGTCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGC  >  NC_000913/4623203‑4623681

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