Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,037 T→C Q27Q (CAA→CAG)  lplA ← lipoate‑‑protein ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,624,0370TC97.4% 115.3 / ‑4.4 39Q27Q (CAA→CAG) lplAlipoate‑‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (15/23);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

AGCGCATTGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAG  >  NC_000913/4623792‑4624277
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
aGCGCATTGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCg                                                                                                                                                                                                                                              <  1:195386/251‑1 (MQ=255)
     aTTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACGCTTTCCACGGGTTCTGCGCTCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGaaaaa                                                                                                                                                                                                                                         <  1:126916/251‑1 (MQ=255)
      ttGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCGGGCGAAAAAt                                                                                                                                                                                                                                        >  2:189357/1‑251 (MQ=255)
        gAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATac                                                                                                                                                                                                                                      <  1:193030/251‑1 (MQ=255)
                 cGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAGTACACtcctccc                                                                                                                                                                                                                             >  2:256462/1‑250 (MQ=255)
                  gACGTGGAGATGGTTTTAGCGTACGCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTGTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGTCGAAAAATACTCTCCTCCAc                                                                                                                                                                                                                            <  2:63248/251‑1 (MQ=255)
                              gTTTTATCGTACGCCTGCTGGCCAGCCAGAAAGGTAAAGCAGGTATGGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTaaa                                                                                                                                                                                                                <  2:25027/251‑1 (MQ=255)
                                        acgccGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCg                                                                                                                                                                                                      <  1:256462/248‑1 (MQ=255)
                                          cccGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTa                                                                                                                                                                                                    <  1:223241/250‑1 (MQ=255)
                                           cccgCTTGCCAGCCATAAAGGTAAAGCAGGGATTGCCGAGATCGGGGAACACCGCGCCGCCCCCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCAGCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACGCCTCCTCCTCCAGGTTAAACCACGGGTCGt                                                                                                                                                                                                     <  2:37029/246‑1 (MQ=255)
                                           ccGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGt                                                                                                                                                                                                     <  1:33938/249‑1 (MQ=255)
                                             ggCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACCCCGCGCCGCCACCGCTACTGCGCCGTGACAGGCGGACGTTATCTTCTTCCATCCGCCGGGCATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCTCTACCGTGTCGGCATTGCGCCAGCGCAACCGATCGCGCTGCGTGCCGGGCATCTGGCGAAAAATACACTCTTCCACCGCAAGGTTATACCACgtggctta                                                                                                                                                                                                    >  1:230059/1‑241 (MQ=255)
                                              gCTTGCCAGCCAGAAAGGTAAAACAGGTAGTGCCGAGATCGTGGAACCCCGCGCCGCCACCGCTACGGCGCCGTGCCAGGCGGACGTTATCTTCTTCCCTCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAagag                                                                                                                                                                                                <  2:151641/251‑1 (MQ=255)
                                                  gccagccaGAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTTCTGCTCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCa                                                                                                                                                                                             <  2:65861/250‑1 (MQ=255)
                                                             aGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAgg                                                                                                                                                                                >  2:205816/1‑251 (MQ=255)
                                                               gTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCCCCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCATAGATG‑AGCAGGCg                                                                                                                                                                              <  1:201172/251‑1 (MQ=255)
                                                                            tGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACAt                                                                                                                                                                  >  1:77225/1‑250 (MQ=255)
                                                                             gCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATaa                                                                                                                                                                <  1:130753/251‑1 (MQ=255)
                                                                                             cccGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgt                                                                                                                                                <  1:240085/250‑1 (MQ=255)
                                                                                             ccagcgcCGCCACCGCTCCGGCGCCGTTCCAGGCGGACGTTAACTTCTGCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTGCCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCAGCTGGCGAAAAATACACTCCGCCACCGCTAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCATGCGTAATGTGGACATAACGATTTCCTTGCTTgt                                                                                                                                                <  2:264601/248‑1 (MQ=255)
                                                                                                        aCCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACCGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGACAGAGATGCACCAG‑CGTAATGTGGCCATAACGATTTCCTTTCTTGCAGTAGTGGGt                                                                                                                                      >  2:8567/1‑250 (MQ=255)
                                                                                                           gCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGGGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGTTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCAct                                                                                                                                  <  1:237355/251‑1 (MQ=255)
                                                                                                            cTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCActc                                                                                                                                 >  1:166925/1‑251 (MQ=255)
                                                                                                                     cgtgccAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTtctc                                                                                                                        >  1:247931/4‑251 (MQ=255)
                                                                                                                           aGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCttttaaaagagtgac                                                                                                                  <  1:89182/251‑12 (MQ=255)
                                                                                                                                          tcttcCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCa                                                                                                    <  1:205816/250‑1 (MQ=255)
                                                                                                                                             tcCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCAcc                                                                                                  <  1:149131/249‑1 (MQ=255)
                                                                                                                                                   cgccgGGTATTACACTCTTTCAACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCCCCGGTTTGCt                                                                                          >  2:59681/1‑251 (MQ=255)
                                                                                                                                                              acacacTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGCAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAATTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGCAATGGGTCACTCTTTTTTTTCGATTTCTTCCTCTTCCGGCACCGGGTGTCTGGCGGtttac                                                                                >  2:150506/1‑247 (MQ=255)
                                                                                                                                                                   ctTTCCACGGGTTCTGCGCGCGACCAATCACTGCCTTGTCGGCAGTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAATCCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGTGTCACTCGTTTTACTCGCTGTCTTCCGCTTCCGGCACCGGTTTGCTGGCGTTTACCAGGa                                                                            <  2:136126/249‑1 (MQ=255)
                                                                                                                                                                                        gACCAATCACTTCCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTTGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGta                                                     <  2:195214/251‑2 (MQ=255)
                                                                                                                                                                                                 ctacCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTTACAGGGGGGGCGATTGCTGCCAGCGGGGCCGTTTAccc                                            >  2:219822/1‑251 (MQ=255)
                                                                                                                                                                                                   acCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAGAGGCGCTTGCTGCCAGCGGGGACGTTTATCCTg                                          >  2:219282/1‑251 (MQ=255)
                                                                                                                                                                                                          cGGCATTGCGACAGAGAAACAGAACGCGCTGCGGGGCGGGCATCTGGCGAAAAATACACTACCCCACCGCCAGGGTAAAACACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTGACCCTGTagcag                                    <  2:226647/250‑1 (MQ=255)
                                                                                                                                                                                                                 gcgcCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTTTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCgg                             <  1:120080/250‑1 (MQ=255)
                                                                                                                                                                                                                 gcgcCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCgg                             <  1:267869/250‑1 (MQ=255)
                                                                                                                                                                                                                   gcCAGCGAAACAGAACGGCCTGCCTGGCGGGCATCTGGCGTAAAATACACTCCTCCCCCGCCCGGTTAAACCACGGGTCGTGAGAGTCAGAGATG‑CGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGCAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACCGGGAAGGCGCTTACTGCCCGCGGGTAAGTTTAccccatcgccgcaggcgcgtt                          >  2:20263/1‑233 (MQ=255)
                                                                                                                                                                                                                        agaAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGGGGGCGGGTccgccc                     >  2:266490/1‑246 (MQ=255)
                                                                                                                                                                                                                                              cGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGGTTGcccgggaa  >  2:178166/1‑242 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
AGCGCATTGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAG  >  NC_000913/4623792‑4624277

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑