Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 125,817 A→G 6.4% 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 6.4% 208.7 / 2.4 76K41K (AAA→AAGaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base A (31/39);  new base G (2/3);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAG  >  NC_000913/125669‑125964
                                                                                                                                                    |                                                                                                                                                   
cgTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCTAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAg                                                                                                                                                    <  2:1286079/150‑1 (MQ=255)
cgTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAGGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAg                                                                                                                                                    >  1:1167779/1‑150 (MQ=255)
cgTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAg                                                                                                                                                    <  1:797064/150‑1 (MQ=255)
cgTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAg                                                                                                                                                    <  1:378044/150‑1 (MQ=255)
    tGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCctc                                                                                                                                                <  2:455493/150‑1 (MQ=255)
    tGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCctc                                                                                                                                                <  2:598993/150‑1 (MQ=255)
    tGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCctc                                                                                                                                                <  2:650727/150‑1 (MQ=255)
      gCGTAAGAGGTAAAAGAATAATGGCTATCGAACTCAAATTCCCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTa                                                                                                                                              <  2:898422/150‑1 (MQ=255)
         tAAGAGGTAAAAGATTTATGGCTATAGAAATCAAAGTTCCGGACATCGGGTGTGATGAAGTTGAAATCACCGATATCCTGGTGAAAGTGGGCGACAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATgg                                                                                                                                           <  2:222406/150‑1 (MQ=17)
         tAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATgg                                                                                                                                           >  1:509884/1‑150 (MQ=255)
         tAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCAGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATgg                                                                                                                                           >  2:738777/1‑150 (MQ=255)
          acgagGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGATATCATGGTCAAAGTGGGCGACAAAGTTGCAGCCGAACAGTCTCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGa                                                                                                                                          <  1:738777/148‑1 (MQ=25)
          aaGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGa                                                                                                                                          <  2:404885/150‑1 (MQ=255)
          aaGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGa                                                                                                                                          <  2:1017749/150‑1 (MQ=255)
            gagGTAATAGAATAATGGCTATCGAAATCAAAGTACAGTACCTCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGTAGACAAAGTTGAATCCGAACAGTCGCTGATCACCGTATAAGCCGACAACGACTCTATGGatg                                                                                                                                        >  1:821559/1‑148 (MQ=25)
             agGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGt                                                                                                                                       <  2:789778/150‑1 (MQ=255)
               gTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTc                                                                                                                                     >  1:1185532/1‑150 (MQ=255)
                tAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTcc                                                                                                                                    >  1:1277226/1‑150 (MQ=255)
                tAAAAGAATAATGGCTATCGAAATAAAAGTACCGGACATCGGGGCTGATGAATTTGAAATCACATAGATAATGGTCAAAGTGGGAGACAAAGTTGAAGAAGAACAGACGCTGATAAACGTAGAAGGCGACAAAGCCTCTATGTATGTTca                                                                                                                                    >  1:1315411/1‑149 (MQ=255)
                 aaaaTAATAATGGCTATCGAAATCAAAGTACCGGACCTCGGGGATGATGAAGTTGAAATCACCGAGATACTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCg                                                                                                                                   >  2:827606/1‑150 (MQ=25)
                 aaaaGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCg                                                                                                                                   <  1:723823/150‑1 (MQ=255)
                     gAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGATGATGAATTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGtctc                                                                                                                               <  2:188497/150‑1 (MQ=255)
                       ataatGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGTCTCCg                                                                                                                             <  2:385882/150‑1 (MQ=255)
                         aatGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCa                                                                                                                           <  2:1153567/150‑1 (MQ=255)
                             gCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAATGCGACAAAGCCTCTATGGAAGATCCGTCTAATCAGGCg                                                                                                                       >  2:442945/1‑150 (MQ=25)
                             gCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCgtttgct                                                                                                                       <  2:509884/150‑7 (MQ=18)
                              cTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCTCCgtttgctg                                                                                                                      >  1:451666/1‑143 (MQ=255)
                               tATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGACGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCggg                                                                                                                     <  2:80503/150‑1 (MQ=255)
                                     aaTCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCCCTGATCACCGAAGCATGCGACAAAGGCTATATGGAAGTTCCGTCGCCGCAGGCGCGTATCGt                                                                                                               >  2:614293/1‑150 (MQ=255)
                                           aGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAaga                                                                                                         <  2:543677/150‑1 (MQ=255)
                                            gTACCGGACATCGGGGCGGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTaatgtg                                                                                                        <  2:497219/150‑5 (MQ=255)
                                              aCCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATc                                                                                                     >  2:1145820/1‑150 (MQ=255)
                                               ccGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGAta                                                                                                     >  2:204701/1‑149 (MQ=255)
                                               ccGGACAGCGGGGCTGAAGAAGTTGAAATCACCGAGATACTTTTCAAAGTTGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATc                                                                                                     <  1:1118782/150‑1 (MQ=25)
                                                cGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCa                                                                                                    <  2:56532/150‑1 (MQ=255)
                                                   aCATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAg                                                                                                 <  2:451666/150‑1 (MQ=11)
                                                     aTCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGAGCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGtc                                                                                               <  2:535391/150‑1 (MQ=255)
                                                           gCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGTATCCGGATGCGGGCAGCGTGAAAGAGATCAAAGtgaacgtg                                                                                         >  2:641794/1‑143 (MQ=14)
                                                           gCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGtt                                                                                         >  1:1013150/1‑150 (MQ=255)
                                                            cTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCATCTATGGAAGTTCCGTCTCCGAATGCGGGAAACGTGAAAGAGATCAAAGtctaagtgg                                                                                        >  2:1057111/1‑144 (MQ=17)
                                                              tatgaAGTTGAAATCATCGAGATCCTGGTCAAAGTGGGCGACAAAGTTTAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCATCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGc                                                                                      <  2:500212/149‑1 (MQ=255)
                                                              tatgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTTGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGc                                                                                      <  2:20600/149‑1 (MQ=255)
                                                              gatgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGc                                                                                      <  2:569111/150‑1 (MQ=255)
                                                                 taaGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGTACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGAt                                                                                   <  2:1315411/149‑1 (MQ=255)
                                                                  aaGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGGAGGCTGGAATCGTTAAAGAGATCAAAGTCTCGGTTGgagaaa                                                                                  >  1:722281/1‑148 (MQ=18)
                                                                    gTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCATCTATGGAAGTTCCGTCTCCGCATGCTGGAACCGTGAAAGAGATCAAAGTCAAAGTTGGAGAaaaa                                                                                >  2:1270289/1‑150 (MQ=14)
                                                                     ttGAAATCACCGAGATCCTGGTCAAAGTGCGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGGTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATaaaa                                                                               <  1:685850/150‑1 (MQ=255)
                                                                       gAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAAcc                                                                             >  1:27719/1‑150 (MQ=255)
                                                                         aaTCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCa                                                                           <  1:1225569/150‑1 (MQ=255)
                                                                          gtaaCCGAGATCCTGGGCAAAGTGGGCGACAAAGGTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAg                                                                          <  2:249397/147‑1 (MQ=255)
                                                                          aTCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAg                                                                          >  2:422322/1‑150 (MQ=255)
                                                                              ccGAGATCCTGGTCAAAGTGGTCGACAAAGTTGAAGCCGAACAGTCGCTGAGCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGGTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCg                                                                      <  1:188258/150‑1 (MQ=255)
                                                                              ccGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTGTGTTGGCGCTTAAACGCAGACCg                                                                      >  2:112135/1‑150 (MQ=255)
                                                                              ccGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGCCCg                                                                      >  1:1141580/1‑150 (MQ=255)
                                                                              ccGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCg                                                                      >  1:2857/1‑150 (MQ=255)
                                                                                gagaTCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGgc                                                                    >  2:372701/1‑150 (MQ=255)
                                                                                    tCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAGCCAAGACCGgcgcgc                                                                >  1:946829/1‑148 (MQ=255)
                                                                                     ccTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACt                                                               >  2:1283639/1‑150 (MQ=255)
                                                                                     ccTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACt                                                               <  1:475703/150‑1 (MQ=255)
                                                                                            aaaGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATg                                                        <  1:690183/150‑1 (MQ=255)
                                                                                                 gggCGACAAAGTTGAAGCCGAGCAGTCGCTGATCACCGTAGAAGGTGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGAtttt                                                   >  1:164358/1‑150 (MQ=255)
                                                                                                 gggCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGAAGGCGGGAAACGTTAAAGAGATCAAAGTAGCTCTTGGCGAACACACCCAGAGCTTCGCACTAAATATGAttat                                                   >  2:434291/1‑148 (MQ=18)
                                                                                                  ggCGACAAATTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTc                                                  <  2:266857/150‑1 (MQ=255)
                                                                                                  ggCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTTGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTc                                                  <  2:128904/150‑1 (MQ=255)
                                                                                                  ggCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTc                                                  >  2:1254634/1‑150 (MQ=255)
                                                                                                     gACAAAGTTGAACACGAACAGTCGCTGATCACCGTAGCAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGCATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGAt                                               <  2:952947/150‑1 (MQ=255)
                                                                                                       aaaaGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTc                                             <  2:444551/149‑1 (MQ=255)
                                                                                                           gTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTccgcc                                         <  2:166645/150‑1 (MQ=255)
                                                                                                           gTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGAGATAAAAACCAGACCGGCGCCCTGATTATGATATTCGATTacgcc                                         >  1:559858/1‑150 (MQ=37)
                                                                                                                      aCAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTgcagc                              <  1:594095/150‑1 (MQ=255)
                                                                                                                       ccggcgcgGATCACCGTAGAAGGCGACAAAGCCGCTATGGTAGTTCCGTCTCCGCAGGCGGGTATAGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATCATGATTTTCGATTCCGCCGACGGTgcagca                             <  2:800315/146‑1 (MQ=255)
                                                                                                                       caggcccTGATCACCGTAGAAGGCCACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGTCGGGAATCGATAAAGCGATCAAAGTCTCTGTTGCCGATAAAACCCAGACCGGCGCACTGATGATGATTTTCGATTCCGCCGACGGTgcagca                             <  1:184069/144‑1 (MQ=255)
                                                                                                                                            ggCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCATGCGGGTATCGTTAAAGAGATCAAAGTCTATGTTGGCGATAAGACCAAGACAGACGCACTGATATTGATTTTCGATTCCGCCTACGGTGAAGCAGACGCTGCACCTGCGCAGGct        >  1:391730/1‑149 (MQ=255)
                                                                                                                                             gCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGGACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCag       >  2:1234664/1‑150 (MQ=255)
                                                                                                                                             gCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCag       >  1:878167/1‑150 (MQ=255)
                                                                                                                                                 cAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCACTGTTGGCGATAAAACCCAGACCCGCGGACTGATTATGAATTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGCCagaag  >  2:460989/1‑149 (MQ=255)
                                                                                                                                                  aaaGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCagaagag  >  1:889750/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAG  >  NC_000913/125669‑125964

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: