Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,388,902 G→T 87.6% 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,9020GT87.6% 973.9 / 135.3 347A558A (GCG→GCT) mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (23/20);  new base T (185/119);  total (208/139)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 348 total aligned reads displayed.

TTGGGTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGG  >  CP009273/4388756‑4389050
                                                                                                                                                  |                                                                                                                                                    
 tGGGTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAg                                                                                                                                                   <  1:1831863/149‑1 (MQ=255)
   ggTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTg                                                                                                                                                 <  1:2221011/149‑1 (MQ=255)
      aTTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGAt                                                                                                                                              >  2:2439526/1‑149 (MQ=255)
       ttGATTTCCAGTCAGATACACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGAtt                                                                                                                                             <  1:174994/149‑1 (MQ=255)
        tGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTg                                                                                                                                            >  2:2967086/1‑149 (MQ=255)
         gATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGc                                                                                                                                           <  2:3329989/149‑1 (MQ=255)
         gATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGc                                                                                                                                           <  2:3290376/149‑1 (MQ=255)
          aTTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCa                                                                                                                                          >  1:1491322/1‑149 (MQ=255)
                 gTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATc                                                                                                                                   >  1:109596/1‑149 (MQ=255)
                 gTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATc                                                                                                                                   >  2:1795706/1‑149 (MQ=255)
                 gTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATc                                                                                                                                   >  2:2400783/1‑149 (MQ=255)
                  tCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCt                                                                                                                                  <  2:513690/149‑1 (MQ=255)
                   cAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCtg                                                                                                                                 >  2:4247303/1‑149 (MQ=255)
                     gATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCtgat                                                                                                                               >  1:3331582/1‑149 (MQ=255)
                      aTGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCtgatg                                                                                                                              <  1:1888871/149‑1 (MQ=255)
                      aTGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCtgatg                                                                                                                              >  1:3509359/1‑149 (MQ=255)
                        gCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAg                                                                                                                            >  1:1155445/1‑149 (MQ=255)
                            aGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGaa                                                                                                                        <  1:4206429/149‑1 (MQ=255)
                              cATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACa                                                                                                                      <  1:2883564/149‑1 (MQ=255)
                              cATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACa                                                                                                                      <  2:454634/149‑1 (MQ=255)
                                   gACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATgcgc                                                                                                                 <  2:2125408/149‑1 (MQ=255)
                                            ggCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCaa                                                                                                        <  2:3603342/149‑1 (MQ=255)
                                            ggCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCaa                                                                                                        <  1:971310/149‑1 (MQ=255)
                                                gTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGc                                                                                                    >  1:2200725/1‑149 (MQ=255)
                                                 tGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGca                                                                                                   <  2:2256490/149‑1 (MQ=255)
                                                   ccTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGcaca                                                                                                 >  2:2809205/1‑149 (MQ=255)
                                                    cTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAg                                                                                                >  2:1947392/1‑149 (MQ=255)
                                                     tttACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCCCAgg                                                                                               >  2:830112/1‑149 (MQ=255)
                                                     tttACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAgg                                                                                               >  2:4387304/1‑149 (MQ=255)
                                                      ttACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGc                                                                                              >  1:426100/1‑149 (MQ=255)
                                                      ttACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGc                                                                                              <  1:1539560/149‑1 (MQ=255)
                                                       tACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGcc                                                                                             <  2:2314106/149‑1 (MQ=255)
                                                        aCCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCa                                                                                            >  2:3420541/1‑149 (MQ=255)
                                                         cccTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCAt                                                                                           <  2:1339824/149‑1 (MQ=255)
                                                         cccTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCAt                                                                                           >  1:2151077/1‑149 (MQ=255)
                                                           cTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATaa                                                                                         <  1:194772/149‑1 (MQ=255)
                                                              aCGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAAccc                                                                                      >  1:3770277/1‑149 (MQ=255)
                                                                 ccAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgc                                                                                   >  1:63674/1‑149 (MQ=255)
                                                                  caacaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgct                                                                                  >  2:1539750/1‑149 (MQ=255)
                                                                  caacaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgct                                                                                  >  1:2629862/1‑149 (MQ=255)
                                                                   aacaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgctg                                                                                 >  2:2156047/1‑149 (MQ=255)
                                                                   aacaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgctg                                                                                 >  1:2476517/1‑149 (MQ=255)
                                                                     gaataTTGACAAATGTTGATTGCTTACCTGATTGGCTAGCTGGCGTAGCAGTACGTATTCGAACCTGTCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGAGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGc                                                                               <  1:4640551/145‑1 (MQ=255)
                                                                      aaaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATAGTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCg                                                                              <  1:3798244/149‑1 (MQ=255)
                                                                       aaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCgg                                                                             >  1:3131849/1‑149 (MQ=255)
                                                                         aTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGAc                                                                           <  2:1254164/149‑1 (MQ=255)
                                                                            tACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTg                                                                        <  2:2380540/149‑1 (MQ=255)
                                                                            tACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTg                                                                        >  1:790442/1‑149 (MQ=255)
                                                                            tACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTg                                                                        >  1:319593/1‑149 (MQ=255)
                                                                             aCAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTgg                                                                       <  1:1184941/149‑1 (MQ=255)
                                                                               aaaTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGaa                                                                     >  1:2427885/1‑149 (MQ=255)
                                                                               aaaTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGaa                                                                     <  2:1157130/149‑1 (MQ=255)
                                                                                aaTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAAc                                                                    <  1:4267430/149‑1 (MQ=255)
                                                                                aaTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAAc                                                                    <  1:2745559/149‑1 (MQ=255)
                                                                                 aTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACg                                                                   >  1:3986500/1‑149 (MQ=255)
                                                                                  tCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACgg                                                                  >  1:4515756/1‑149 (MQ=255)
                                                                                  tCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACgg                                                                  >  1:1131126/1‑149 (MQ=255)
                                                                                     tGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTa                                                               >  1:3697923/1‑149 (MQ=255)
                                                                                       aTTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATg                                                             >  1:3583522/1‑149 (MQ=255)
                                                                                         tCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTc                                                           <  1:1464365/149‑1 (MQ=255)
                                                                                             gAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCa                                                       >  2:2526741/1‑149 (MQ=255)
                                                                                             gAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCa                                                       >  2:2294583/1‑149 (MQ=255)
                                                                                             gAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCa                                                       >  2:679702/1‑149 (MQ=255)
                                                                                                   aTAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTtgt                                                 >  1:851338/1‑149 (MQ=255)
                                                                                                      ggCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaa                                              <  1:4173872/149‑1 (MQ=255)
                                                                                                       gCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaaa                                             <  2:3120953/149‑1 (MQ=255)
                                                                                                        cTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaaaa                                            <  1:2828820/149‑1 (MQ=255)
                                                                                                         tACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAc                                           >  2:3747278/1‑149 (MQ=255)
                                                                                                             tGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgccg                                       >  1:1066368/1‑149 (MQ=255)
                                                                                                              ggCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgccgc                                      >  2:3690439/1‑149 (MQ=255)
                                                                                                              ggCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgccgc                                      >  2:3822950/1‑149 (MQ=255)
                                                                                                                 gAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCgg                                   <  1:4571330/149‑1 (MQ=255)
                                                                                                                  aaGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCggg                                  >  2:2779333/1‑149 (MQ=255)
                                                                                                                  aaGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCggg                                  >  1:2725298/1‑149 (MQ=255)
                                                                                                                      aGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggtggt                              >  2:1414824/1‑149 (MQ=255)
                                                                                                                      aGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggtggt                              >  2:4673578/1‑149 (MQ=255)
                                                                                                                       gTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTc                             >  2:65549/1‑149 (MQ=255)
                                                                                                                       gTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTc                             >  2:235723/1‑149 (MQ=255)
                                                                                                                       gTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTc                             <  1:2295257/149‑1 (MQ=255)
                                                                                                                          cGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGt                          <  2:1603749/149‑1 (MQ=255)
                                                                                                                            tATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTa                        >  2:1657620/1‑149 (MQ=255)
                                                                                                                             aTTCGAACCTGGCAATAGTGCGCAGGGGATGGCACGAAATCTGATGAGCGAACATGCGCAGGGGTCAATGGCACAGGCCATACCCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTAc                       <  1:1014436/149‑1 (MQ=255)
                                                                                                                              ttCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACa                      >  1:4371639/1‑149 (MQ=255)
                                                                                                                               tCGACCCTGGCAATATTGCTCTGTGGTTTGCACGAAATCTGTTGAGCGAACAGGCGCTGTGGTCAATGGCACAGGCCATACCCCTGCGGGCGGACGTGGAACGGTTATGGCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACaa                     <  1:3960014/149‑1 (MQ=255)
                                                                                                                               tCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACaa                     <  2:1716775/149‑1 (MQ=255)
                                                                                                                                cGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGACCGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAAt                    <  1:1630619/149‑1 (MQ=255)
                                                                                                                                cGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAAt                    >  1:1975434/1‑149 (MQ=255)
                                                                                                                                cGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAAt                    >  1:310786/1‑149 (MQ=255)
                                                                                                                                cGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAAt                    >  1:527989/1‑149 (MQ=255)
                                                                                                                                 gAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATc                   <  1:728257/149‑1 (MQ=255)
                                                                                                                                    ccTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGt                <  2:3060396/149‑1 (MQ=255)
                                                                                                                                      tGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTg              >  2:1097140/1‑149 (MQ=255)
                                                                                                                                      tGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTg              <  1:1146350/149‑1 (MQ=255)
                                                                                                                                      tGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTg              >  2:898708/1‑149 (MQ=255)
                                                                                                                                       ggCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGa             >  2:4617613/1‑149 (MQ=255)
                                                                                                                                        gCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAt            >  1:1317784/1‑149 (MQ=255)
                                                                                                                                           atatTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTa         >  2:3211564/1‑149 (MQ=255)
                                                                                                                                           atatTGCGCTGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACTGGCCATAACCCTGGTGGCGGTGGTGGAACGGTTTTGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTa         <  1:3822950/149‑1 (MQ=255)
                                                                                                                                               tGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATc     <  2:4185657/149‑1 (MQ=255)
                                                                                                                                                  gcAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCgg  >  2:3975338/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: ‑