Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 125,724 C→T 5.7% I10I (ATC→ATT)  aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913125,7240CT 5.7% 228.6 / 2.4 70I10I (ATC→ATT) aceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base C (28/38);  new base T (2/2);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGGCGAACTGGCTAAACGTGGCGAAATCGATAAGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAA  >  NC_000913/125578‑125868
                                                                                                                                                  |                                                                                                                                                
gggCGAACTGGCTAAACGTGGCGAAATCGATAAGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCgg                                                                                                                                                <  1:622116/149‑1 (MQ=255)
  gCGAACTGGCTAAACGTGGCGAAATCGATAAGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCgggg                                                                                                                                              <  1:191604/149‑1 (MQ=255)
           cTAAACGTGGCGAAATCGATAAGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAg                                                                                                                                     <  1:11645/149‑1 (MQ=255)
            tAAACGTGGCGAAATCGATAAGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAg                                                                                                                                     >  3:470637/1‑148 (MQ=255)
            tAAACGTGGCGAAATCGATAAGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAg                                                                                                                                     <  4:470637/148‑1 (MQ=255)
             aaaCGTGGCGAAATCGATAAGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGtt                                                                                                                                   <  2:438741/149‑1 (MQ=255)
               aCGTGGCGAAATCGATAAGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGa                                                                                                                                 <  3:593908/149‑1 (MQ=255)
                     cGAAATCGATAAGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCTAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCAc                                                                                                                           >  2:63201/1‑149 (MQ=255)
                     cGAAATCGATAAGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCAc                                                                                                                           >  1:190367/1‑149 (MQ=255)
                     cGAAATCGATAAGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCAc                                                                                                                           <  2:359611/149‑1 (MQ=255)
                      gAAATCGATAAGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCAcc                                                                                                                          >  3:484523/1‑149 (MQ=255)
                      gAAATCGATAAGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCAcc                                                                                                                          <  2:307630/149‑1 (MQ=255)
                        aaTCGATAAGAAAGTGGTTGCTGACGCAATCACCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCga                                                                                                                        <  1:179538/149‑1 (MQ=255)
                             aTAAGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATcc                                                                                                                   <  1:186649/149‑1 (MQ=255)
                               aaGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTg                                                                                                                 <  3:120973/149‑1 (MQ=255)
                                aGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTgg                                                                                                                <  1:525845/149‑1 (MQ=255)
                                aGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTgg                                                                                                                <  3:357012/149‑1 (MQ=255)
                                     gttgTTGCTCACTCAATCGCCCATATCAACATCGATGCAGATAAAGCTAACCCGTGGCTGGCGCAAGAGGTAAAAGAATCATGGCTATCGAAATCAAAGTACCGGATATCGGGGCTGATGAAATTGAAATCACCTAGATCCTGGTCaaa                                                                                                           <  4:301082/146‑1 (MQ=255)
                                        gTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTg                                                                                                        <  4:509407/149‑1 (MQ=255)
                                        gTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTg                                                                                                        >  2:128492/1‑149 (MQ=255)
                                             tGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGa                                                                                                   >  1:401552/1‑149 (MQ=255)
                                             tGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGa                                                                                                   <  3:43925/149‑1 (MQ=255)
                                              gACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGAc                                                                                                  >  4:416469/1‑149 (MQ=255)
                                                 gCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACaaa                                                                                               <  3:424231/149‑1 (MQ=255)
                                                    aTCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGtt                                                                                            >  4:185589/1‑149 (MQ=255)
                                                     tCGCCAAATTCACCATCGCTGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTg                                                                                           >  4:376246/1‑149 (MQ=255)
                                                         cAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGc                                                                                                     >  1:215905/1‑135 (MQ=255)
                                                         cAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGc                                                                                                     <  2:215905/135‑1 (MQ=255)
                                                               cAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGaa                                                                                         <  4:129178/141‑1 (MQ=255)
                                                               cAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGaa                                                                                         >  3:129178/1‑141 (MQ=255)
                                                                aaCATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAg                                                                                >  3:577980/1‑149 (MQ=255)
                                                                aaCATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAACCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAg                                                                                >  4:503165/1‑149 (MQ=255)
                                                                   aTCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCg                                                                             >  4:534470/1‑149 (MQ=255)
                                                                   aTCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCg                                                                             >  2:201379/1‑149 (MQ=255)
                                                                     cGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCt                                                                           <  3:416469/149‑1 (MQ=255)
                                                                      gATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATTGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTg                                                                          >  2:390509/1‑149 (MQ=255)
                                                                      gATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATTGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTg                                                                          >  2:14325/1‑149 (MQ=255)
                                                                      gATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTg                                                                          >  2:192834/1‑149 (MQ=255)
                                                                       aTGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCtaa                                                                         <  1:295396/149‑3 (MQ=255)
                                                                       aTGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGa                                                                         <  4:484523/149‑1 (MQ=255)
                                                                            gATAAAGTTAACTCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGACGCCGAACAGTCGCTGATCAcc                                                                    >  2:203279/1‑149 (MQ=255)
                                                                                aaGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAg                                                                >  4:395398/1‑149 (MQ=255)
                                                                                 aGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATcc                                                                                                                   <  1:428985/97‑1 (MQ=255)
                                                                                 aGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATcc                                                                                                                   >  2:428985/1‑97 (MQ=255)
                                                                                  gTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGaa                                                              >  3:242282/1‑149 (MQ=255)
                                                                                             tCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGc                                                   <  2:213966/149‑1 (MQ=255)
                                                                                             tCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGc                                                   <  3:376694/149‑1 (MQ=255)
                                                                                               tGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCct                                                 <  4:577980/149‑1 (MQ=255)
                                                                                                ggCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCctc                                                <  3:503165/149‑1 (MQ=255)
                                                                                                ggCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCctc                                                <  4:242282/149‑1 (MQ=255)
                                                                                                  cGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTa                                              <  1:63201/149‑1 (MQ=255)
                                                                                                   gTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATTGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTAt                                             <  1:14325/149‑1 (MQ=255)
                                                                                                   gTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATTGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTAt                                             <  1:390509/149‑1 (MQ=255)
                                                                                                           tAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTc                                     >  2:612780/1‑149 (MQ=255)
                                                                                                            aaaaGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTcc                                    >  2:365469/1‑149 (MQ=255)
                                                                                                            aaaaGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTcc                                    >  3:136488/1‑149 (MQ=255)
                                                                                                                 aataatGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGtctc                               <  4:136488/149‑1 (MQ=255)
                                                                                                                         cTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCg                       >  4:193901/1‑149 (MQ=255)
                                                                                                                                 aTCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGt               >  4:416468/1‑149 (MQ=255)
                                                                                                                                     aaGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTaaa           <  4:222194/149‑1 (MQ=255)
                                                                                                                                      aGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAag          <  2:401552/149‑1 (MQ=255)
                                                                                                                                       gTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAaga         <  3:185589/149‑1 (MQ=255)
                                                                                                                                         aCCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAgaga       >  3:189048/1‑149 (MQ=255)
                                                                                                                                         aCCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAgaga       >  3:281795/1‑149 (MQ=255)
                                                                                                                                         aCCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAgaga       <  3:534470/149‑1 (MQ=255)
                                                                                                                                          ccGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGAt      <  1:561121/149‑1 (MQ=255)
                                                                                                                                           cGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATc     >  1:597500/1‑149 (MQ=255)
                                                                                                                                              aCATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCaaa  <  2:190367/149‑1 (MQ=255)
                                                                                                                                              aCATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCaaa  <  1:612780/149‑1 (MQ=255)
                                                                                                                                              aCATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCaaa  <  3:416468/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGGCGAACTGGCTAAACGTGGCGAAATCGATAAGAAAGTGGTTGCTGACGCAATCGCCAAATTCAACATCGATGCAGATAAAGTTAACCCGCGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAA  >  NC_000913/125578‑125868

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