Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 753,205 C→A 6.1% Q21K (CAG→AAG)  sdhB → succinate dehydrogenase, FeS subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273753,2050CA 6.1% 1077.7 / 31.8 277Q21K (CAG→AAG) sdhBsuccinate dehydrogenase, FeS subunit
Reads supporting (aligned to +/- strand):  ref base C (119/141);  new base A (13/4);  total (132/145)
Fisher's exact test for biased strand distribution p-value = 2.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 282 total aligned reads displayed.

CCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAG  >  CP009273/753056‑753354
                                                                                                                                                     |                                                                                                                                                     
ccATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGc                                                                                                                                                       <  2:2512873/150‑1 (MQ=255)
 cATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCAGTCCCGCCGAAGATTCGTACTGACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCa                                                                                                                                                      <  2:3278913/150‑1 (MQ=255)
     aCGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGAt                                                                                                                                                  <  1:762571/150‑1 (MQ=255)
      cgcgCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGAtt                                                                                                                                                 <  2:3038886/150‑1 (MQ=255)
           cGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACAcc                                                                                                                                            >  2:2148754/1‑150 (MQ=255)
            gAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACAccc                                                                                                                                           <  2:1018905/150‑1 (MQ=255)
                  tCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAg                                                                                                                                     <  1:574862/150‑1 (MQ=255)
                   cAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGc                                                                                                                                    <  1:2429311/150‑1 (MQ=255)
                   cAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGc                                                                                                                                    >  2:2893202/1‑150 (MQ=255)
                     aCATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCgg                                                                                                                                  <  1:1435123/150‑1 (MQ=255)
                       aTGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGAt                                                                                                                                >  2:3583030/1‑150 (MQ=255)
                       aTGGAACCGAAACTGCGCCCGGCAGTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGAt                                                                                                                                <  1:3262791/150‑1 (MQ=255)
                          gAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGaa                                                                                                                             >  1:42389/1‑150 (MQ=255)
                             ccGAAACTGCGCCCGGCAGTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGt                                                                                                                          <  1:4148986/150‑1 (MQ=255)
                              cGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTc                                                                                                                         <  1:71812/150‑1 (MQ=255)
                              cGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTc                                                                                                                         <  2:4237616/150‑1 (MQ=255)
                              cGAAACTGCGCCCGGCAGTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTc                                                                                                                         <  1:2737399/150‑1 (MQ=255)
                               gAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTcg                                                                                                                        >  2:70169/1‑150 (MQ=255)
                                 aaCTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTcgcg                                                                                                                      >  1:1925409/1‑150 (MQ=255)
                                  aCTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGa                                                                                                                     >  2:1989/1‑150 (MQ=255)
                                       gcCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGAAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACatga                                                                                                                <  2:2595067/150‑1 (MQ=255)
                                        cccGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACatgat                                                                                                               <  2:363520/150‑1 (MQ=255)
                                        cccGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACatgat                                                                                                               >  2:3378319/1‑150 (MQ=255)
                                         ccGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACatgatg                                                                                                              >  1:1556757/1‑150 (MQ=255)
                                             cATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGAtgctgc                                                                                                          >  1:3941388/1‑150 (MQ=255)
                                             cATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGAtgctgc                                                                                                          >  1:3941418/1‑150 (MQ=255)
                                                tCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTgg                                                                                                       <  2:1695493/150‑1 (MQ=255)
                                                  ccgccgAAGATTCGTACTTACTACTGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGAt                                                                                                     <  1:2953721/150‑1 (MQ=255)
                                                  ccgccgAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGAt                                                                                                     <  2:3202320/150‑1 (MQ=255)
                                                   cgccgAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATg                                                                                                    >  1:537558/1‑150 (MQ=255)
                                                   cgccgAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATg                                                                                                    <  1:3622500/150‑1 (MQ=255)
                                                    gccgAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATgc                                                                                                   <  2:3774852/150‑1 (MQ=255)
                                                        aaGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCtt                                                                                               >  2:1459608/1‑150 (MQ=255)
                                                        aaGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCtt                                                                                               <  1:2199638/150‑1 (MQ=255)
                                                        aaGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCtt                                                                                               <  2:3424293/150‑1 (MQ=255)
                                                           aTTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATc                                                                                            >  2:547041/1‑150 (MQ=255)
                                                            ttCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATcc                                                                                           <  2:3799165/150‑1 (MQ=255)
                                                             tCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCa                                                                                          <  1:3128795/150‑1 (MQ=255)
                                                              cGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAg                                                                                         <  1:440537/150‑1 (MQ=255)
                                                              cGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAg                                                                                         <  2:1900828/150‑1 (MQ=255)
                                                              cGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAg                                                                                         <  1:2872290/150‑1 (MQ=255)
                                                              cGTACTGACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAg                                                                                         <  2:2469734/150‑1 (MQ=255)
                                                                 acttactAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTa                                                                                      <  2:3951249/150‑1 (MQ=255)
                                                                 acttactAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTa                                                                                      >  2:4067806/1‑150 (MQ=255)
                                                                 acttactAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTa                                                                                      <  1:1187746/150‑1 (MQ=255)
                                                                 acttactAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTa                                                                                      >  1:2317973/1‑150 (MQ=255)
                                                                 acttactAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTa                                                                                      <  1:170492/150‑1 (MQ=255)
                                                                   ttactAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGAAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaaa                                                                                    >  1:3868068/1‑150 (MQ=255)
                                                                   ttactAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGAAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaaa                                                                                    >  2:2397271/1‑150 (MQ=255)
                                                                    tactAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaaaa                                                                                   <  2:3506435/150‑1 (MQ=255)
                                                                     actAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAag                                                                                  <  2:1094875/150‑1 (MQ=255)
                                                                       tAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAagag                                                                                >  1:1954871/1‑150 (MQ=255)
                                                                        aaTGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAgaga                                                                               <  1:1767639/150‑1 (MQ=255)
                                                                          tGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGaaa                                                                             <  1:4015315/150‑1 (MQ=255)
                                                                          tGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGaaa                                                                             <  1:2833534/150‑1 (MQ=255)
                                                                          tGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGaaa                                                                             <  2:1323215/150‑1 (MQ=255)
                                                                            cGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGa                                                                           >  1:2048782/1‑150 (MQ=255)
                                                                             ggAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGAt                                                                          <  1:2578743/150‑1 (MQ=255)
                                                                             ggAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGAt                                                                          <  1:405744/150‑1 (MQ=255)
                                                                                       aaaTGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGACGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGt                                                                >  2:1780344/1‑150 (MQ=255)
                                                                                       aaaTGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGt                                                                >  2:1547316/1‑150 (MQ=255)
                                                                                       aaaTGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGt                                                                >  2:158393/1‑150 (MQ=255)
                                                                                           gagaCTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGtt                                                            >  2:2691391/1‑150 (MQ=255)
                                                                                            agaCTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTc                                                           >  2:2988187/1‑150 (MQ=255)
                                                                                            agaCTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTc                                                           >  1:4306802/1‑150 (MQ=255)
                                                                                             gaCTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTcc                                                          >  2:1162206/1‑150 (MQ=255)
                                                                                              aCTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTccg                                                         <  2:596571/150‑1 (MQ=255)
                                                                                                  gAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCcgccgc                                                     >  1:548878/1‑150 (MQ=255)
                                                                                                   aGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCt                                                    >  1:1583472/1‑150 (MQ=255)
                                                                                                   aGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCt                                                    >  1:3771082/1‑150 (MQ=255)
                                                                                                    gTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGACGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTc                                                   >  2:1201513/1‑150 (MQ=255)
                                                                                                        tCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGc                                               <  1:2804484/150‑1 (MQ=255)
                                                                                                          aaTTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCg                                             <  1:1081862/150‑1 (MQ=255)
                                                                                                              tATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGaa                                         <  1:2778339/150‑1 (MQ=255)
                                                                                                              tATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGaa                                         <  1:4194341/150‑1 (MQ=255)
                                                                                                              tATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGaa                                         <  2:2342534/150‑1 (MQ=255)
                                                                                                               aTCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAg                                        <  1:1304987/150‑1 (MQ=255)
                                                                                                                   cTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGAGTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGgtgt                                    <  2:3105595/150‑1 (MQ=255)
                                                                                                                    tataACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGgtgtg                                   >  2:3682268/1‑150 (MQ=255)
                                                                                                                     ataACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACCTGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGgtgtgt                                  <  2:1938113/150‑1 (MQ=255)
                                                                                                                     ataACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGgtgtgt                                  <  1:3746650/150‑1 (MQ=255)
                                                                                                                       aaCCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGc                                >  1:4358959/1‑150 (MQ=255)
                                                                                                                          ccGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGt                             >  1:2277792/1‑150 (MQ=255)
                                                                                                                            ggATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTc                           >  1:2065118/1‑150 (MQ=255)
                                                                                                                                gTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGAc                       <  2:203533/150‑1 (MQ=255)
                                                                                                                                 ttGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACg                      <  1:3843546/150‑1 (MQ=255)
                                                                                                                                  tgatgaTGCTCCGCGTATGAAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACgg                     >  2:2227082/1‑150 (MQ=255)
                                                                                                                                   gatgatGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGt                    <  2:377272/150‑1 (MQ=255)
                                                                                                                                   gatgatGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGt                    <  1:3652847/150‑1 (MQ=255)
                                                                                                                                    atgatgCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTc                   >  2:4341651/1‑150 (MQ=255)
                                                                                                                                    atgatgCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTc                   >  2:4204559/1‑150 (MQ=255)
                                                                                                                                     tgatgCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCt                  <  1:3680911/150‑1 (MQ=255)
                                                                                                                                     tgatgCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCt                  <  1:3654931/150‑1 (MQ=255)
                                                                                                                                      gatgCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCAGACGGTCTg                 <  1:3453835/150‑1 (MQ=255)
                                                                                                                                          cTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACa             >  1:3784493/1‑150 (MQ=255)
                                                                                                                                             cgcgTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGa          <  1:3359247/150‑1 (MQ=255)
                                                                                                                                               cgTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAAc        <  1:241745/150‑1 (MQ=255)
                                                                                                                                                 tATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGTAAGATCCCTGCCTGTCGTTCCGCCGCTCCTGCCGTGTAGGTGGGTGCGGTTCCGACGGTCTGTACATGATCgg      <  2:680184/150‑1 (MQ=255)
                                                                                                                                                     cAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAg  >  1:1249517/1‑150 (MQ=255)
                                                                                                                                                     cAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAg  <  2:1556757/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAG  >  CP009273/753056‑753354

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: