Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 126,032 T→G I113S (ATC→AGC)  aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913126,0320TG98.3% 200.4 / ‑3.6 58I113S (ATC→AGC) aceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (24/33);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

CCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAG  >  NC_000913/125889‑126171
                                                                                                                                               |                                                                                                                                           
ccAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGc                                                                                                                                      <  2:650544/151‑1 (MQ=255)
ccAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGc                                                                                                                                      >  2:316719/1‑151 (MQ=255)
ccAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGc                                                                                                                                      <  1:390321/151‑1 (MQ=255)
        gcgcACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGt                                                                                                                              <  1:388986/151‑1 (MQ=255)
         cgcACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGtt                                                                                                                             <  2:179305/151‑1 (MQ=255)
            aCTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGa                                                                                                                           >  2:363693/1‑150 (MQ=255)
             cTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCGCCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGa                                                                                                                           >  2:686875/1‑149 (MQ=255)
               gATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTg                                                                                                                       <  1:425645/151‑1 (MQ=255)
                  tATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGAcc                                                                                                                    >  2:531909/1‑151 (MQ=255)
                      aTTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGaa                                                                                                                 >  1:57684/1‑150 (MQ=255)
                         ttCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATc                                                                                                              <  1:453124/150‑1 (MQ=255)
                          tCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCt                                                                                                            >  1:394271/1‑151 (MQ=255)
                           cGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCt                                                                                                            <  1:377991/150‑1 (MQ=255)
                            gATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTgg                                                                                                          <  1:805688/151‑1 (MQ=255)
                             aTTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGt                                                                                                         <  1:208059/151‑1 (MQ=255)
                               tCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTg                                                                                                        <  1:311230/150‑1 (MQ=255)
                                  gccgACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGaa                                                                                                      >  2:772568/1‑149 (MQ=255)
                                     gACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTg                                                                                                 >  2:585524/1‑151 (MQ=255)
                                     gACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTg                                                                                                 <  1:282802/151‑1 (MQ=255)
                                       cGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGc                                                                                               <  2:651370/151‑1 (MQ=255)
                                           gcagcaGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGAt                                                                                            <  1:85203/150‑1 (MQ=255)
                                           gcagcaGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATa                                                                                           <  2:589601/151‑1 (MQ=255)
                                              gcagACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATaa                                                                                          >  1:926109/1‑149 (MQ=255)
                                               cagACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGt                                                                                       <  2:457100/151‑1 (MQ=255)
                                                agACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGtt                                                                                      >  1:438364/1‑151 (MQ=255)
                                                 gACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTg                                                                                     <  2:205576/151‑1 (MQ=255)
                                                   cGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGATGGCGATAAAGTTGa                                                                                    >  1:891400/1‑150 (MQ=255)
                                                        cACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGa                                                                              <  2:736982/151‑1 (MQ=255)
                                                         aCCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCGGCGACGAAGTTGAAGTGACCGAAATCCTGGAGAAAGGTGGCGATAAAGTTGAAGCTGa                                                                              >  1:630626/1‑150 (MQ=255)
                                                         aCCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGaa                                                                             <  1:825167/151‑1 (MQ=255)
                                                           cTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATGCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAAc                                                                            >  1:199419/1‑150 (MQ=255)
                                                             gCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACa                                                                           <  2:99908/149‑1 (MQ=255)
                                                               tCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTc                                                                        >  2:654095/1‑150 (MQ=255)
                                                                cAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTc                                                                        <  1:684364/149‑1 (MQ=255)
                                                                      gaagagaagaAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCa                                                                <  2:49976/151‑1 (MQ=255)
                                                                              gaaagaaGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGaa                                                        <  2:3134/151‑1 (MQ=255)
                                                                                  gaagCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCg                                                    <  2:389950/151‑1 (MQ=255)
                                                                                    agcagcTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGAc                                                  >  2:305243/1‑151 (MQ=255)
                                                                                     gcagcTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGAc                                                  >  2:887161/1‑150 (MQ=255)
                                                                                           ccGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGc                                             <  1:903119/149‑1 (MQ=255)
                                                                                            cGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTc                                          >  2:490954/1‑151 (MQ=255)
                                                                                             ggCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCt                                         >  2:582617/1‑151 (MQ=255)
                                                                                                   agcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGa                                    >  2:799377/1‑150 (MQ=255)
                                                                                                   agcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGa                                    >  2:512388/1‑150 (MQ=255)
                                                                                                          gCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCgg                            >  1:238613/1‑151 (MQ=255)
                                                                                                          gCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCgg                            >  2:879409/1‑151 (MQ=255)
                                                                                                           cGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGc                           >  2:882961/1‑151 (MQ=255)
                                                                                                                gcggcgGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTcc                        <  2:578752/149‑1 (MQ=255)
                                                                                                                gcggcgGCAAAAGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGt                      <  2:130632/151‑1 (MQ=255)
                                                                                                                         aaaGACGTTAACGTTCCGGATAGCGGGAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCa             <  2:698256/151‑1 (MQ=255)
                                                                                                                         aaaGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTCACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCa             <  1:565159/151‑1 (MQ=255)
                                                                                                                           aGACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAcc           <  2:256707/151‑1 (MQ=255)
                                                                                                                            gACGTTAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCg          >  1:87649/1‑151 (MQ=255)
                                                                                                                                 tAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGaa      <  1:542972/150‑1 (MQ=255)
                                                                                                                                 tAACGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGaa      <  1:363679/150‑1 (MQ=255)
                                                                                                                                  aaCGTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCTCATTCCTGGTGAAAGTTGGCGATAAACTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACAGTGaaa     <  1:138422/150‑1 (MQ=18)
                                                                                                                                     gTTCGGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAGATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAagag  <  1:304910/150‑1 (MQ=255)
                                                                                                                                     gTTCCGGATAGCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAag    >  2:11733/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                           
CCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAG  >  NC_000913/125889‑126171

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