Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 125,817 A→G 7.6% K41K (AAA→AAG)  aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913125,8170AG 7.6% 238.5 / 5.9 90K41K (AAA→AAG) aceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base A (39/43);  new base G (4/3);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGA  >  NC_000913/125670‑125965
                                                                                                                                                   |                                                                                                                                                    
gTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAg                                                                                                                                                     <  2:625289/149‑1 (MQ=255)
 tCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGc                                                                                                                                                    >  2:380057/1‑149 (MQ=255)
 tCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTAGCTGATCACAGTAGAAGGCGACAAAGc                                                                                                                                                    >  1:907266/1‑149 (MQ=255)
   tGGCGTTAGAGGTAAAAGAATAATGGCGATCGAAATCAAAGTACCGGACCTCGGGGCGGATGTAGGTGAAATACCCGACATCCTGGTCAAAGAGGGCGACAAAGTTGAAGCCCTACAGTCGCTGATCACCGTAGCAGGCGACAAAGccg                                                                                                                                                  >  2:477824/1‑148 (MQ=255)
   tGGCGTAAGAGGTAAACGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAACTCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCct                                                                                                                                                  <  1:354392/149‑1 (MQ=255)
    ggCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGGTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCTACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCctc                                                                                                                                                 <  2:1487920/149‑1 (MQ=255)
    ggCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCctc                                                                                                                                                 >  2:1536897/1‑149 (MQ=255)
       gTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCCAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTAt                                                                                                                                              >  2:1221465/1‑149 (MQ=255)
           gagGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGaa                                                                                                                                          <  2:605713/149‑1 (MQ=255)
           gagGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGaa                                                                                                                                          <  2:1097238/149‑1 (MQ=255)
           gagGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGaa                                                                                                                                          <  1:225401/149‑1 (MQ=255)
             ggTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGt                                                                                                                                        >  1:191766/1‑149 (MQ=255)
               tAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCTCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTc                                                                                                                                      >  2:60455/1‑149 (MQ=255)
               tAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTc                                                                                                                                      >  1:216106/1‑149 (MQ=255)
               tAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTc                                                                                                                                      >  2:267390/1‑149 (MQ=255)
                aaaaGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTcc                                                                                                                                     >  2:25355/1‑149 (MQ=255)
                 aaaGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGCCAAAGTTGAAGCCGAACAGGCGCTGATCACCGTAGAAGGCGACAACGCCTCTATGGAAGTACCg                                                                                                                                    >  2:876852/1‑149 (MQ=255)
                  aaGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGt                                                                                                                                   <  2:1433709/149‑1 (MQ=255)
                       taatGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCg                                                                                                                              <  1:242507/149‑1 (MQ=255)
                        aagggCTATCGAACTCAAAGTACAGGACATCGGGGCTGATTATCTTGAAATCACCGAGATCCAGGTCAAAGTGCGCGACAAAGTTTAACACGAACAGTCGCTGATCACCGTAGACGGCGACAAAGCCGCTATGTAAGTTCCGTCTCCGc                                                                                                                             <  2:1352184/146‑1 (MQ=25)
                              tATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCgg                                                                                                                       <  2:1117436/149‑1 (MQ=255)
                              cagcGAAATCAAAGTACCGGACCTCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCgg                                                                                                                       <  2:227365/146‑1 (MQ=255)
                                  gAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTAt                                                                                                                   >  2:187399/1‑149 (MQ=255)
                                   aaaTCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCCGGCGGGTATc                                                                                                                  >  2:295016/1‑149 (MQ=255)
                                       cAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTa                                                                                                              >  2:1524100/1‑149 (MQ=255)
                                        aaaGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTaa                                                                                                             >  2:193384/1‑149 (MQ=255)
                                         aaGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTaat                                                                                                            <  1:267390/149‑2 (MQ=255)
                                         aaGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTaaa                                                                                                            <  2:144564/149‑1 (MQ=255)
                                         aaGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGATGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTAAGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTaaa                                                                                                            >  1:978555/1‑149 (MQ=255)
                                             aCCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGATCAGAAGGCTGGGATCGTTAAagata                                                                                                        >  2:1329344/1‑147 (MQ=21)
                                             aCCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAgaga                                                                                                        >  2:87253/1‑149 (MQ=255)
                                              ccGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGCCAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGAt                                                                                                       >  2:855585/1‑149 (MQ=255)
                                               cGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCg                                                                                                                 >  1:357231/1‑138 (MQ=255)
                                               cGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCg                                                                                                                 <  2:357231/138‑1 (MQ=255)
                                                  aCATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCaaa                                                                                                   <  1:564448/149‑1 (MQ=255)
                                                  aCATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCaaa                                                                                                   <  2:464431/149‑1 (MQ=255)
                                                          gCTTATGAAGGTGAACTCACCGAGCTCCTGGTCAAAGTGGGAGACAAAGTTGAAGCCGAACAGGCGCTGATCACCGTAGAAGGCGACAAAGCCGCTAGGGAAGTGCCGGCTCCGCAGGCGGGTATCGTTAACTAGATCAAAGTCTCTGt                                                                                           <  2:1122928/149‑1 (MQ=21)
                                                          gCTGATGAAGTTGAACTCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGt                                                                                           <  2:1310800/149‑1 (MQ=255)
                                                            tgatgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTg                                                                                         >  2:1220149/1‑149 (MQ=255)
                                                             gatgaAGTTGAAATCACCGAGATCCTGGGCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGCTGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTgg                                                                                        <  1:1121487/149‑1 (MQ=255)
                                                              atgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCAGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGa                                                                                                                                           >  2:1107055/1‑97 (MQ=255)
                                                              atgaAGTTGAAATCACCGACATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGa                                                                                                                                           <  1:1107055/97‑1 (MQ=255)
                                                                  aggtGAAATCACCGAGATCCTGGTCAAAGTGGTCGACAATGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCTACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATa                                                                                   <  2:1046514/146‑1 (MQ=255)
                                                                    ttGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATaaa                                                                                 >  2:24916/1‑149 (MQ=255)
                                                                      gAAATCCCCGAGATCCTGGTCAAAGTGGGCGACAAAGGTGAAGCGGCTCAGGCGCTGATCACCGTAGCAGGCGACAAGGCGTCGATGGAAGTTCCGGCTCGGGTGGCGGGGAGCGGGAAAGAGAGCCAAGTGGCGGTTGGTGACaaggg                                                                               >  2:719983/1‑146 (MQ=255)
                                                                      gAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGCTCAAAGTCTCTGTTGGCGATAAAAc                                                                               >  2:664976/1‑149 (MQ=255)
                                                                      gAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAAc                                                                               >  1:542538/1‑149 (MQ=255)
                                                                       aaaTCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAAAATTCGCTGATCACCGTAGAAGGCGACTAAGGCTCTATGGAAGATCCTTCTCCGCAGGAGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGGTatcacc                                                                              >  1:1400643/1‑144 (MQ=25)
                                                                        aaTCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAAccc                                                                             <  1:87253/149‑1 (MQ=255)
                                                                           cACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCGGATCACCGTAGAAGGCGACTAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGTGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGa                                                                          <  2:1213685/149‑1 (MQ=255)
                                                                             ccGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGCGCTCAACGTCTCGGGTGGCGAGAAACCCCCGAcc                                                                        >  1:1571568/1‑149 (MQ=21)
                                                                              cGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGaa                                                                                                                                          <  2:732118/82‑1 (MQ=255)
                                                                              cGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGaa                                                                                                                                          >  1:732118/1‑82 (MQ=255)
                                                                                agttCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCGGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGgc                                                                     <  1:1207038/146‑1 (MQ=255)
                                                                                agaTCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGgc                                                                     >  2:48378/1‑149 (MQ=255)
                                                                                 gaTCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCTTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGTTCAAAGTCTCTGTTGGCGATAAAACCCAGACCGgcg                                                                    <  1:90967/149‑1 (MQ=255)
                                                                                  aTCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGAGCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGgcgc                                                                   <  2:1057413/149‑1 (MQ=255)
                                                                                    ccTGGGCACAGTGGGCGACAAAGTTGAAGCCGAACAGGCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCAc                                                                 <  2:873290/149‑1 (MQ=255)
                                                                                    ccTCGTCAAAGTGGGCGACAAAGTTGAACCCGAACAGTCGCTGATCACCGTAGAAGTCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCAc                                                                 <  2:1276498/149‑1 (MQ=255)
                                                                                     cTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGCCCGGCGCACt                                                                >  2:739803/1‑149 (MQ=255)
                                                                                     cTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACt                                                                >  2:138381/1‑149 (MQ=255)
                                                                                             aGTGGGCGACAAAGTTTAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGa                                                        <  1:477824/149‑1 (MQ=255)
                                                                                             aGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCgg                                                                                                                                   <  1:73215/74‑2 (MQ=255)
                                                                                             aGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCgg                                                                                                                                   >  2:73215/1‑73 (MQ=255)
                                                                                                  gCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTc                                                   <  2:216106/149‑1 (MQ=255)
                                                                                                    gACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACGGGCGCACTGATGATGATTTTCGa                                                 >  1:542420/1‑149 (MQ=255)
                                                                                                          gTTGAAGCCGAACTGTCGCTGATCTCCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTATAGAGATCATATTCTCTGTTGTGGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTccgc                                           <  2:1365313/149‑1 (MQ=255)
                                                                                                          gTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAATGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTccgc                                           >  2:1154549/1‑149 (MQ=255)
                                                                                                           ttGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTccgcc                                          <  1:1221465/149‑1 (MQ=255)
                                                                                                                  cGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTg                                   <  1:734506/149‑1 (MQ=255)
                                                                                                                    aaCAGTCGCTGATCAACGTAGAACGCTACAAACCCTCTATGGAACTACCGTCTCCGCAGGCGGGTATCGTTAAAGAGATGAAAGTCTCTGATGGCGATAAAACCCAGACCGTCGCAATGATGATGATTTTCGATTCCGCCGACGGTgca                                 >  1:300286/1‑149 (MQ=255)
                                                                                                                       aGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGAAAAAACCCAGACCCGCGCACTGATTATGATTTTAGATTCCTCCGACGGTgcagca                              >  2:1018705/1‑149 (MQ=255)
                                                                                                                           gCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACg                          >  1:467665/1‑149 (MQ=255)
                                                                                                                           gCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACg                          <  2:191766/149‑1 (MQ=255)
                                                                                                                           gCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACg                          >  2:1174765/1‑149 (MQ=255)
                                                                                                                           gCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACg                          >  2:561632/1‑149 (MQ=255)
                                                                                                                            cTGATCACCGTAGAAGTCGAAAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGTTCAAAGTCTCTGTTGTCGATACAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGc                         <  2:907266/149‑1 (MQ=255)
                                                                                                                               aTCACCGGAGAAGGCGACAAAGCCTCTATGGCAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGc                      <  2:429634/149‑1 (MQ=255)
                                                                                                                                tCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGAATTTCGATTCCGCCGACGGTGCAGCAGACGCTGCa                     <  1:380057/149‑1 (MQ=255)
                                                                                                                                  aCCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCACTGTTGGCGATAAAACCCAGACCCGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCAcc                   >  1:1303324/1‑149 (MQ=255)
                                                                                                                                    cGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTg                 >  1:101512/1‑149 (MQ=255)
                                                                                                                                     gTAGAAGGCTAAAAAGCCTCTATGCAAGTTCCGTCTCCGCAGTCGTTTATCGTTAAAGAGATCAAAGTCTCTTTTGGCGATAAAACCCAGACCGGCGCACTGATTATGAATTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGc                <  1:876852/149‑1 (MQ=255)
                                                                                                                                       aGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGCCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTc              <  1:295016/149‑1 (MQ=255)
                                                                                                                                         aaGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAg            <  1:1495489/149‑1 (MQ=255)
                                                                                                                                         aaGGCGACAAAGCCGCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTCAAGCGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTACGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAg            <  1:129982/149‑1 (MQ=255)
                                                                                                                                          aGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAgg           >  2:1213011/1‑149 (MQ=255)
                                                                                                                                           ggCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGc          <  1:932896/149‑1 (MQ=255)
                                                                                                                                              gACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATAATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCaga       <  1:138381/149‑1 (MQ=255)
                                                                                                                                              gACAAAGCCTCTATGGAAGATCCGTCTCCGCATGCGGGAATCGTTTAAGAGAACAAAGTCTCTCTTGCCGATAACACCCAGAACGGCGCACTTATTATGATTTTCGATTCCGCCGACCGTGCAGCAGACGCTGCCCCTGCTGAGGCaga       >  1:1518707/1‑149 (MQ=255)
                                                                                                                                                cAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCagaag     <  2:542538/149‑1 (MQ=255)
                                                                                                                                                 aaaGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCagaaga    >  2:583996/1‑149 (MQ=255)
                                                                                                                                                   aGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCagaagaga  <  1:48378/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGA  >  NC_000913/125670‑125965

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