Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,987,104 G→T 8.4% 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,1040GT 8.4% 1072.8 / 42.0 289R106S (CGC→AGC) cyaYiron‑dependent inhibitor of iron‑sulfur cluster formation; frataxin; iron‑binding and oxidizing protein
Reads supporting (aligned to +/- strand):  ref base G (121/141);  new base T (12/15);  total (133/156)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 292 total aligned reads displayed.

GATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCG  >  CP009273/3986955‑3987253
                                                                                                                                                     |                                                                                                                                                     
       gTGAACAGGTGATTCCGTTCCGCATAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAactg                                                                                                                                                <  2:3105052/150‑1 (MQ=255)
       gTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAactg                                                                                                                                                >  2:3202746/1‑150 (MQ=255)
        tGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAactga                                                                                                                                               <  1:1549480/150‑1 (MQ=255)
            cAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGt                                                                                                                                           <  1:3249762/150‑1 (MQ=255)
             aGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGtt                                                                                                                                          >  2:390977/1‑150 (MQ=255)
             aGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCAGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGtt                                                                                                                                          >  2:3455177/1‑150 (MQ=255)
               gTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTc                                                                                                                                        <  1:3873420/150‑1 (MQ=255)
                  atacagttccgCACAAAATCTATCGGTAACATTCCGCCAGCCAAACAGTATCACGATACGCAGCTATTACACCAAAATTTATCGAGATGACCCAGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTATCGGAAACTGACTGTATCAca                                                                                                                                     >  2:194972/6‑149 (MQ=255)
                  aTTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCAcc                                                                                                                                     <  2:3280604/150‑1 (MQ=255)
                  aTTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCAcc                                                                                                                                     >  2:562353/1‑150 (MQ=255)
                  aTTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCAcc                                                                                                                                     >  2:45729/1‑150 (MQ=255)
                   ttccgttccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCAccc                                                                                                                                    >  2:2601226/1‑150 (MQ=255)
                    tccgttccgCACAAACTCTATCGGTAACATGCCGCCTGCCAATCAGGATCACAATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCg                                                                                                                                   <  1:3147858/150‑1 (MQ=255)
                      cgttccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGcc                                                                                                                                 >  2:4162563/1‑150 (MQ=255)
                      cgttccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGcc                                                                                                                                 >  1:2306714/1‑150 (MQ=255)
                         tccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCctgc                                                                                                                              >  1:3555933/1‑150 (MQ=255)
                         tccgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCctgc                                                                                                                              >  1:1684649/1‑150 (MQ=255)
                           cgCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCctgctg                                                                                                                            <  1:86568/150‑1 (MQ=255)
                             cacaAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCTGAAACTGACTGTTTCACCCGCCTGCTGCg                                                                                                                          <  1:2793382/150‑1 (MQ=255)
                              acaAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGt                                                                                                                         <  1:4030601/150‑1 (MQ=255)
                              acaAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGt                                                                                                                         <  2:82499/150‑1 (MQ=255)
                                aaaaTCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTcg                                                                                                                       <  1:1611272/150‑1 (MQ=255)
                                     cTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCgccgcc                                                                                                                  <  1:4265003/150‑1 (MQ=255)
                                             aaCATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAg                                                                                                          <  1:1775032/150‑1 (MQ=255)
                                             aaCATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAg                                                                                                          >  2:1504897/1‑150 (MQ=255)
                                               cATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCa                                                                                                        >  1:526812/1‑150 (MQ=255)
                                                  gccgccTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTATTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCACTCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGATCCAGCAAAt                                                                                                     >  1:2234748/1‑150 (MQ=255)
                                                          ccAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAg                                                                                             >  2:225007/1‑150 (MQ=255)
                                                           cAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAgg                                                                                            >  1:3419130/1‑150 (MQ=255)
                                                            aaTCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGt                                                                                           >  2:2020356/1‑150 (MQ=255)
                                                              tCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGttt                                                                                         <  1:2842686/150‑1 (MQ=255)
                                                               cAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTc                                                                                        <  1:3156571/150‑1 (MQ=255)
                                                               cAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTc                                                                                        <  2:1425365/150‑1 (MQ=255)
                                                                aGGATCACGATACGCCNCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCNGCTGCGTCTCCGCCTGTTCCAGCAAATCCCAGAAGGTTTcg                                                                                       >  2:3591767/1‑150 (MQ=255)
                                                                aGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTcg                                                                                       >  1:1049909/1‑150 (MQ=255)
                                                                  gATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTcgcc                                                                                     <  2:1695487/150‑1 (MQ=255)
                                                                    tCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTcgccgc                                                                                   >  1:1231114/1‑150 (MQ=255)
                                                                    tCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTcgccgc                                                                                   >  1:1341321/1‑150 (MQ=255)
                                                                     cACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCt                                                                                  >  2:213698/1‑150 (MQ=255)
                                                                     cACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCt                                                                                  <  2:3715243/150‑1 (MQ=255)
                                                                     cACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCt                                                                                  <  1:4375272/150‑1 (MQ=255)
                                                                      aCGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTg                                                                                 <  1:3216005/150‑1 (MQ=255)
                                                                      aCGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTg                                                                                 <  1:1482877/150‑1 (MQ=255)
                                                                       cGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGc                                                                                <  2:1980028/150‑1 (MQ=255)
                                                                        gATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCTGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCg                                                                               <  2:2522324/150‑1 (MQ=255)
                                                                        gATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCg                                                                               >  1:562531/1‑150 (MQ=255)
                                                                         aTACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGa                                                                              >  1:3474206/1‑150 (MQ=255)
                                                                          tACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGAt                                                                             >  2:1383613/1‑150 (MQ=255)
                                                                            cgccgcTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATca                                                                           >  2:4404542/1‑150 (MQ=255)
                                                                            cgccgcTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATca                                                                           >  2:3333242/1‑150 (MQ=255)
                                                                            cgccgcTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATca                                                                           >  2:922127/1‑150 (MQ=255)
                                                                            cgccgcTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATca                                                                           >  2:2458528/1‑150 (MQ=255)
                                                                            cgccgcTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATca                                                                           <  1:1128738/150‑1 (MQ=255)
                                                                             gccgcTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATcac                                                                          <  1:3982126/150‑1 (MQ=255)
                                                                               cgcTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACaa                                                                        <  1:3825468/150‑1 (MQ=255)
                                                                                gcTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACaaa                                                                       >  1:3884595/1‑150 (MQ=255)
                                                                                  tATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCTGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGCAGGTTTCGCCGCTGCGATCACAAATc                                                                     >  2:112346/1‑150 (MQ=255)
                                                                                  tATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATc                                                                     <  2:496367/150‑1 (MQ=255)
                                                                                    ttACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCa                                                                   >  2:2902930/1‑150 (MQ=255)
                                                                                     tACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCTGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCAc                                                                  <  1:3942520/150‑1 (MQ=255)
                                                                                     tACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCAc                                                                  >  2:2193908/1‑150 (MQ=255)
                                                                                        aGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCa                                                               >  1:3188041/1‑150 (MQ=255)
                                                                                         gCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCAt                                                              <  2:3881314/150‑1 (MQ=255)
                                                                                              atTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCTGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCt                                                         >  2:2152671/1‑150 (MQ=255)
                                                                                                 tttCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTc                                                      <  2:1531735/150‑1 (MQ=255)
                                                                                                           aaCGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAAt                                            <  1:2621413/150‑1 (MQ=255)
                                                                                                           aaCGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAAt                                            >  2:2773860/1‑150 (MQ=255)
                                                                                                             cGTGCCGGAAAGCGATGCTTATCAGGCATGCAATATTAGCGGAAACTGACTGTTTCACCCGCATGCTGCGTCGCAGACTGTTCCAGCAAATCCCATAATGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTAAAAATgg                                          <  2:407487/150‑1 (MQ=255)
                                                                                                             cGTGCCGGAAAGCGAGGCTTATCCGGCTTGCAATCTTAGCGGATACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATgg                                          <  1:2567293/150‑1 (MQ=255)
                                                                                                                 ccGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAgc                                      <  1:650532/150‑1 (MQ=255)
                                                                                                                   ggAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAgccg                                    >  1:3022659/1‑150 (MQ=255)
                                                                                                                     aaaGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAgccgcc                                  >  1:2156394/1‑150 (MQ=255)
                                                                                                                      aaGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGccc                                 >  1:617104/1‑150 (MQ=255)
                                                                                                                      aaGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGccc                                 >  1:2784542/1‑150 (MQ=255)
                                                                                                                       aGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCt                                >  2:3426918/1‑150 (MQ=255)
                                                                                                                         cGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGt                              >  2:3045973/1‑150 (MQ=255)
                                                                                                                          caGGCTTATCCGGCATGCAAACATAGCGGAAACAGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAAACCAAGAAGGTTACGCCGCTGCGATCACAAATCCACTCATCGCCATTCAGGTCAAAATGGTAGCCGCCCTgat                             >  2:2351639/2‑148 (MQ=255)
                                                                                                                           aGGCTTATCCGGCATGCAATCTTAGCTGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGttt                            <  1:3098772/150‑1 (MQ=255)
                                                                                                                           aGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGttt                            >  1:1915351/1‑150 (MQ=255)
                                                                                                                           aGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGttt                            <  1:390977/150‑1 (MQ=255)
                                                                                                                            ggCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTg                           <  2:1670408/150‑1 (MQ=255)
                                                                                                                              cTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGt                         >  2:3491738/1‑150 (MQ=255)
                                                                                                                              cTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGt                         <  2:3938530/150‑1 (MQ=255)
                                                                                                                              cTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGt                         <  2:3691301/150‑1 (MQ=255)
                                                                                                                                tATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTg                       <  1:2400693/150‑1 (MQ=255)
                                                                                                                                  tCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTgcc                     <  1:2864683/150‑1 (MQ=255)
                                                                                                                                    cGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTgccag                   <  1:2106331/150‑1 (MQ=255)
                                                                                                                                    cGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTgccag                   <  1:4162563/150‑1 (MQ=255)
                                                                                                                                      gCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTgccagcc                 <  1:2728044/150‑1 (MQ=255)
                                                                                                                                       cATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTgccagcca                <  2:3876978/150‑1 (MQ=255)
                                                                                                                                        aTGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCAt               <  2:2507340/150‑1 (MQ=255)
                                                                                                                                         tGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATa              <  1:996507/150‑1 (MQ=255)
                                                                                                                                          gCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATAc             >  2:2278562/1‑150 (MQ=255)
                                                                                                                                          gCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATAc             <  2:4285145/150‑1 (MQ=255)
                                                                                                                                           cAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATAcc            <  2:1100064/150‑1 (MQ=255)
                                                                                                                                            aaTCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCt           >  2:38108/1‑150 (MQ=255)
                                                                                                                                                  aGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCa     <  2:4511121/150‑1 (MQ=255)
                                                                                                                                                  aGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCa     >  1:184807/1‑150 (MQ=255)
                                                                                                                                                   gCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAg    <  1:3359399/150‑1 (MQ=255)
                                                                                                                                                   gCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAg    >  1:1283632/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCGCCGCCTGTTCCAGCAAATCCCAGAAGGTTTCGCCGCTGCGATCACAAATCCACTCATCGCCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCG  >  CP009273/3986955‑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: