Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,388,902 G→T 53.7% A558A (GCG→GCT mutL → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,388,9020GT53.7% 2.2 / inf 338A558A (GCG→GCTmutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (87/69);  new base T (100/82);  total (187/151)
Fisher's exact test for biased strand distribution p-value = 9.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.70e-01

Only 100 of 340 total aligned reads displayed.

TTGGGTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGG  >  CP009273/4388756‑4389050
                                                                                                                                                  |                                                                                                                                                    
      aTTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGAt                                                                                                                                              >  1:4525180/1‑149 (MQ=255)
       ttGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGAtt                                                                                                                                             <  1:3670013/149‑1 (MQ=255)
       ttGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTGACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGAtt                                                                                                                                             <  2:4806493/149‑1 (MQ=255)
         gATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGc                                                                                                                                           >  2:697543/1‑149 (MQ=255)
          aTTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCa                                                                                                                                          <  2:2622690/149‑1 (MQ=255)
              ccAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGaa                                                                                                                                      <  2:232703/149‑1 (MQ=255)
                  tCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCt                                                                                                                                  <  1:3324529/149‑1 (MQ=255)
                       tGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCtgatga                                                                                                                             <  1:3442822/149‑1 (MQ=255)
                         cacaGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGc                                                                                                                           <  2:4559933/149‑1 (MQ=255)
                            aGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGaa                                                                                                                        <  1:3247156/149‑1 (MQ=255)
                             gCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAAc                                                                                                                       <  2:3972209/149‑1 (MQ=255)
                                   gACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATgcgc                                                                                                                 >  1:1682713/1‑149 (MQ=255)
                                    aCCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCa                                                                                                                >  2:3689430/1‑149 (MQ=255)
                                     ccATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAg                                                                                                               <  1:1548839/149‑1 (MQ=255)
                                      cATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGt                                                                                                              >  2:2870385/1‑149 (MQ=255)
                                      cATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGt                                                                                                              >  2:268072/1‑149 (MQ=255)
                                       aTCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTg                                                                                                             <  1:2779743/149‑1 (MQ=255)
                                       aTCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTg                                                                                                             >  2:2795338/1‑149 (MQ=255)
                                          aGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTc                                                                                                          >  1:4162443/1‑149 (MQ=255)
                                          aGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTc                                                                                                          <  2:2541501/149‑1 (MQ=255)
                                            ggCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCaa                                                                                                        >  2:3503765/1‑149 (MQ=255)
                                             gCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAAt                                                                                                       <  2:2427695/149‑1 (MQ=255)
                                              cAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATg                                                                                                      >  1:3146594/1‑149 (MQ=255)
                                               aGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATgg                                                                                                     <  1:4645437/149‑1 (MQ=255)
                                                  gCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGcac                                                                                                  <  2:1196370/149‑1 (MQ=255)
                                                  gCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGcac                                                                                                  >  1:2492328/1‑149 (MQ=255)
                                                  gCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGcac                                                                                                  >  1:4373040/1‑149 (MQ=255)
                                                       tACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGcc                                                                                             <  2:2483731/149‑1 (MQ=255)
                                                        aCCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCa                                                                                            >  2:589756/1‑149 (MQ=255)
                                                        aCCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCa                                                                                            >  1:97/1‑149 (MQ=255)
                                                        aCCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCa                                                                                            >  1:62/1‑149 (MQ=255)
                                                          ccTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATa                                                                                          >  1:3772423/1‑149 (MQ=255)
                                                               cGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCct                                                                                     <  1:806142/149‑1 (MQ=255)
                                                                gCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctg                                                                                    >  1:1469490/1‑149 (MQ=255)
                                                                 ccAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgc                                                                                   >  2:4181174/1‑149 (MQ=255)
                                                                  caacaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgct                                                                                  <  1:4457783/149‑1 (MQ=255)
                                                                  caacaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgct                                                                                  >  1:902174/1‑149 (MQ=255)
                                                                  caacaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgct                                                                                  >  2:4908448/1‑149 (MQ=255)
                                                                  caacaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgct                                                                                  >  2:1789957/1‑149 (MQ=255)
                                                                     caaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGc                                                                               <  2:4701986/149‑1 (MQ=255)
                                                                                  tCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACgg                                                                  >  1:2457863/1‑149 (MQ=255)
                                                                                      gATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTAt                                                              <  1:4079163/149‑1 (MQ=255)
                                                                                       aTTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATg                                                             >  1:3017413/1‑149 (MQ=255)
                                                                                        ttCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGt                                                            >  2:2075079/1‑149 (MQ=255)
                                                                                         tCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTc                                                           >  2:5095121/1‑149 (MQ=255)
                                                                                          ccTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTcc                                                          <  2:4299065/149‑1 (MQ=255)
                                                                                            tGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGc                                                        >  1:2313327/1‑149 (MQ=255)
                                                                                            tGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACTGTTATGTCCGc                                                        >  1:2965023/1‑149 (MQ=255)
                                                                                            tGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGc                                                        >  1:2885721/1‑149 (MQ=255)
                                                                                              aaCTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCaa                                                      >  1:1800012/1‑149 (MQ=255)
                                                                                                cTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACt                                                    <  2:3625169/149‑1 (MQ=255)
                                                                                                 tGATAGGCTACCTGGCGTAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACtt                                                   >  2:4390946/1‑149 (MQ=255)
                                                                                                 tGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAGCATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACtt                                                   >  2:1799071/1‑149 (MQ=255)
                                                                                                 tGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACtt                                                   <  2:3918476/149‑1 (MQ=255)
                                                                                                  gATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTtg                                                  >  1:640315/1‑149 (MQ=255)
                                                                                                      ggCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaa                                              >  2:4976792/1‑149 (MQ=255)
                                                                                                      ggCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaa                                              >  1:2118682/1‑149 (MQ=255)
                                                                                                      ggCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaa                                              >  1:3826993/1‑149 (MQ=255)
                                                                                                       gCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaaa                                             <  2:1185324/149‑1 (MQ=255)
                                                                                                        cTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaaaa                                            >  2:4384946/1‑149 (MQ=255)
                                                                                                        cTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaaaa                                            <  1:1355775/149‑1 (MQ=255)
                                                                                                        cTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaaaa                                            >  1:3161183/1‑149 (MQ=255)
                                                                                                        cTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaaaa                                            >  1:1720389/1‑149 (MQ=255)
                                                                                                        cTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaaaa                                            >  1:561463/1‑149 (MQ=255)
                                                                                                          aCCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAATcg                                          >  2:2745602/1‑149 (MQ=255)
                                                                                                          aCCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcg                                          >  2:1285179/1‑149 (MQ=255)
                                                                                                           ccTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgc                                         <  2:831389/149‑1 (MQ=255)
                                                                                                            cTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgcc                                        >  2:2171545/1‑149 (MQ=255)
                                                                                                                 gAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCgg                                   >  1:1004055/1‑149 (MQ=255)
                                                                                                                 gAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCgg                                   >  1:3581093/1‑149 (MQ=255)
                                                                                                                   aGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggt                                 >  2:2647473/1‑149 (MQ=255)
                                                                                                                   aGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggt                                 >  2:4557152/1‑149 (MQ=255)
                                                                                                                   aGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggt                                 >  1:2637755/1‑149 (MQ=255)
                                                                                                                   aGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggt                                 >  2:1621932/1‑149 (MQ=255)
                                                                                                                   aGCAGGCCGTATTCGACCCTGGCAATATGGCGCAGTGGAGTGCACGAAATCTGAGGCGCGAACAGGCGCCGTGGTCAAGGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTAGGGCCGCAACTGGTGAAAACGCCGCCGggt                                 <  2:1324663/149‑1 (MQ=255)
                                                                                                                    gCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggtg                                >  1:4546941/1‑149 (MQ=255)
                                                                                                                    gCAGTCCGTATTCGAACCTGGAAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCACAACTTGTGAAAACGCCGCCGggtg                                <  2:4604548/149‑1 (MQ=255)
                                                                                                                     cAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggtgg                               <  2:1230710/149‑1 (MQ=255)
                                                                                                                      aGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggtggt                              >  1:4938543/1‑149 (MQ=255)
                                                                                                                          cGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGt                          <  2:138704/149‑1 (MQ=255)
                                                                                                                           gTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGtt                         <  1:4697635/149‑1 (MQ=255)
                                                                                                                           gTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGtt                         <  2:3805252/149‑1 (MQ=255)
                                                                                                                           gTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGtt                         <  2:3371799/149‑1 (MQ=255)
                                                                                                                           gTATTCGAACCGGGCAATAGTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGtt                         <  2:4645661/149‑1 (MQ=255)
                                                                                                                             aTTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTAc                       >  1:4041553/1‑149 (MQ=255)
                                                                                                                             aTTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTAc                       >  1:5171028/1‑149 (MQ=255)
                                                                                                                             aTTCGAACCTGCCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTAc                       <  2:2271957/149‑1 (MQ=255)
                                                                                                                              ttCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACa                      <  1:2855321/149‑1 (MQ=255)
                                                                                                                               tCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACaa                     <  1:2814990/149‑1 (MQ=255)
                                                                                                                                cGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAAt                    >  1:2560845/1‑149 (MQ=255)
                                                                                                                                cGAACCGGGCAATATTGCTCCGTGGATTGCACGAAATCTGATGAGCGAACAGGCGCAGTGGTCAATGGCACAGGCCATAACCCTGGGGGCGGACGTGGAACGGTTTTGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAAt                    <  1:5032771/149‑1 (MQ=255)
                                                                                                                                  aaCCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCt                  >  1:4207516/1‑149 (MQ=255)
                                                                                                                                   aCCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTg                 <  2:3094872/149‑1 (MQ=255)
                                                                                                                                   aCCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTTAAAACGCCGCCGGGTGGTCTGTTACAATCTg                 <  1:4976163/149‑1 (MQ=255)
                                                                                                                                       ggCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCAATAACACTGCTGGCTTACGTGGAACGTTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGa             <  2:156048/149‑1 (MQ=255)
                                                                                                                                        gCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAt            <  2:3132956/149‑1 (MQ=255)
                                                                                                                                            tatTGCTCTGTGTATTGCCCGAAATCTGATGTGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGGTGGCGGACGTGGAAGGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTAc        <  1:2376597/149‑1 (MQ=255)
                                                                                                                                             atTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACa       <  2:4158384/149‑1 (MQ=255)
                                                                                                                                               tGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATc     >  1:3508987/1‑149 (MQ=255)
                                                                                                                                                  gcAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCgg  >  2:911599/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTGGGTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGG  >  CP009273/4388756‑4389050

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: