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,8360CA97.3% 697.3 / 8.6 224A47A (GCC→GCAcsdAcysteine sulfinate desulfinase
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base A (96/122);  total (98/126)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 224 total aligned reads displayed.

GTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAAC  >  CP009273/2936691‑2936980
                                                                                                                                                 |                                                                                                                                                
   ccATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGaaa                                                                                                                                             >  2:1917482/1‑148 (MQ=255)
   ccATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGaaa                                                                                                                                             >  2:1021680/1‑148 (MQ=255)
   ccATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGaaa                                                                                                                                             <  2:3298725/148‑1 (MQ=255)
    cATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAAc                                                                                                                                            >  2:2761054/1‑148 (MQ=255)
    cATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAAc                                                                                                                                            >  1:152120/1‑148 (MQ=255)
     aTGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACg                                                                                                                                           <  1:1035223/148‑1 (MQ=255)
     aTGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACg                                                                                                                                           >  2:3513886/1‑148 (MQ=255)
      tGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGt                                                                                                                                          <  1:2354255/148‑1 (MQ=255)
      tGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGt                                                                                                                                          <  2:873931/148‑1 (MQ=255)
        aaCGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTcc                                                                                                                                        <  2:884982/148‑1 (MQ=255)
             ttttAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGc                                                                                                                                   >  2:1405781/1‑148 (MQ=255)
                tAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGc                                                                                                                                <  2:172084/148‑1 (MQ=255)
                 aaTCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGcc                                                                                                                               <  2:3021518/148‑1 (MQ=255)
                   tCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAg                                                                                                                             >  1:814302/1‑148 (MQ=255)
                    cccGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGt                                                                                                                            <  1:2156231/148‑1 (MQ=255)
                       gcgcAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTg                                                                                                                         <  1:3492977/148‑1 (MQ=255)
                       gcgcAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTg                                                                                                                         <  2:431464/148‑1 (MQ=255)
                        cgcAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGc                                                                                                                        >  1:3131716/1‑148 (MQ=255)
                          cAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCg                                                                                                                      >  2:1635111/1‑148 (MQ=255)
                           aGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGa                                                                                                                     <  1:2946688/148‑1 (MQ=255)
                            gTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGaa                                                                                                                    >  2:432451/1‑148 (MQ=255)
                            gTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGaa                                                                                                                    <  2:3091308/148‑1 (MQ=255)
                             tttCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAg                                                                                                                   >  1:3151156/1‑148 (MQ=255)
                              ttCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGc                                                                                                                  >  2:3054301/1‑148 (MQ=255)
                              ttCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGc                                                                                                                  >  2:2806083/1‑148 (MQ=255)
                               tCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGcc                                                                                                                 <  2:1284647/148‑1 (MQ=255)
                               tCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGcc                                                                                                                 <  2:546037/148‑1 (MQ=255)
                                 gcgcCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCa                                                                                                               <  1:451883/148‑1 (MQ=255)
                                     ccAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGc                                                                                                           <  2:775599/148‑1 (MQ=255)
                                     ccAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGc                                                                                                           <  2:2142241/148‑1 (MQ=255)
                                         tttCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCCACAGTTTTACAGTCTGAGCGCAGGAAACGTCCAACGCAGCCAGGTTGCCGAAGCCCAACGCCTGa                                                                                                       >  2:3201456/1‑148 (MQ=255)
                                         tttCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGa                                                                                                       >  2:1990269/1‑148 (MQ=255)
                                              cGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACcgcg                                                                                                  <  1:581200/148‑1 (MQ=255)
                                                 aCTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGt                                                                                               >  1:1266328/1‑148 (MQ=255)
                                                 aCTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGt                                                                                               >  1:2606984/1‑148 (MQ=255)
                                                  cTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGtt                                                                                              <  1:988105/148‑1 (MQ=255)
                                                   tACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTa                                                                                             >  1:161707/1‑148 (MQ=255)
                                                     cAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTTCCGAAGCCCAACGCCTGACCGCGCGTTATg                                                                                           >  1:1548145/1‑148 (MQ=255)
                                                     cAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATg                                                                                           <  2:3532921/148‑1 (MQ=255)
                                                     cAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATg                                                                                           <  2:400294/148‑1 (MQ=255)
                                                     cAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATg                                                                                           >  1:972739/1‑148 (MQ=255)
                                                      aGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGa                                                                                          <  1:1907704/148‑1 (MQ=255)
                                                       ggATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGaa                                                                                         >  1:1394105/1‑148 (MQ=255)
                                                        gATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAg                                                                                        >  2:845905/1‑148 (MQ=255)
                                                        gATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAg                                                                                        <  1:3262706/148‑1 (MQ=255)
                                                          tGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCt                                                                                      <  2:504053/148‑1 (MQ=255)
                                                             gggCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCa                                                                                   <  1:1977221/148‑1 (MQ=255)
                                                                      tctcGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGaaa                                                                          <  1:2971332/148‑1 (MQ=255)
                                                                       ctcGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAg                                                                         <  2:125408/148‑1 (MQ=255)
                                                                         cGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTg                                                                       <  1:739537/148‑1 (MQ=255)
                                                                          gACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTgg                                                                      <  1:1341465/148‑1 (MQ=255)
                                                                             aGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGgcgc                                                                   >  2:3000318/1‑148 (MQ=255)
                                                                               cgccgcGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCaa                                                                 >  2:574808/1‑148 (MQ=255)
                                                                                  cgcgACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTa                                                              <  1:2494860/148‑1 (MQ=255)
                                                                                  cgcgACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTa                                                              <  2:1483245/148‑1 (MQ=255)
                                                                                     gACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTg                                                           <  2:3056962/148‑1 (MQ=255)
                                                                                      aCCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGa                                                          >  1:230640/1‑148 (MQ=255)
                                                                                       ccGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCACAATTACTGaa                                                         >  2:3356258/1‑148 (MQ=255)
                                                                                        cgcgCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAAt                                                        <  1:1601934/148‑1 (MQ=255)
                                                                                        cgcgCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAAt                                                        <  1:761534/148‑1 (MQ=255)
                                                                                          cgcTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGc                                                      <  1:516226/148‑1 (MQ=255)
                                                                                               aaaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCgg                                                 <  1:2199100/148‑1 (MQ=255)
                                                                                               aaaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCgg                                                 >  2:2911738/1‑148 (MQ=255)
                                                                                               aaaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCgg                                                 <  1:2199124/148‑1 (MQ=255)
                                                                                               aaaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTCGCGCAATTACTGAATGCACCgg                                                 >  2:1139409/1‑148 (MQ=255)
                                                                                                aaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGga                                                >  1:2733170/1‑148 (MQ=255)
                                                                                                  ccTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGgatg                                              >  1:2081456/1‑148 (MQ=255)
                                                                                                   cTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGgatga                                             <  1:1174945/148‑1 (MQ=255)
                                                                                                   cTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGgatga                                             <  1:3363805/148‑1 (MQ=255)
                                                                                                    tGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGgatgat                                            >  1:3022467/1‑148 (MQ=255)
                                                                                                     gAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATa                                           <  2:2267040/148‑1 (MQ=255)
                                                                                                     gAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATa                                           <  2:2466300/148‑1 (MQ=255)
                                                                                                          cGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACt                                      <  1:2002192/148‑1 (MQ=255)
                                                                                                           gTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTa                                     <  1:1917482/148‑1 (MQ=255)
                                                                                                           gTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTa                                     <  2:2359573/148‑1 (MQ=255)
                                                                                                           gTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTa                                     <  1:3514910/148‑1 (MQ=255)
                                                                                                           gTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTa                                     <  1:2206774/148‑1 (MQ=255)
                                                                                                                tGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTc                                <  2:161707/148‑1 (MQ=255)
                                                                                                                   aGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACNCCTGACCGCGCGTTATGANGCTGCACGAGAGAAAGTGGCGCAATTACTGNATGCACCGGATGATAAAACTATCGTCTgg                             <  2:480043/148‑1 (MQ=255)
                                                                                                                   aGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGTATGATAAAACTATCGTCTgg                             <  1:432446/148‑1 (MQ=255)
                                                                                                                    gCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGa                            >  1:1995285/1‑148 (MQ=255)
                                                                                                                    gCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGa                            >  2:200402/1‑148 (MQ=255)
                                                                                                                    gCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGa                            >  1:649380/1‑148 (MQ=255)
                                                                                                                        cccAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGATAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGAcgcg                        >  1:1954884/1‑148 (MQ=255)
                                                                                                                        cccAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGAcgcg                        >  2:415036/1‑148 (MQ=255)
                                                                                                                          cAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGAcgcgcg                      >  1:3224362/1‑148 (MQ=255)
                                                                                                                             cAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGca                   <  2:2670498/148‑1 (MQ=255)
                                                                                                                              aGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGcac                  >  2:1074757/1‑148 (MQ=255)
                                                                                                                                ttttACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGcacca                >  2:25886/1‑148 (MQ=255)
                                                                                                                                ttttACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGcacca                >  2:3239258/1‑148 (MQ=255)
                                                                                                                                 tttACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGcaccac               >  2:322246/1‑148 (MQ=255)
                                                                                                                                    aCAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGa            <  1:2792311/148‑1 (MQ=255)
                                                                                                                                     cAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCATCGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGaa           <  2:1202194/148‑1 (MQ=255)
                                                                                                                                     cAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGaa           >  2:117028/1‑148 (MQ=255)
                                                                                                                                        tCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATcc        >  2:1400057/1‑148 (MQ=255)
                                                                                                                                         cttaGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCa       >  2:628833/4‑148 (MQ=255)
                                                                                                                                          tGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCAt      <  2:1505327/148‑1 (MQ=255)
                                                                                                                                            aGCGCAGGAAACGTCCATCGCAGCCAGTGTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCa    <  1:3431342/148‑1 (MQ=255)
                                                                                                                                              cgcAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAAc  <  2:3156418/148‑1 (MQ=255)
                                                                                                                                              cgcAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAAc  <  1:2021182/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAAC  >  CP009273/2936691‑2936980

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: