Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,936,836 C→A 100% A47A (GCC→GCA csdA → cysteine sulfinate desulfinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,936,8360CA99.7% 990.6 / ‑4.5 303A47A (GCC→GCAcsdAcysteine sulfinate desulfinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (119/183);  total (119/184)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 303 total aligned reads displayed.

GAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAACAT  >  CP009273/2936689‑2936982
                                                                                                                                                   |                                                                                                                                                  
    aCCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAAc                                                                                                                                              >  2:870333/1‑150 (MQ=255)
    aCCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAAc                                                                                                                                              <  2:32093/150‑1 (MQ=255)
    aCCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAAc                                                                                                                                              <  1:1556152/150‑1 (MQ=255)
    aCCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTGTACAGTCTGAGCGCAGGAAAc                                                                                                                                              >  2:350201/1‑150 (MQ=255)
     ccATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACg                                                                                                                                             <  1:1222386/150‑1 (MQ=255)
     ccATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACg                                                                                                                                             <  2:4782817/150‑1 (MQ=255)
        tGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTcc                                                                                                                                          >  2:1185877/1‑150 (MQ=255)
        tGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTcc                                                                                                                                          <  2:5029740/150‑1 (MQ=255)
         gAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCa                                                                                                                                         <  1:4415770/150‑1 (MQ=255)
          aaCGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCAt                                                                                                                                        >  2:430857/1‑150 (MQ=255)
             gTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGANACGTCCATCGc                                                                                                                                     <  1:2758004/150‑1 (MQ=255)
             gTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGc                                                                                                                                     <  2:2203736/150‑1 (MQ=255)
               ttttAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAg                                                                                                                                   <  1:4723728/150‑1 (MQ=255)
                tttAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGc                                                                                                                                  >  1:212533/1‑150 (MQ=255)
                 ttAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGcc                                                                                                                                 >  2:758034/1‑150 (MQ=255)
                  tAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCa                                                                                                                                <  2:331165/150‑1 (MQ=255)
                     tCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGtt                                                                                                                             <  2:3454568/150‑1 (MQ=255)
                      cccGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGttt                                                                                                                            <  2:4474586/150‑1 (MQ=255)
                      cccGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGttt                                                                                                                            <  1:3329696/150‑1 (MQ=255)
                             aGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAg                                                                                                                     <  2:2493877/150‑1 (MQ=255)
                             aGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAg                                                                                                                     >  2:4277525/1‑150 (MQ=255)
                               tttCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGcc                                                                                                                   <  1:2573916/150‑1 (MQ=255)
                               tttCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGcc                                                                                                                   <  2:2013720/150‑1 (MQ=255)
                                   gcgcCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAAc                                                                                                               <  1:1495923/150‑1 (MQ=255)
                                    cgcCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACg                                                                                                              >  1:4180317/1‑150 (MQ=255)
                                     gcCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGc                                                                                                             >  1:2948719/1‑150 (MQ=255)
                                     gcCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGAGCTTAAACATGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGc                                                                                                             >  1:4851798/1‑150 (MQ=255)
                                      cccAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGcc                                                                                                            >  2:3897375/1‑150 (MQ=255)
                                        cAGGTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTg                                                                                                          <  2:4897260/150‑1 (MQ=255)
                                          gTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGAc                                                                                                        >  1:5125900/1‑150 (MQ=255)
                                                  ctcTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGtt                                                                                                <  1:4983924/148‑1 (MQ=255)
                                                  cACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGtt                                                                                                <  2:2663274/150‑1 (MQ=255)
                                                  cACGACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGtt                                                                                                <  2:2114055/150‑1 (MQ=255)
                                                     tACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCCACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTGTGCCGAAGCCCAACGCCTGACCGCGCGTTATg                                                                                             <  2:2293200/150‑1 (MQ=255)
                                                     tACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATg                                                                                             >  2:4278860/1‑150 (MQ=255)
                                                     tACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATg                                                                                             <  2:1159015/150‑1 (MQ=255)
                                                      aCAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGa                                                                                            <  2:1629512/150‑1 (MQ=255)
                                                           ctGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGACGCTg                                                                                       <  2:837686/149‑1 (MQ=255)
                                                            tGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGc                                                                                      >  2:2391473/1‑150 (MQ=255)
                                                                ggCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACga                                                                                  <  1:244081/150‑1 (MQ=255)
                                                                ggCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACga                                                                                  >  1:429762/1‑150 (MQ=255)
                                                                   gTCTATCTCGACAGCGCCGCGACCGCTCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACgagag                                                                               <  2:2749015/150‑1 (MQ=255)
                                                                   gTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACgagag                                                                               >  1:182537/1‑150 (MQ=255)
                                                                    tCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACgagaga                                                                              <  2:4113329/150‑1 (MQ=255)
                                                                    tCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACgagaga                                                                              <  1:2999504/150‑1 (MQ=255)
                                                                    tCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACgagaga                                                                              <  2:4690131/150‑1 (MQ=255)
                                                                     cTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGaa                                                                             >  1:4675286/1‑150 (MQ=255)
                                                                      tATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCTCGCGTTATGAAGCTGCACGAGAGaaa                                                                            >  1:4461615/1‑150 (MQ=255)
                                                                       aTCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAg                                                                           <  2:4045860/150‑1 (MQ=255)
                                                                       aTCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTGTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAg                                                                           <  2:1409714/150‑1 (MQ=255)
                                                                       aTCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCGGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCCGACCGCGCGTTATGAAGCTGCACGAGAGAAAg                                                                           >  1:4786405/1‑150 (MQ=255)
                                                                         ctcGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTg                                                                         <  1:1146878/150‑1 (MQ=255)
                                                                         ctcGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTg                                                                         <  1:4891505/150‑1 (MQ=255)
                                                                         ctcGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTg                                                                         <  2:3787414/150‑1 (MQ=255)
                                                                         ctcGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTg                                                                         <  1:3124428/150‑1 (MQ=255)
                                                                         ctcGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTg                                                                         <  1:2113477/150‑1 (MQ=255)
                                                                               aGCGCCGCGACCGCGATTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCaa                                                                   <  2:3770665/150‑1 (MQ=255)
                                                                                 cgccgcGACCGCGCTTAACCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAAtt                                                                 <  1:73198/150‑1 (MQ=255)
                                                                                 cgccgcGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAAtt                                                                 <  1:3909053/150‑1 (MQ=255)
                                                                                  gccgcGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTa                                                                <  2:3734783/150‑1 (MQ=255)
                                                                                    cgcgACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACt                                                              <  1:1313112/150‑1 (MQ=255)
                                                                                    cgcgACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTAACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACt                                                              <  1:4386444/150‑1 (MQ=255)
                                                                                     gcgACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTg                                                             <  1:2280256/150‑1 (MQ=255)
                                                                                        aCCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAAt                                                          <  1:85433/150‑1 (MQ=255)
                                                                                        aCCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAAt                                                          <  1:4643161/150‑1 (MQ=255)
                                                                                         ccGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATg                                                         <  1:4000937/150‑1 (MQ=255)
                                                                                         ccGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATg                                                         <  1:373076/150‑1 (MQ=255)
                                                                                           gcgcTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCa                                                       <  2:2171898/150‑1 (MQ=255)
                                                                                            cgcTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCAc                                                      <  1:5080683/150‑1 (MQ=255)
                                                                                             gcTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCAcc                                                     >  2:775191/1‑150 (MQ=255)
                                                                                             gcTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCAcc                                                     <  2:4470113/150‑1 (MQ=255)
                                                                                                 aaaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGgat                                                 <  1:1029013/150‑1 (MQ=255)
                                                                                                  aaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGgatg                                                <  2:177598/150‑1 (MQ=255)
                                                                                                  aaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGgatg                                                <  1:1356595/150‑1 (MQ=255)
                                                                                                  aaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGgatg                                                <  2:4190414/150‑1 (MQ=255)
                                                                                                     cTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATa                                             <  1:4269338/150‑1 (MQ=255)
                                                                                                        aaGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATaaaa                                          <  2:1574141/150‑1 (MQ=255)
                                                                                                         aGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAAc                                         <  1:97798/150‑1 (MQ=255)
                                                                                                            cGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTAt                                      >  1:899925/1‑150 (MQ=255)
                                                                                                            cGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTAt                                      >  1:3035847/1‑150 (MQ=255)
                                                                                                             gTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATc                                     <  1:2221958/150‑1 (MQ=255)
                                                                                                                   gAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTgg                               >  2:5285111/1‑150 (MQ=255)
                                                                                                                   gAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTgg                               >  2:4857425/1‑150 (MQ=255)
                                                                                                                    aaGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGa                              >  2:5203169/1‑150 (MQ=255)
                                                                                                                     aGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGAc                             >  1:3252762/1‑150 (MQ=255)
                                                                                                                       ccACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGAcgc                           >  1:1093003/1‑150 (MQ=255)
                                                                                                                        cACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGAcgcg                          >  2:4835323/1‑150 (MQ=255)
                                                                                                                            cAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCACTTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGc                      >  1:79171/1‑150 (MQ=255)
                                                                                                                             aaCAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGca                     >  1:11145/1‑150 (MQ=255)
                                                                                                                                aGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGcacca                  >  1:3331447/1‑150 (MQ=255)
                                                                                                                                  ttttACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGATAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACt                <  1:613737/150‑1 (MQ=255)
                                                                                                                                  ttttACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACt                <  1:2539282/150‑1 (MQ=255)
                                                                                                                                  ttttACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACt                <  2:2916408/150‑1 (MQ=255)
                                                                                                                                    ttACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGa              <  1:3470854/150‑1 (MQ=255)
                                                                                                                                      aCAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAAt            <  2:3022946/150‑1 (MQ=255)
                                                                                                                                        aGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATcc          <  2:3427764/150‑1 (MQ=255)
                                                                                                                                           cTGAGCGCAGGAAACGTCCATCGAAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAATCTGCACGAGCGAAAGTGGCGCAATTACTGAAAGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGATTCCATc       >  2:4992983/1‑150 (MQ=255)
                                                                                                                                            tGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCa      >  1:1438364/1‑150 (MQ=255)
                                                                                                                                             gAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCaa     >  2:17942/1‑150 (MQ=255)
                                                                                                                                               gcgcAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAACa   <  1:4112566/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAACAT  >  CP009273/2936689‑2936982

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: