Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 753,549 G→T 95.0% Q135H (CAG→CAT sdhB → succinate dehydrogenase, FeS subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273753,5490GT95.0% 653.3 / 24.1 230Q135H (CAG→CATsdhBsuccinate dehydrogenase, FeS subunit
Reads supporting (aligned to +/- strand):  ref base G (4/8);  new base T (104/114);  total (108/122)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 231 total aligned reads displayed.

AGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATT  >  CP009273/753404‑753693
                                                                                                                                                 |                                                                                                                                                
 gaagaTTGTGATTCGCCCGCTGCCATGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATg                                                                                                                                               <  1:3424334/148‑1 (MQ=255)
 gaagaTTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATg                                                                                                                                               >  2:1546042/1‑148 (MQ=255)
  aagaTTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGc                                                                                                                                              <  2:345358/148‑1 (MQ=255)
   agaTTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGcc                                                                                                                                             <  2:1635442/148‑1 (MQ=255)
      ttGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCaga                                                                                                                                          >  1:2251559/1‑148 (MQ=255)
      ttGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCaga                                                                                                                                          >  1:3023051/1‑148 (MQ=255)
       tgtgATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCagag                                                                                                                                         <  2:2446600/148‑1 (MQ=255)
            ttCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAgcg                                                                                                                                    <  2:1995203/148‑1 (MQ=255)
                  cGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGaaaa                                                                                                                              <  1:3052161/148‑1 (MQ=255)
                        cAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGa                                                                                                                        <  2:470560/148‑1 (MQ=255)
                        cAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGa                                                                                                                        <  2:1912401/148‑1 (MQ=255)
                         aGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGAc                                                                                                                       >  1:3055134/1‑148 (MQ=255)
                           gTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACgg                                                                                                                     >  1:431462/1‑148 (MQ=255)
                           gTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACgg                                                                                                                     >  1:2760565/1‑148 (MQ=255)
                           gTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACgg                                                                                                                     >  1:563944/1‑148 (MQ=255)
                           gTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACgg                                                                                                                     >  1:2352613/1‑148 (MQ=255)
                            tttACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACggg                                                                                                                    <  2:3195765/148‑1 (MQ=255)
                                    tGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGa                                                                                                            >  1:986994/1‑148 (MQ=255)
                                     gATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGaa                                                                                                           <  1:1309399/148‑1 (MQ=255)
                                      aTCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAAt                                                                                                          >  1:191113/1‑148 (MQ=255)
                                      aTCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAAt                                                                                                          <  2:1255959/148‑1 (MQ=255)
                                       tCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATg                                                                                                         <  1:2952536/148‑1 (MQ=255)
                                        ccGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGt                                                                                                        >  2:774778/1‑148 (MQ=255)
                                         cgcgATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTa                                                                                                       <  1:345128/148‑1 (MQ=255)
                                         cgcgATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTa                                                                                                       <  1:1429458/148‑1 (MQ=255)
                                            gATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATtc                                                                                                    <  1:1621186/148‑1 (MQ=255)
                                             aTTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATtct                                                                                                   <  1:1938440/148‑1 (MQ=255)
                                              tttGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATtctc                                                                                                  >  2:2740739/1‑148 (MQ=255)
                                                tggtggTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTg                                                                                                <  1:3363241/148‑1 (MQ=255)
                                                     gtAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCAt                                                                                           <  2:969906/148‑1 (MQ=255)
                                                     gtAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCAt                                                                                           >  1:1537925/1‑148 (MQ=255)
                                                     gtAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCAt                                                                                           <  1:967922/148‑1 (MQ=255)
                                                        gACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCAtgtt                                                                                        <  2:3253063/148‑1 (MQ=255)
                                                         aCATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCAtgttg                                                                                       >  1:1876336/1‑148 (MQ=255)
                                                            tGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTc                                                                                    >  2:159605/1‑148 (MQ=255)
                                                               gACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAAc                                                                                 >  1:1864440/1‑148 (MQ=255)
                                                               gACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAAc                                                                                 >  1:2397005/1‑148 (MQ=255)
                                                                aCAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAAcc                                                                                >  1:1036612/1‑148 (MQ=255)
                                                                   aTTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACctct                                                                             >  2:620313/1‑148 (MQ=255)
                                                                   aTTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACctct                                                                             >  2:1430329/1‑148 (MQ=255)
                                                                    ttCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACACAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCtt                                                                            <  2:3065001/148‑1 (MQ=255)
                                                                      cTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGt                                                                          <  2:2833507/148‑1 (MQ=255)
                                                                      cTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGt                                                                          <  1:2743693/148‑1 (MQ=255)
                                                                       tATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTc                                                                         <  2:395028/148‑1 (MQ=255)
                                                                       tATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTc                                                                         >  2:3491906/1‑148 (MQ=255)
                                                                       tATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTc                                                                         <  2:2739799/148‑1 (MQ=255)
                                                                          gcgcAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGt                                                                      >  1:230062/1‑148 (MQ=255)
                                                                               atatGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTc                                                                 >  2:1164556/1‑148 (MQ=255)
                                                                               atatGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTc                                                                 <  1:2064371/148‑1 (MQ=255)
                                                                                 atGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCtg                                                               >  2:1052582/1‑148 (MQ=255)
                                                                                  tGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCtgg                                                              >  2:1105477/1‑148 (MQ=255)
                                                                                  tGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCtgg                                                              <  2:3222458/148‑1 (MQ=255)
                                                                                   gagaAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCtggt                                                             <  2:1439291/148‑1 (MQ=255)
                                                                                   gagaAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCtggt                                                             <  2:511359/148‑1 (MQ=255)
                                                                                   gagaAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCtggt                                                             <  2:3139615/148‑1 (MQ=255)
                                                                                    agaAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCtggtg                                                            >  1:2895593/1‑148 (MQ=255)
                                                                                      aaaaTTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGa                                                          >  1:1859738/1‑148 (MQ=255)
                                                                                       aaaTTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGaa                                                         >  2:362154/1‑148 (MQ=255)
                                                                                        aaTTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAAt                                                        >  1:706101/1‑148 (MQ=255)
                                                                                        aaTTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAAt                                                        >  1:2733415/1‑148 (MQ=255)
                                                                                           tAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATccc                                                     >  1:15779/1‑148 (MQ=255)
                                                                                           tAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATccc                                                     >  1:3494392/1‑148 (MQ=255)
                                                                                           tAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATccc                                                     >  2:1769871/1‑148 (MQ=255)
                                                                                            aaGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCg                                                    <  1:1204636/148‑1 (MQ=255)
                                                                                            aaGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCg                                                    >  2:2502575/1‑148 (MQ=255)
                                                                                                 ttACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATCTGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAg                                               <  1:755154/148‑1 (MQ=255)
                                                                                                  tACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGt                                              >  2:3114512/1‑148 (MQ=255)
                                                                                                  tACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGt                                              >  2:2090054/1‑148 (MQ=255)
                                                                                                  tACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGt                                              >  2:294727/1‑148 (MQ=255)
                                                                                                      tGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTAt                                          <  1:814466/148‑1 (MQ=255)
                                                                                                           aataatGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGcc                                     <  1:3256992/148‑1 (MQ=255)
                                                                                                            ataatGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGTAATCCCGATAAGTTTATCGGccc                                    >  2:2317816/1‑148 (MQ=255)
                                                                                                               atGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGc                                 >  2:1507263/1‑148 (MQ=255)
                                                                                                                tGGACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCa                                >  2:519255/1‑148 (MQ=255)
                                                                                                                 ggACAAAATCTGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAg                               >  2:1546963/1‑148 (MQ=255)
                                                                                                                  gACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAgg                              >  1:1046802/1‑148 (MQ=255)
                                                                                                                  gACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAgg                              <  2:2251559/148‑1 (MQ=255)
                                                                                                                  gACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAgg                              <  2:1698628/148‑1 (MQ=255)
                                                                                                                  gACAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAgg                              <  2:3023051/148‑1 (MQ=255)
                                                                                                                  gACAAAATCCGCCAGCTCGCGAGCATGTACATAGGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAgg                              <  2:207386/148‑1 (MQ=255)
                                                                                                                   aCAAAATCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGc                             <  2:1705231/148‑1 (MQ=255)
                                                                                                                     aaaaTCCGCCAGCTCGCGAGCATTTACATCTGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCtt                           <  1:1391978/148‑1 (MQ=255)
                                                                                                                     aaaaTCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCtt                           <  1:2460927/148‑1 (MQ=255)
                                                                                                                     aaaaTCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCtt                           >  2:2660541/1‑148 (MQ=255)
                                                                                                                      aaaTCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTg                          >  2:1985210/1‑148 (MQ=255)
                                                                                                                      aaaTCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTg                          <  2:1760021/148‑1 (MQ=255)
                                                                                                                       aagccGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGt                         <  1:1138535/145‑1 (MQ=255)
                                                                                                                       aaTCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGt                         >  1:20319/1‑148 (MQ=255)
                                                                                                                         tCCGCCAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTa                       >  1:2807332/1‑148 (MQ=255)
                                                                                                                             ccAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCgg                   >  2:1342896/1‑148 (MQ=255)
                                                                                                                              cAGCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGc                  <  2:2895641/148‑1 (MQ=255)
                                                                                                                                gCTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCat                >  1:543085/1‑148 (MQ=255)
                                                                                                                                 cTCGCGAGCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCata               >  2:1985856/1‑148 (MQ=255)
                                                                                                                                        gCATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTc        <  1:696777/148‑1 (MQ=255)
                                                                                                                                         cATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTcc       <  1:3287033/148‑1 (MQ=255)
                                                                                                                                         cATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTcc       <  1:3286977/148‑1 (MQ=255)
                                                                                                                                         cATTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTcc       <  1:299292/148‑1 (MQ=255)
                                                                                                                                          aTTTACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTGGTTAGCGGCAGAGCGTTTCCt      >  2:3365151/1‑148 (MQ=255)
                                                                                                                                           tttACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCtg     <  2:706101/148‑1 (MQ=255)
                                                                                                                                             tACATATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCtgat   <  2:2352613/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATT  >  CP009273/753404‑753693

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: