Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,987,104 G→T 12.3% R106S (CGC→AGC)  cyaY ← iron‑dependent inhibitor of iron‑sulfur cluster formation; frataxin; iron‑binding and oxidizing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,987,1040GT12.3% 984.7 / 74.1 281R106S (CGC→AGC) cyaYiron‑dependent inhibitor of iron‑sulfur cluster formation; frataxin; iron‑binding and oxidizing protein
Reads supporting (aligned to +/- strand):  ref base G (130/116);  new base T (25/10);  total (155/126)
Fisher's exact test for biased strand distribution p-value = 4.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 288 total aligned reads displayed.

ATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCG  >  CP009273/3986956‑3987253
                                                                                                                                                    |                                                                                                                                                     
aTGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCgg                                                                                                                                                      <  2:1028949/150‑1 (MQ=255)
     cGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAact                                                                                                                                                 <  1:4418158/150‑1 (MQ=255)
      gTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAactg                                                                                                                                                <  1:306742/150‑1 (MQ=255)
      gTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAactg                                                                                                                                                >  1:2346933/1‑150 (MQ=255)
        gAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAactgac                                                                                                                                              >  2:4335576/1‑150 (MQ=255)
          aCAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAactgactg                                                                                                                                            >  2:1639895/1‑150 (MQ=255)
            aGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGtt                                                                                                                                          <  1:3889306/150‑1 (MQ=255)
            aGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGtt                                                                                                                                          >  1:903004/1‑150 (MQ=255)
             ggTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCTGAAACTGACTGttt                                                                                                                                         <  2:1283221/150‑1 (MQ=255)
                gATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCAc                                                                                                                                      <  2:2038511/150‑1 (MQ=255)
                 aTTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCAcc                                                                                                                                     <  2:4489665/150‑1 (MQ=255)
                  ttccgttccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCAccc                                                                                                                                    <  2:1253134/150‑1 (MQ=255)
                  ttccgttccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCAccc                                                                                                                                    <  1:2423867/150‑1 (MQ=255)
                   tccgttccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCg                                                                                                                                   >  2:2652986/1‑150 (MQ=255)
                      gttccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCct                                                                                                                                >  1:1967038/1‑150 (MQ=255)
                      gttccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCct                                                                                                                                >  2:4252675/1‑150 (MQ=255)
                          cgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCTGAAACTGACTGTTTCACCCGCctgctg                                                                                                                            <  2:1510928/150‑1 (MQ=255)
                          cgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCctgctg                                                                                                                            >  1:67041/1‑150 (MQ=255)
                             acaAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGt                                                                                                                         <  1:3959329/150‑1 (MQ=255)
                               aaaaTCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTcg                                                                                                                       >  1:1183868/1‑150 (MQ=255)
                                 aaTCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTcgcc                                                                                                                     <  1:2180245/150‑1 (MQ=255)
                                   tCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTcgccgc                                                                                                                   <  1:2646363/150‑1 (MQ=255)
                                     tATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCCATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCt                                                                                                                 <  2:789696/150‑1 (MQ=255)
                                     tATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCt                                                                                                                 <  1:2271093/150‑1 (MQ=255)
                                       tCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGt                                                                                                               >  2:693102/1‑150 (MQ=255)
                                          gTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTcc                                                                                                            >  1:1917164/1‑150 (MQ=255)
                                           tAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCTGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCa                                                                                                           >  2:1479225/1‑150 (MQ=255)
                                              cATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCa                                                                                                        >  2:2006211/1‑150 (MQ=255)
                                              cATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCa                                                                                                        <  2:1441570/150‑1 (MQ=255)
                                                  ccgccTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCTGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATc                                                                                                    >  1:3650826/1‑150 (MQ=255)
                                                  ccgccTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATc                                                                                                    >  1:3956499/1‑150 (MQ=255)
                                                   cgccTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATcc                                                                                                   <  2:3888848/150‑1 (MQ=255)
                                                    gccTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATccc                                                                                                  >  2:2718528/1‑150 (MQ=255)
                                                       tGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCAGCCTGTTCCAGCAAATCCCAGa                                                                                               >  2:3455004/1‑150 (MQ=255)
                                                        gCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGaa                                                                                              >  1:2346849/1‑150 (MQ=255)
                                                         ccAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAg                                                                                             <  1:2962543/150‑1 (MQ=255)
                                                          cAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAgg                                                                                            >  1:1276944/1‑150 (MQ=255)
                                                             tCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGttt                                                                                         <  1:4451677/150‑1 (MQ=255)
                                                               aGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCTGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTcg                                                                                       <  1:471380/150‑1 (MQ=255)
                                                                   tCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTcgccgc                                                                                   >  2:4403643/1‑150 (MQ=255)
                                                                      cGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGc                                                                                >  2:1792811/1‑150 (MQ=255)
                                                                       gATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCTGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCg                                                                               >  2:3777677/1‑150 (MQ=255)
                                                                       gATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCg                                                                               <  1:1266248/150‑1 (MQ=255)
                                                                          aCGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGCTc                                                                            >  2:2474352/1‑150 (MQ=255)
                                                                           cgccgcTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATca                                                                           >  1:1721558/1‑150 (MQ=255)
                                                                            gccgcTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATcac                                                                          >  2:4599339/1‑150 (MQ=255)
                                                                                cTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCTGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAAt                                                                      >  1:2528033/1‑150 (MQ=255)
                                                                                cTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAAt                                                                      >  1:3559973/1‑150 (MQ=255)
                                                                                      cAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTc                                                                >  1:2461261/1‑150 (MQ=255)
                                                                                         cAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATc                                                             <  1:4335576/150‑1 (MQ=255)
                                                                                           atatTTTTCGTGATGATCGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCGTAGCGGAAACTGACGGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCCTCGc                                                           >  2:4183242/1‑150 (MQ=255)
                                                                                             atTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCTGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCt                                                         >  2:3152957/1‑150 (MQ=255)
                                                                                             atTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCt                                                         <  1:2710553/150‑1 (MQ=255)
                                                                                             atTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCt                                                         <  2:1137678/150‑1 (MQ=255)
                                                                                             atTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCt                                                         <  2:893178/150‑1 (MQ=255)
                                                                                               ttttCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCttt                                                       >  1:3020693/1‑150 (MQ=255)
                                                                                               ttttCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCttt                                                       >  1:3584174/1‑150 (MQ=255)
                                                                                                tttCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTc                                                      >  1:1056377/1‑150 (MQ=255)
                                                                                                tttCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTc                                                      >  1:3379598/1‑150 (MQ=255)
                                                                                                    gTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGt                                                  <  1:847350/150‑1 (MQ=255)
                                                                                                     tgatgaACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTc                                                 >  2:4014190/1‑150 (MQ=255)
                                                                                                     tgatgaACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTc                                                 <  2:3297017/150‑1 (MQ=255)
                                                                                                     tgatgaACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTc                                                 <  1:3185844/150‑1 (MQ=255)
                                                                                                     tgatgaACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTc                                                 <  2:3085928/150‑1 (MQ=255)
                                                                                                     tgatgaACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTc                                                 <  1:4075177/150‑1 (MQ=255)
                                                                                                      gatgaACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCa                                                <  1:101031/150‑1 (MQ=255)
                                                                                                      gatgaACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCa                                                <  1:3602681/150‑1 (MQ=255)
                                                                                                         gaACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCTGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCaaaa                                             >  1:1084648/1‑150 (MQ=255)
                                                                                                         gaACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCaaaa                                             >  1:3686358/1‑150 (MQ=255)
                                                                                                         gaACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCaaaa                                             >  1:4463903/1‑150 (MQ=255)
                                                                                                         gaACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCaaaa                                             >  1:2768465/1‑150 (MQ=255)
                                                                                                         gaACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCaaaa                                             >  1:515172/1‑150 (MQ=255)
                                                                                                         gaACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCaaaa                                             <  1:385361/150‑1 (MQ=255)
                                                                                                           aCGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATg                                           <  2:3142315/150‑1 (MQ=255)
                                                                                                                 cGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCTGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAgcc                                     >  2:65177/1‑150 (MQ=255)
                                                                                                                 cGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAgcc                                     >  2:2367415/1‑150 (MQ=255)
                                                                                                                  ggAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAgccg                                    >  1:2387230/1‑150 (MQ=255)
                                                                                                                     aaGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGccc                                 >  1:1519487/1‑150 (MQ=255)
                                                                                                                      aGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCt                                >  2:2773178/1‑150 (MQ=255)
                                                                                                                       gCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTg                               >  2:210616/1‑150 (MQ=255)
                                                                                                                           ggCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTg                           <  1:3032456/150‑1 (MQ=255)
                                                                                                                           ggCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTg                           <  1:4575659/150‑1 (MQ=255)
                                                                                                                            gCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTgg                          <  2:3630509/150‑1 (MQ=255)
                                                                                                                              ttATCCGGCATGCAATCTTAGCTGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGtt                        >  1:1309331/1‑150 (MQ=255)
                                                                                                                              ttATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGtt                        <  2:3855285/150‑1 (MQ=255)
                                                                                                                               tATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTg                       >  2:1854759/1‑150 (MQ=255)
                                                                                                                                 tCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTgcc                     <  2:3215966/150‑1 (MQ=255)
                                                                                                                                 tCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTgcc                     <  2:2324317/150‑1 (MQ=255)
                                                                                                                                  ccGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTgcca                    <  1:573498/150‑1 (MQ=255)
                                                                                                                                     gCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTgccagcc                 >  2:3439423/1‑150 (MQ=255)
                                                                                                                                      cATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTgccagcca                >  1:2785219/1‑150 (MQ=255)
                                                                                                                                       aTGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCAt               <  1:1539117/150‑1 (MQ=255)
                                                                                                                                        tGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATa              <  2:1302294/150‑1 (MQ=255)
                                                                                                                                                 aGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCa     <  2:2747380/150‑1 (MQ=255)
                                                                                                                                                 aGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCa     >  1:1123770/1‑150 (MQ=255)
                                                                                                                                                  gCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAg    <  2:2241894/150‑1 (MQ=255)
                                                                                                                                                   cGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGTGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGc   <  1:1759459/150‑1 (MQ=255)
                                                                                                                                                   cGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGc   >  1:749511/1‑150 (MQ=255)
                                                                                                                                                    ggAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCg  >  1:2094536/1‑150 (MQ=255)
                                                                                                                                                    ggAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCg  <  2:63247/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCG  >  CP009273/3986956‑3987253

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: