Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 125,820 C→T 8.1% A42A (GCC→GCT aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913125,8200CT 8.1% 228.3 / 8.2 88A42A (GCC→GCTaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base C (37/43);  new base T (4/3);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAG  >  NC_000913/125673‑125967
                                                                                                                                                   |                                                                                                                                                   
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:1097238/149‑1 (MQ=255)
        gagGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGaa                                                                                                                                            <  2:605713/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                                                                                                                                        >  2:267390/1‑149 (MQ=255)
            tAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTc                                                                                                                                        >  1:216106/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                                                                                                                   <  2:357231/138‑1 (MQ=255)
                                            cGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCg                                                                                                                   >  1:357231/1‑138 (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                            >  2:561632/1‑149 (MQ=255)
                                                                                                                        gCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACg                            >  2:1174765/1‑149 (MQ=255)
                                                                                                                        gCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACg                            <  2:191766/149‑1 (MQ=255)
                                                                                                                        gCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACg                            >  1:467665/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)
                                                                                                                                                 gCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAACGTCTCTGTTGGCTCTAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCagaagagaa   <  2:488673/149‑1 (MQ=255)
                                                                                                                                                 gCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCAc                      >  1:1259306/1‑130 (MQ=255)
                                                                                                                                                 gCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCAc                      <  2:1259306/130‑1 (MQ=255)
                                                                                                                                                  ccTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGGCGCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCagaagagaag  <  1:187399/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAG  >  NC_000913/125673‑125967

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: