Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,873,087 G→T 100% V15V (GTC→GTA ygbT ← multifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,873,0870GT94.6% 949.4 / 37.0 280V15V (GTC→GTAygbTmultifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme
Reads supporting (aligned to +/- strand):  ref base G (1/14);  new base T (134/131);  total (135/145)
Fisher's exact test for biased strand distribution p-value = 8.08e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 280 total aligned reads displayed.

AGCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGC  >  CP009273/2872940‑2873233
                                                                                                                                                   |                                                                                                                                                  
aGCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAc                                                                                                                                                  <  2:1348968/150‑1 (MQ=255)
 gCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAcg                                                                                                                                                 <  1:360569/150‑1 (MQ=255)
   tGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAcgcg                                                                                                                                               >  1:3067008/1‑150 (MQ=255)
   tGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATAcgcg                                                                                                                                               >  2:2636289/1‑150 (MQ=255)
     cATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGAt                                                                                                                                             <  1:3741048/150‑1 (MQ=255)
     cATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGAt                                                                                                                                             <  1:1349629/150‑1 (MQ=255)
        gCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCtt                                                                                                                                          <  1:848891/150‑1 (MQ=255)
         cGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCttt                                                                                                                                         <  1:2039057/150‑1 (MQ=255)
               ccGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTgg                                                                                                                                   >  2:4712537/1‑150 (MQ=255)
                 gtgtACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGGGGaa                                                                                                                                 >  2:1187704/1‑150 (MQ=255)
                     aCCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATggg                                                                                                                             >  1:537238/1‑150 (MQ=255)
                      ccAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGa                                                                                                                            <  2:941491/150‑1 (MQ=255)
                        aGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGAtt                                                                                                                          >  1:3101070/1‑150 (MQ=255)
                          gTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTaa                                                                                                                        >  1:4490239/1‑150 (MQ=255)
                            tCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAgg                                                                                                                      >  1:1689878/1‑150 (MQ=255)
                            tCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAgg                                                                                                                      <  2:130816/150‑1 (MQ=255)
                             ccAGCATGATGCAGGCAACCGAGCCAACAGGAATCTGAGTGCGGAGCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAggg                                                                                                                     <  2:2467577/150‑1 (MQ=255)
                             ccAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAggg                                                                                                                     <  2:2055825/150‑1 (MQ=255)
                              cAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAgggg                                                                                                                    >  2:3656453/1‑150 (MQ=255)
                                 cATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAg                                                                                                                 <  2:721705/150‑1 (MQ=255)
                                  atgatgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGc                                                                                                                >  1:306594/1‑150 (MQ=255)
                                  atgatgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGc                                                                                                                <  1:3631911/150‑1 (MQ=255)
                                     atgCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATCTTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAg                                                                                                             <  2:318453/150‑1 (MQ=255)
                                          ggCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAt                                                                                                        <  2:1310733/150‑1 (MQ=255)
                                          ggCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAt                                                                                                        <  2:3450503/150‑1 (MQ=255)
                                             aaCCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAtttt                                                                                                     <  2:679898/150‑1 (MQ=255)
                                             aaCCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCAtttt                                                                                                     <  1:484491/150‑1 (MQ=255)
                                               ccGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTttat                                                                                                   <  2:595074/150‑1 (MQ=255)
                                                 gAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTttatta                                                                                                 >  1:1740430/1‑150 (MQ=255)
                                                 gAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTttatta                                                                                                 >  1:744690/1‑150 (MQ=255)
                                                  aGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTac                                                                                                <  1:2215732/150‑1 (MQ=255)
                                                    ccAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATAGTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTacac                                                                                              <  2:2165756/150‑1 (MQ=255)
                                                     cAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATTGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACAcc                                                                                             <  1:619384/150‑1 (MQ=255)
                                                        cAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCa                                                                                          <  1:2937043/150‑1 (MQ=255)
                                                           gAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATc                                                                                       <  2:1794761/150‑1 (MQ=255)
                                                              tatGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATcaca                                                                                    <  1:4712537/150‑1 (MQ=255)
                                                               atGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAg                                                                                   >  1:949169/1‑150 (MQ=255)
                                                                  aGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTgg                                                                                >  2:2836492/1‑150 (MQ=255)
                                                                  aGTGCGGATCCCTGTCTTGACGATAAGTACAAACGCGCCATCTATTACATCGATCTGCACATATAGCAGAAAGATCATGGATACGCGATATTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATATTATTACACCTCAATCACAGTgg                                                                                >  2:4163423/1‑150 (MQ=255)
                                                                   gTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGa                                                                               >  1:2342154/1‑150 (MQ=255)
                                                                         aTCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCaaa                                                                         >  2:3480982/1‑150 (MQ=255)
                                                                         aTCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCaaa                                                                         >  2:1005311/1‑150 (MQ=255)
                                                                          tCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAg                                                                        >  1:3953035/1‑150 (MQ=255)
                                                                             cTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATa                                                                     <  1:1531026/150‑1 (MQ=255)
                                                                                     tCGAGAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATAGTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCca                                                             <  2:4438236/150‑1 (MQ=255)
                                                                                      cGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCcac                                                            <  1:2898865/150‑1 (MQ=255)
                                                                                        aTAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACAt                                                          >  1:3231391/1‑150 (MQ=255)
                                                                                        aTAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACAt                                                          <  2:3944923/150‑1 (MQ=255)
                                                                                         tAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATc                                                         >  1:2963983/1‑150 (MQ=255)
                                                                                         tAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATc                                                         >  1:3428173/1‑150 (MQ=255)
                                                                                         tAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATc                                                         >  2:4049053/1‑150 (MQ=255)
                                                                                          aaGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATcc                                                        <  1:4048396/150‑1 (MQ=255)
                                                                                          aaGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATcc                                                        <  1:4223254/150‑1 (MQ=255)
                                                                                           aGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATccc                                                       >  1:3592190/1‑150 (MQ=255)
                                                                                           aGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATccc                                                       <  1:2816011/150‑1 (MQ=255)
                                                                                            gTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCa                                                      <  2:3881948/150‑1 (MQ=255)
                                                                                             tACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCAt                                                     >  1:1453848/1‑150 (MQ=255)
                                                                                               cAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCg                                                   >  2:3313756/1‑150 (MQ=255)
                                                                                                aaaCGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGa                                                  >  1:2232825/1‑150 (MQ=255)
                                                                                                 aaCGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGAt                                                 >  2:4542179/1‑150 (MQ=255)
                                                                                                 aaCGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGAt                                                 <  1:4110890/150‑1 (MQ=255)
                                                                                                 aaCGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGAt                                                 <  1:2793993/150‑1 (MQ=255)
                                                                                                    gcgcCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTa                                              <  1:2568837/150‑1 (MQ=255)
                                                                                                         aTCTATTACATCGATCTGCCCATATTGCAGAACGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTgg                                         <  1:2622731/150‑1 (MQ=255)
                                                                                                             aTTACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCa                                     <  2:1585386/150‑1 (MQ=255)
                                                                                                              ttACATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCaa                                    >  2:2751690/1‑150 (MQ=255)
                                                                                                                aCATCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATa                                  >  1:1872086/1‑150 (MQ=255)
                                                                                                                  aTCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATAcc                                <  2:530196/150‑1 (MQ=255)
                                                                                                                   tCGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCACTCACAGTGGAGCCAACGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCt                               <  2:3101070/150‑1 (MQ=255)
                                                                                                                    cGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCtt                              <  2:4002721/150‑1 (MQ=255)
                                                                                                                    cGATCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCtt                              <  2:2782056/150‑1 (MQ=255)
                                                                                                                       tCTGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCt                           <  2:3551201/150‑1 (MQ=255)
                                                                                                                         tGCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGt                         >  1:4478351/1‑150 (MQ=255)
                                                                                                                          gCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTa                        <  1:113697/150‑1 (MQ=255)
                                                                                                                          gCCCATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTa                        <  1:1993851/150‑1 (MQ=255)
                                                                                                                           cccATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTAc                       <  1:4303401/150‑1 (MQ=255)
                                                                                                                           cccATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTAc                       >  1:460957/1‑150 (MQ=255)
                                                                                                                           cccATATTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAACAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTAc                       <  2:3351826/150‑1 (MQ=255)
                                                                                                                               tatTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGa                   >  2:2539101/1‑150 (MQ=255)
                                                                                                                                atTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGAt                  >  1:1034060/1‑150 (MQ=255)
                                                                                                                                atTGCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGAt                  >  1:2068090/1‑150 (MQ=255)
                                                                                                                                 ttGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATc                 <  2:3198950/150‑1 (MQ=255)
                                                                                                                                   gCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTa               <  2:717946/150‑1 (MQ=255)
                                                                                                                                   gCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTa               <  1:214063/150‑1 (MQ=255)
                                                                                                                                   gCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTa               <  2:4659752/150‑1 (MQ=255)
                                                                                                                                   gCAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTa               <  2:1118950/150‑1 (MQ=255)
                                                                                                                                    cAGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTAt              <  1:510543/150‑1 (MQ=255)
                                                                                                                                     aGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTAtt             >  1:3675072/1‑150 (MQ=255)
                                                                                                                                     aGAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTAtt             <  1:1227756/150‑1 (MQ=255)
                                                                                                                                      gAAAGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTa            >  2:1920613/1‑150 (MQ=255)
                                                                                                                                       aaaGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTaa           <  2:3769895/150‑1 (MQ=255)
                                                                                                                                       aaaGATCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTaa           >  2:1314986/1‑150 (MQ=255)
                                                                                                                                           aTCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCt       <  2:227275/150‑1 (MQ=255)
                                                                                                                                           aTCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCt       <  2:1486568/150‑1 (MQ=255)
                                                                                                                                           aTCATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCt       >  2:2364005/1‑150 (MQ=255)
                                                                                                                                             cATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTgg     >  1:1605764/1‑150 (MQ=255)
                                                                                                                                             cATGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTgg     <  2:1826696/150‑1 (MQ=255)
                                                                                                                                              aTGGATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGAGAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGgc    <  1:242672/150‑1 (MQ=255)
                                                                                                                                                ggATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGgcgc  >  2:2342552/1‑150 (MQ=255)
                                                                                                                                                ggATACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGgcgc  >  2:907975/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGCTGCATGCGAAACCCGTGTACCAGGTTCCAGCATGATGCAGGCAACCGAGCCAACAGGAATATGAGTGCGGATCCCTGTCTTGTCGATAAGTACAAACGCGCCATCTATTACATCGATCTGCCCATATTGCAGAAAGATCATGGAGACGCGATCTTTGAGTGGAATGGGATTAAGGGGAAGCCAGGTCATTTTATTACACCTCAATCACAGTGGAGCCAAAGATAGCAAGCCACATCCCATCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGC  >  CP009273/2872940‑2873233

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: