Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 862,542 A→T 15.6% E173D (GAA→GAT)  iaaA → Isoaspartyl peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273862,5420AT15.6% 738.2 / 117.6 256E173D (GAA→GAT) iaaAIsoaspartyl peptidase
Reads supporting (aligned to +/- strand):  ref base A (98/118);  new base T (26/14);  total (124/132)
Fisher's exact test for biased strand distribution p-value = 2.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 269 total aligned reads displayed.

GCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGC  >  CP009273/862394‑862691
                                                                                                                                                    |                                                                                                                                                     
gCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGaaa                                                                                                                                                      <  1:3655156/150‑1 (MQ=255)
 cGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCCGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGataa                                                                                                                                                     >  1:4163932/1‑150 (MQ=255)
 cGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGaaaa                                                                                                                                                     >  2:1068817/1‑150 (MQ=255)
 cGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGaaaa                                                                                                                                                     >  2:3351627/1‑150 (MQ=255)
  gAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGaaaaa                                                                                                                                                    >  1:1011367/1‑150 (MQ=255)
   aaGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGATAAAc                                                                                                                                                   >  1:2865053/1‑150 (MQ=255)
   aaGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAAc                                                                                                                                                   <  1:964528/150‑1 (MQ=255)
         cAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGaaaa                                                                                                                                             >  1:1378825/1‑150 (MQ=255)
         cAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGaaaa                                                                                                                                             >  1:1351824/1‑150 (MQ=255)
         cAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGaaaa                                                                                                                                             >  1:3636822/1‑150 (MQ=255)
         cAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGaaaa                                                                                                                                             >  1:3636504/1‑150 (MQ=255)
         cAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGaaaa                                                                                                                                             >  2:3481911/1‑150 (MQ=255)
            aaaaTTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATgg                                                                                                                                          >  1:879515/1‑150 (MQ=255)
                   tGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGt                                                                                                                                   <  2:3330633/150‑1 (MQ=255)
                     cGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTgg                                                                                                                                 <  2:3461840/150‑1 (MQ=255)
                      gTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTggg                                                                                                                                <  1:648162/150‑1 (MQ=255)
                       tttGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTgggg                                                                                                                               <  1:1824660/150‑1 (MQ=255)
                       tttGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTgggg                                                                                                                               <  2:377968/150‑1 (MQ=255)
                       tttGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTgggg                                                                                                                               <  1:3302156/150‑1 (MQ=255)
                        ttGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTggggg                                                                                                                              >  2:48344/1‑150 (MQ=255)
                                  cATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGtt                                                                                                                    <  1:3978455/150‑1 (MQ=255)
                                    tGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTgg                                                                                                                  <  1:2585869/150‑1 (MQ=255)
                                     ggAGCGCGTCTCGCCGGAGATTTTCTCCAAGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGTCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGAGGGGGCCGTGGCGTTGGa                                                                                                                 >  1:214116/1‑150 (MQ=255)
                                        gcgcGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGATAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGAttt                                                                                                              <  2:711903/150‑1 (MQ=255)
                                        gcgcGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGAttt                                                                                                              <  2:4268038/150‑1 (MQ=255)
                                         cgcgTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTa                                                                                                             <  1:2273042/150‑1 (MQ=255)
                                           cgTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAACCTGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGa                                                                                                           <  1:109154/150‑1 (MQ=255)
                                           cgTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACGGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGa                                                                                                           <  1:2375174/150‑1 (MQ=255)
                                              ctcGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGATAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACgg                                                                                                        >  1:2948439/1‑150 (MQ=255)
                                                   cGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAAtt                                                                                                   >  1:3034084/1‑150 (MQ=255)
                                                     gagaTTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTg                                                                                                 <  2:826724/150‑1 (MQ=255)
                                                          tttCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTggcggc                                                                                            >  2:3942401/1‑150 (MQ=255)
                                                           ttCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCa                                                                                           <  2:2014294/150‑1 (MQ=255)
                                                             ctcCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGATAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGc                                                                                         >  2:2154121/1‑150 (MQ=255)
                                                             ctcCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGc                                                                                         >  1:2668977/1‑150 (MQ=255)
                                                               ccACGTCTTTGCGTTATGAACAACTCCTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCa                                                                                       <  2:3109544/150‑1 (MQ=255)
                                                               ccACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCa                                                                                       <  2:1123238/150‑1 (MQ=255)
                                                                   gTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGATAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTc                                                                                   >  2:2996414/1‑150 (MQ=255)
                                                                        tGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAg                                                                              >  2:1968761/1‑150 (MQ=255)
                                                                        tGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAg                                                                              >  2:1404049/1‑150 (MQ=255)
                                                                        tGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAg                                                                              <  2:2505150/150‑1 (MQ=255)
                                                                        tGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAg                                                                              >  2:271371/1‑150 (MQ=255)
                                                                          cGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGATAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGc                                                                            <  2:2206855/150‑1 (MQ=255)
                                                                            ttATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCgg                                                                          <  1:1068817/150‑1 (MQ=255)
                                                                            ttATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCgg                                                                          >  1:3068378/1‑150 (MQ=255)
                                                                              aTGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGaa                                                                        >  2:937923/1‑150 (MQ=255)
                                                                              aTGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGaa                                                                        >  2:616142/1‑150 (MQ=255)
                                                                               tGAACAACGACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCAACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAAt                                                                       <  2:4105077/150‑1 (MQ=255)
                                                                                gAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATg                                                                      <  2:1187223/150‑1 (MQ=255)
                                                                                   caacTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGAcc                                                                   >  1:3217468/1‑150 (MQ=255)
                                                                                     actactGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCaa                                                                 <  1:3114692/150‑1 (MQ=255)
                                                                                      ctactGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGATAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAAt                                                                >  1:4283634/1‑150 (MQ=255)
                                                                                      ctactGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAAt                                                                <  1:1244562/150‑1 (MQ=255)
                                                                                       tactGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGCCCATAGCGGTGCGCCACTGGATGAAAAACAGAAAAGGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATa                                                               <  1:1177850/150‑1 (MQ=255)
                                                                                       tactGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATa                                                               <  1:399895/150‑1 (MQ=255)
                                                                                       tacgGGCAGCGCGCAAGGAAGGGGCACCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATa                                                               <  1:3866370/150‑1 (MQ=255)
                                                                                        actGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATaa                                                              <  1:144616/150‑1 (MQ=255)
                                                                                          tGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAAt                                                            >  2:2570940/1‑150 (MQ=255)
                                                                                           ggCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGATAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAAtt                                                           >  1:589218/1‑150 (MQ=255)
                                                                                            gCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTa                                                          <  2:2912334/150‑1 (MQ=255)
                                                                                              aGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGATAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTAcc                                                        >  1:3392596/1‑150 (MQ=255)
                                                                                              aGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTAcc                                                        >  1:3591328/1‑150 (MQ=255)
                                                                                               gcgcgcAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTAccc                                                       <  2:4462485/150‑1 (MQ=255)
                                                                                                   gcAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGAc                                                   <  2:193499/150‑1 (MQ=255)
                                                                                                    cAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACg                                                  >  2:822119/1‑150 (MQ=255)
                                                                                                      aggaaggGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAg                                                <  2:3562229/150‑1 (MQ=255)
                                                                                                         aaggGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTg                                             <  2:1559850/150‑1 (MQ=255)
                                                                                                            gggCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCg                                          <  2:2371312/150‑1 (MQ=255)
                                                                                                             ggCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGATAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGa                                         >  2:3593728/1‑150 (MQ=255)
                                                                                                              gCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGAt                                        <  1:3364310/150‑1 (MQ=255)
                                                                                                               cAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATa                                       <  1:2035136/150‑1 (MQ=255)
                                                                                                                aaCCGTCCTCGACCATAGCGGTGCGCCACTGGATGATAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAg                                      <  2:3851321/150‑1 (MQ=255)
                                                                                                                aaCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGACTGACCAATAAATTACCCGGACGAGTTGGCGATAg                                      >  2:4438405/1‑150 (MQ=255)
                                                                                                                 aCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGt                                     <  2:1765388/150‑1 (MQ=255)
                                                                                                                 aCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGt                                     <  2:1603245/150‑1 (MQ=255)
                                                                                                                 aCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGt                                     <  2:4040059/150‑1 (MQ=255)
                                                                                                                     tCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCt                                 >  2:367074/1‑150 (MQ=255)
                                                                                                                       cTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTa                               >  2:2166461/1‑150 (MQ=255)
                                                                                                                        tCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAg                              <  2:3300916/150‑1 (MQ=255)
                                                                                                                          gACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTg                            >  2:2681815/1‑150 (MQ=255)
                                                                                                                          gACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTg                            >  1:3765703/1‑150 (MQ=255)
                                                                                                                             cATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGt                         <  1:2473715/150‑1 (MQ=255)
                                                                                                                              aTAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTg                        <  2:962730/150‑1 (MQ=255)
                                                                                                                              aTAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTg                        <  2:1725407/150‑1 (MQ=255)
                                                                                                                               tAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGc                       >  2:1901008/1‑150 (MQ=255)
                                                                                                                                aGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGcc                      <  1:362035/150‑1 (MQ=255)
                                                                                                                                       cgcCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCt               <  2:4502649/150‑1 (MQ=255)
                                                                                                                                       cgcCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCt               <  2:4502817/150‑1 (MQ=255)
                                                                                                                                       cgcCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCt               >  2:179424/1‑150 (MQ=255)
                                                                                                                                        gcCACTGGATGATAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTa              >  1:1722380/1‑150 (MQ=255)
                                                                                                                                        gcCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTa              >  1:3677518/1‑150 (MQ=255)
                                                                                                                                          cACTGGATGATAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACg            >  2:779294/1‑150 (MQ=255)
                                                                                                                                              ggATGATAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCaa        >  2:1514187/1‑150 (MQ=255)
                                                                                                                                               gATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAAt       <  2:4189315/150‑1 (MQ=255)
                                                                                                                                                aTGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATa      >  1:410248/1‑150 (MQ=255)
                                                                                                                                                aTGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATa      <  1:3164719/150‑1 (MQ=255)
                                                                                                                                                 tGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACGGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATaa     <  1:2372349/150‑1 (MQ=255)
                                                                                                                                                 tGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATaa     >  2:3951082/1‑150 (MQ=255)
                                                                                                                                                   aaaaaCAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACg   >  1:4412197/1‑150 (MQ=255)
                                                                                                                                                    aaaaCAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCCATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCGTAGTGGGTGCCGGATGCTACGCCAATAACGc  >  1:443329/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGC  >  CP009273/862394‑862691

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: ‑