Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 372,581 G→T 100% A118A (GCG→GCT yaiL → DUF2058 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273372,5810GT97.5% 1026.4 / 11.4 322A118A (GCG→GCTyaiLDUF2058 family protein
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base T (168/146);  total (171/151)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 322 total aligned reads displayed.

AAAGCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCC  >  CP009273/372435‑372727
                                                                                                                                                  |                                                                                                                                                  
aaaGCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTAtt                                                                                                                                                 >  2:487717/1‑150 (MQ=255)
aaaGCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTAtt                                                                                                                                                 >  2:3568191/1‑150 (MQ=255)
 aaGCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTg                                                                                                                                                >  2:5004801/1‑150 (MQ=255)
  aGCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGc                                                                                                                                               >  2:3994473/1‑150 (MQ=255)
     tCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCg                                                                                                                                            <  2:886758/150‑1 (MQ=255)
     tCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCg                                                                                                                                            <  1:3587289/150‑1 (MQ=255)
          tGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCttgt                                                                                                                                       <  1:1325405/150‑1 (MQ=255)
           gAAGCAGCTTATTGAAATGAACCGAATCACCATTTCCAATGGCGATATTGGTTTTAACTTCTCTGACGGCAATCTGATTAAGTAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCttgtt                                                                                                                                      <  1:2540609/150‑1 (MQ=255)
              gcagcTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGt                                                                                                                                   >  2:112549/1‑150 (MQ=255)
              gcagcTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGt                                                                                                                                   <  1:4099997/150‑1 (MQ=255)
               cagcTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTc                                                                                                                                  <  2:5035953/150‑1 (MQ=255)
               cagcTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTc                                                                                                                                  <  1:2318302/150‑1 (MQ=255)
                agcTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCg                                                                                                                                 <  1:2258934/150‑1 (MQ=255)
                 gcTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGa                                                                                                                                >  2:1175464/1‑150 (MQ=255)
                 gcTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGa                                                                                                                                <  1:1438327/150‑1 (MQ=255)
                  cTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGAt                                                                                                                               >  1:3037811/1‑150 (MQ=255)
                  cTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGAt                                                                                                                               >  1:1361852/1‑150 (MQ=255)
                    tATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATaa                                                                                                                             <  1:2269632/150‑1 (MQ=255)
                     aTTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAAc                                                                                                                            >  2:1259072/1‑150 (MQ=255)
                       tGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACaa                                                                                                                          >  1:1573851/1‑150 (MQ=255)
                         aaaTGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATa                                                                                                                        <  2:476623/150‑1 (MQ=255)
                         aaaTGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATa                                                                                                                        >  2:1435444/1‑150 (MQ=255)
                         aaaTGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATa                                                                                                                        <  1:4667714/150‑1 (MQ=255)
                          aaTGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAg                                                                                                                       >  1:647882/1‑150 (MQ=255)
                          aaTGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAg                                                                                                                       <  2:3537142/150‑1 (MQ=255)
                            tGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCg                                                                                                                     >  2:4048957/1‑150 (MQ=255)
                             gAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGa                                                                                                                    >  2:1590963/1‑150 (MQ=255)
                             gAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGa                                                                                                                    >  2:670103/1‑150 (MQ=255)
                               aCCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAg                                                                                                                  <  2:3824834/150‑1 (MQ=255)
                                 cGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGt                                                                                                                >  1:2576737/1‑150 (MQ=255)
                                 cGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGt                                                                                                                <  2:2390235/150‑1 (MQ=255)
                                 cGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGt                                                                                                                >  1:96040/1‑150 (MQ=255)
                                  gAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTg                                                                                                               <  1:4665772/150‑1 (MQ=255)
                                    aTCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGaa                                                                                                             >  1:5097790/1‑150 (MQ=255)
                                      cACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATa                                                                                                           <  1:479664/150‑1 (MQ=255)
                                       aCCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATAc                                                                                                          >  2:208473/1‑150 (MQ=255)
                                        ccATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACg                                                                                                         <  1:1652702/150‑1 (MQ=255)
                                         cATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGc                                                                                                        >  2:4168116/1‑150 (MQ=255)
                                         cATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGc                                                                                                        >  1:5120154/1‑150 (MQ=255)
                                          aTTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGcc                                                                                                       <  2:1788342/150‑1 (MQ=255)
                                            tGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCat                                                                                                     >  1:4301655/1‑150 (MQ=255)
                                             gCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCatt                                                                                                    <  2:3039378/150‑1 (MQ=255)
                                             gCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCatt                                                                                                    <  2:4598465/150‑1 (MQ=255)
                                              ccAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCatta                                                                                                   >  2:594971/1‑150 (MQ=255)
                                              ccAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCatta                                                                                                   >  1:823759/1‑150 (MQ=255)
                                               cAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCattat                                                                                                  <  1:1069550/150‑1 (MQ=255)
                                                aaTGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCattatt                                                                                                 >  2:3108098/1‑150 (MQ=255)
                                                aaTGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCattatt                                                                                                 >  2:1920031/1‑150 (MQ=255)
                                                aaTGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCattatt                                                                                                 <  1:5042953/150‑1 (MQ=255)
                                                 aTGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTc                                                                                                >  2:3224363/1‑150 (MQ=255)
                                                 aTGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTc                                                                                                >  2:2982111/1‑150 (MQ=255)
                                                   ggCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTccc                                                                                              <  1:2814210/150‑1 (MQ=255)
                                                   ggCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTccc                                                                                              <  1:560478/150‑1 (MQ=255)
                                                      gATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGcc                                                                                           <  2:4351165/150‑1 (MQ=255)
                                                       atatTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCa                                                                                          <  2:1001097/150‑1 (MQ=255)
                                                        tatTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAg                                                                                         <  2:3234915/150‑1 (MQ=255)
                                                         atTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGc                                                                                        >  1:321994/1‑150 (MQ=255)
                                                          ttGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCg                                                                                       <  2:4946009/150‑1 (MQ=255)
                                                           tGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGt                                                                                      <  1:2919684/150‑1 (MQ=255)
                                                               tttAACTTCACTGACGGGAATCTGATTATGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGcc                                                                                  <  2:3910670/150‑1 (MQ=255)
                                                                 tAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGa                                                                                <  2:2897502/150‑1 (MQ=255)
                                                                 tAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGa                                                                                <  2:738790/150‑1 (MQ=255)
                                                                  aaCTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGAt                                                                               >  1:1120312/1‑150 (MQ=255)
                                                                    cTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATaa                                                                             >  2:1878089/1‑150 (MQ=255)
                                                                      tCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATaaaa                                                                           >  1:4889238/1‑150 (MQ=255)
                                                                          tGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGc                                                                       >  2:5082371/1‑150 (MQ=255)
                                                                           gACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCt                                                                      >  2:2372050/1‑150 (MQ=255)
                                                                                 aaTCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAgcgc                                                                >  1:204934/1‑150 (MQ=255)
                                                                                  aTCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGcgcg                                                               <  1:2454683/150‑1 (MQ=255)
                                                                                     tGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATg                                                            <  2:4101256/150‑1 (MQ=255)
                                                                                            gaagaTTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCcagcag                                                     >  1:1001432/1‑150 (MQ=255)
                                                                                             aagaTTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGt                                                    >  2:1053587/1‑150 (MQ=255)
                                                                                               gaTTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGtat                                                  >  1:336730/1‑150 (MQ=255)
                                                                                               gaTTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGtat                                                  >  2:4285412/1‑150 (MQ=255)
                                                                                                  ttttGTCGATAAGCTCACACAGGCGCAGTTGATTACTGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAgt                                               <  1:4409067/150‑1 (MQ=255)
                                                                                                  ttttGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAgt                                               <  1:4652376/150‑1 (MQ=255)
                                                                                                      gTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTg                                           <  1:3688420/150‑1 (MQ=255)
                                                                                                      gTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTg                                           >  1:1597995/1‑150 (MQ=255)
                                                                                                        cGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGca                                         <  2:2837742/150‑1 (MQ=255)
                                                                                                         gATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGcac                                        >  1:682968/1‑150 (MQ=255)
                                                                                                           tAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAg                                      <  2:1363432/150‑1 (MQ=255)
                                                                                                            aaGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAgc                                     <  2:449482/150‑1 (MQ=255)
                                                                                                               cTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAgcgcg                                  >  2:3479735/1‑150 (MQ=255)
                                                                                                                tCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAgcgcgc                                 >  2:4455515/1‑150 (MQ=255)
                                                                                                                   cacaGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGAGCGCTCa                              >  2:2619770/1‑150 (MQ=255)
                                                                                                                    acaGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAg                             >  1:2433523/1‑150 (MQ=255)
                                                                                                                     caGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAgc                            >  2:2138335/1‑150 (MQ=255)
                                                                                                                      aGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAgcg                           <  2:2445357/150‑1 (MQ=255)
                                                                                                                        gcgcAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGcc                         <  1:3259237/150‑1 (MQ=255)
                                                                                                                        gcgcAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGcc                         <  1:1801872/150‑1 (MQ=255)
                                                                                                                           cAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGaa                      <  2:2816524/150‑1 (MQ=255)
                                                                                                                            aGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAag                     >  1:459115/1‑150 (MQ=255)
                                                                                                                                 aTTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAg                >  1:812578/1‑150 (MQ=255)
                                                                                                                                  ttAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAgg               <  2:2328994/150‑1 (MQ=255)
                                                                                                                                    aaTGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGAt             >  2:2550792/1‑150 (MQ=255)
                                                                                                                                        gtcgtcTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAg         <  1:4160406/150‑1 (MQ=255)
                                                                                                                                         tcgtcTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAga        >  2:3021855/1‑150 (MQ=255)
                                                                                                                                           gtcTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAgacg      <  2:3495406/150‑1 (MQ=255)
                                                                                                                                            tcTGGCTAGTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAgacga     <  1:3108098/150‑1 (MQ=255)
                                                                                                                                             cTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAgacgac    >  1:3535964/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAAGCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCC  >  CP009273/372435‑372727

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: