Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,122,781 C→A 31.5% R80S (CGC→AGC)  metF → 5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,122,7810CA31.5% 523.8 / 228.4 269R80S (CGC→AGC) metF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (82/102);  new base A (44/41);  total (126/143)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 275 total aligned reads displayed.

TCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGAAGTGGTAAGCCAGAA  >  CP009273/4122632‑4122930
                                                                                                                                                     |                                                                                                                                                     
tCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATc                                                                                                                                                       >  1:1684727/1‑150 (MQ=255)
tCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATc                                                                                                                                                       >  1:2109288/1‑150 (MQ=255)
  ccgccgCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGc                                                                                                                                                     <  2:3713932/150‑1 (MQ=255)
  ccgccgCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGc                                                                                                                                                     <  1:777829/150‑1 (MQ=255)
     ccgcGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACt                                                                                                                                                  <  2:2887726/150‑1 (MQ=255)
     ccgcGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACt                                                                                                                                                  <  2:2887738/150‑1 (MQ=255)
        cgTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGt                                                                                                                                               <  1:2814614/150‑1 (MQ=255)
        cgTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGt                                                                                                                                               <  1:2248606/150‑1 (MQ=255)
         gTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTc                                                                                                                                              <  2:3702470/150‑1 (MQ=255)
           aCCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTg                                                                                                                                            >  1:4147269/1‑150 (MQ=255)
           aCCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTg                                                                                                                                            <  1:7002/150‑1 (MQ=255)
            ccAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTgg                                                                                                                                           >  2:3150039/1‑150 (MQ=255)
             cAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGa                                                                                                                                          <  1:192956/150‑1 (MQ=255)
             cAGAGAAATGGAGCAGACCCTGTGTAACTCCATCGATCGCCTTAGCAGCATGAAACCGAAGATTGAATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACAGTACGCACAGCATTATTAAACGCATTAAAGATCGCACTGGTCAGGa                                                                                                                                          >  2:2825242/1‑150 (MQ=255)
              aGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGaa                                                                                                                                         >  1:451918/1‑150 (MQ=255)
               gTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGacg                                                                                                                                        >  2:1573477/1‑148 (MQ=255)
               gTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAg                                                                                                                                        >  2:2413879/1‑150 (MQ=255)
                 gAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCg                                                                                                                                      >  2:323933/1‑150 (MQ=255)
                 gAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCg                                                                                                                                      <  2:3605112/150‑1 (MQ=255)
                 gAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCg                                                                                                                                      <  2:3850996/150‑1 (MQ=255)
                   aaTGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGc                                                                                                                                    <  1:1230978/150‑1 (MQ=255)
                      tgAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCAcc                                                                                                                                 <  2:1160074/149‑1 (MQ=255)
                      ggAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCAcc                                                                                                                                 <  2:2558158/150‑1 (MQ=255)
                      ggAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCAcc                                                                                                                                 >  2:298000/1‑150 (MQ=255)
                        aGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGc                                                                                                                               >  2:230917/1‑150 (MQ=255)
                         gCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCa                                                                                                                              <  2:2102941/150‑1 (MQ=255)
                              cccTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTa                                                                                                                         >  1:1633102/1‑150 (MQ=255)
                               ccTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTAc                                                                                                                        >  2:13905/1‑150 (MQ=255)
                                cTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACt                                                                                                                       <  1:276605/150‑1 (MQ=255)
                                cTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACt                                                                                                                       <  1:4415473/150‑1 (MQ=255)
                                 tgtgGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACtt                                                                                                                      <  1:1740611/150‑1 (MQ=255)
                                     gAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCAt                                                                                                                  <  2:3302763/150‑1 (MQ=255)
                                       aCTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTg                                                                                                                <  1:1028615/150‑1 (MQ=255)
                                         tCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGAt                                                                                                              >  1:2161984/1‑150 (MQ=255)
                                            atcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCg                                                                                                           >  1:3843679/1‑150 (MQ=255)
                                               gatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACg                                                                                                        >  2:928518/1‑150 (MQ=255)
                                                 tcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGcc                                                                                                      >  2:645581/1‑150 (MQ=255)
                                                 tcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGcc                                                                                                      <  2:1214349/150‑1 (MQ=255)
                                                       tAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACATGCATTGATGCGACGCCcgacga                                                                                                >  1:1572148/1‑150 (MQ=255)
                                                        agcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAg                                                                                               >  1:3635654/1‑150 (MQ=255)
                                                        agcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAg                                                                                               >  1:4096363/1‑150 (MQ=255)
                                                        agcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAg                                                                                               >  2:732306/1‑150 (MQ=255)
                                                        agcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAg                                                                                               >  1:39453/1‑150 (MQ=255)
                                                        agcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAg                                                                                               >  2:1726340/1‑150 (MQ=255)
                                                        agcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAg                                                                                               >  2:192726/1‑150 (MQ=255)
                                                          cagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCt                                                                                             <  2:2771476/150‑1 (MQ=255)
                                                          cagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCt                                                                                             <  1:4377421/150‑1 (MQ=255)
                                                           agcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTg                                                                                            >  1:3156195/1‑150 (MQ=255)
                                                           agcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTg                                                                                            >  2:2681034/1‑150 (MQ=255)
                                                           agcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTg                                                                                            >  2:4568457/1‑150 (MQ=255)
                                                               tGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCa                                                                                        <  1:1417085/150‑1 (MQ=255)
                                                               tGAAACCGAAGTATGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCa                                                                                        >  1:4576792/1‑150 (MQ=255)
                                                                gAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCAc                                                                                       >  1:1301620/1‑150 (MQ=255)
                                                                  aaCCGAAGTTTGTATCGGTGACCTATGGCGCGAACACCGGCGAGCGCGACCGTACGCACAGCATTAATAAAGGCATTAAATATCGCACTGGTCTGGAAGCGTCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCa                                                                                     >  2:756298/1‑150 (MQ=255)
                                                                   aCCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCAt                                                                                    <  1:67686/150‑1 (MQ=255)
                                                                   aCCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCAt                                                                                    <  1:99324/150‑1 (MQ=255)
                                                                      gAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGc                                                                                 >  1:1964046/1‑150 (MQ=255)
                                                                        aGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCAc                                                                               <  1:624688/150‑1 (MQ=255)
                                                                           ttGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCAcgcg                                                                            <  2:4454368/150‑1 (MQ=255)
                                                                            tGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGa                                                                           >  2:2177690/1‑150 (MQ=255)
                                                                             gTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGac                                                                          <  2:2363073/150‑1 (MQ=255)
                                                                             gTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGac                                                                          <  2:3572552/150‑1 (MQ=255)
                                                                                     gACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGaa                                                                  <  2:2013921/150‑1 (MQ=255)
                                                                                        cTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATaa                                                               <  1:1424780/150‑1 (MQ=255)
                                                                                         tATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAAc                                                              <  2:2563682/150‑1 (MQ=255)
                                                                                              cgcgAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTAt                                                         <  1:489919/150‑1 (MQ=255)
                                                                                                  aaCTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGt                                                     <  1:1836851/150‑1 (MQ=255)
                                                                                                     tCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCat                                                  <  2:4017596/150‑1 (MQ=255)
                                                                                                      ccGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCata                                                 <  2:4396655/150‑1 (MQ=255)
                                                                                                        ggCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTCTTCGTCATATc                                               >  1:3221105/1‑150 (MQ=255)
                                                                                                        ggCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATc                                               <  1:2796557/150‑1 (MQ=255)
                                                                                                        ggCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATc                                               <  1:1368203/150‑1 (MQ=255)
                                                                                                          cGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCAGATCGt                                             >  2:4219843/1‑150 (MQ=255)
                                                                                                           gAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTg                                            >  2:4327098/1‑150 (MQ=255)
                                                                                                            aGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTgg                                           >  1:87508/1‑150 (MQ=255)
                                                                                                             gcgcGACCGTACGCAAAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGgc                                          <  1:3458971/150‑1 (MQ=255)
                                                                                                                cgACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCt                                       >  1:3021705/1‑150 (MQ=255)
                                                                                                                cgACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCt                                       >  2:3919228/1‑150 (MQ=255)
                                                                                                                 gACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTg                                      <  2:2071671/150‑1 (MQ=255)
                                                                                                                 gACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTg                                      >  1:2895133/1‑150 (MQ=255)
                                                                                                                           acaGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATc                            <  1:1206864/150‑1 (MQ=255)
                                                                                                                            caGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCt                           >  2:867570/1‑150 (MQ=255)
                                                                                                                              gCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTgc                         <  2:1189100/150‑1 (MQ=255)
                                                                                                                               cATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTgcc                        <  2:3066114/150‑1 (MQ=255)
                                                                                                                                attattAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTgccg                       <  2:1633102/150‑1 (MQ=255)
                                                                                                                                attattAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTgccg                       >  2:1418463/1‑150 (MQ=255)
                                                                                                                                attattAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTgccg                       <  1:4296013/150‑1 (MQ=255)
                                                                                                                                 ttattaAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTgccgc                      >  1:49616/1‑150 (MQ=255)
                                                                                                                                   attaAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGccgccg                    >  1:236671/1‑150 (MQ=255)
                                                                                                                                    ttaAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCATACTTGCATTTATGCTACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGACAAATCGTGGAGCTGCGTGGAGATCTGCCGCCgg                   >  1:2624372/1‑150 (MQ=255)
                                                                                                                                     taAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCggg                  >  2:4295688/1‑150 (MQ=255)
                                                                                                                                     taAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCggg                  >  2:1044475/1‑150 (MQ=255)
                                                                                                                                       aaGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGaa                <  2:763300/150‑1 (MQ=255)
                                                                                                                                           cATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGAAGTgg            <  2:1826574/150‑1 (MQ=255)
                                                                                                                                            aTTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGAAGTGGt           >  2:4123329/1‑150 (MQ=255)
                                                                                                                                             ttAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGAAGTGGTa          <  1:2433481/150‑1 (MQ=255)
                                                                                                                                             ttAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGAAGTGGTa          >  2:2707301/1‑150 (MQ=255)
                                                                                                                                                aaGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGAAGTGGTAAGc       <  1:3463914/150‑1 (MQ=255)
                                                                                                                                                  gATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGAAGTGGTAAGCCa     >  1:2929926/1‑150 (MQ=255)
                                                                                                                                                  gATCGCACTGGTCAGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCAAAGCACGCGACTACGGGAATAACGGTATACGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGAAGTGGTAAGCCa     >  2:2994469/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGAAGTGGTAAGCCAGAA  >  CP009273/4122632‑4122930

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: