Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,122,781 C→A 16.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,7810CA16.5% 881.5 / 105.1 276R80S (CGC→AGC) metF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (111/119);  new base A (34/12);  total (145/131)
Fisher's exact test for biased strand distribution p-value = 1.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 284 total aligned reads displayed.

TCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGAAGTGGTAAGCCAGAA  >  CP009273/4122632‑4122930
                                                                                                                                                     |                                                                                                                                                     
 cccGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCg                                                                                                                                                      >  2:2385433/1‑150 (MQ=255)
   cgccgcGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCa                                                                                                                                                    <  1:3834666/150‑1 (MQ=255)
    gccgcGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCAc                                                                                                                                                   <  1:538735/150‑1 (MQ=255)
       gcgTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTgg                                                                                                                                                <  1:1657525/150‑1 (MQ=255)
        cgTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGt                                                                                                                                               <  2:2877913/150‑1 (MQ=255)
        cgTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGt                                                                                                                                               >  2:710942/1‑150 (MQ=255)
         gTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTc                                                                                                                                              >  1:1291612/1‑150 (MQ=255)
          tACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCt                                                                                                                                             >  1:2015374/1‑150 (MQ=255)
           aCCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTg                                                                                                                                            <  1:3286194/150‑1 (MQ=255)
            cctgTGAAATGGAGCAGACACTGTGGATCTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTgg                                                                                                                                           <  2:1709906/147‑1 (MQ=255)
             cAGTGAAATGGAGCAGACCCTGTGGACCTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGa                                                                                                                                          >  2:323391/1‑150 (MQ=255)
             cAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGa                                                                                                                                          <  1:3961432/150‑1 (MQ=255)
             cAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGa                                                                                                                                          <  1:745358/150‑1 (MQ=255)
                tGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACTGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGc                                                                                                                                       <  1:466774/150‑1 (MQ=255)
                 gAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCg                                                                                                                                      <  2:449119/150‑1 (MQ=255)
                      ggAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCAcc                                                                                                                                 <  1:230847/150‑1 (MQ=255)
                      ggAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCAcc                                                                                                                                 <  2:3513018/150‑1 (MQ=255)
                       gAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCg                                                                                                                                <  2:4155113/150‑1 (MQ=255)
                          cAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCAt                                                                                                                             >  2:253903/1‑150 (MQ=255)
                           aGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATc                                                                                                                            >  2:1387874/1‑150 (MQ=255)
                            gACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCt                                                                                                                           <  2:3166637/150‑1 (MQ=255)
                            gACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCt                                                                                                                           >  1:2150022/1‑150 (MQ=255)
                               ccTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTAc                                                                                                                        >  2:637489/1‑150 (MQ=255)
                                cTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGCATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACt                                                                                                                       <  1:88424/150‑1 (MQ=255)
                                    ggAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCa                                                                                                                   <  2:4618106/150‑1 (MQ=255)
                                     gAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCAt                                                                                                                  >  1:3152593/1‑150 (MQ=255)
                                      aaCTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCAtt                                                                                                                 <  1:1634619/150‑1 (MQ=255)
                                      aaCTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCAtt                                                                                                                 >  1:1540596/1‑150 (MQ=255)
                                         tCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGAt                                                                                                              >  2:3019741/1‑150 (MQ=255)
                                         tCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGAt                                                                                                              <  1:4093451/150‑1 (MQ=255)
                                          ccATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATg                                                                                                             <  2:179492/150‑1 (MQ=255)
                                          ccATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATg                                                                                                             <  2:1356623/150‑1 (MQ=255)
                                           cATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGc                                                                                                            >  2:15329/1‑150 (MQ=255)
                                             tcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGa                                                                                                          >  2:310275/1‑150 (MQ=255)
                                             tcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGa                                                                                                          <  2:605975/150‑1 (MQ=255)
                                                 tcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTCTGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGcc                                                                                                      <  2:4576524/150‑1 (MQ=255)
                                                 tcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGcc                                                                                                      >  2:980285/1‑150 (MQ=255)
                                                  cgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGccc                                                                                                     >  2:386349/1‑150 (MQ=255)
                                                  cgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGccc                                                                                                     <  2:1662353/150‑1 (MQ=255)
                                                   gCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCcg                                                                                                    <  1:3973486/150‑1 (MQ=255)
                                                   gCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCcg                                                                                                    >  1:1892620/1‑150 (MQ=255)
                                                    ccTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCcga                                                                                                   <  2:4375658/150‑1 (MQ=255)
                                                        agcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAg                                                                                               >  1:4422891/1‑150 (MQ=255)
                                                           agcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTg                                                                                            <  1:1078852/150‑1 (MQ=255)
                                                            gcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTgc                                                                                           <  1:3421610/150‑1 (MQ=255)
                                                             ccTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTgcg                                                                                          <  2:1845519/150‑1 (MQ=255)
                                                              cTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTgcgc                                                                                         <  1:2745447/150‑1 (MQ=255)
                                                              cTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTgcgc                                                                                         >  1:4525925/1‑150 (MQ=255)
                                                                  aaCCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCa                                                                                     <  2:2764330/150‑1 (MQ=255)
                                                                   aCCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCAt                                                                                    >  1:3802754/1‑150 (MQ=255)
                                                                    ccGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCAtt                                                                                   >  1:2602254/1‑150 (MQ=255)
                                                                      gAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGc                                                                                 >  1:4222271/1‑150 (MQ=255)
                                                                       aaGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCa                                                                                >  1:1143509/1‑150 (MQ=255)
                                                                        aGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCAc                                                                               >  2:556636/1‑150 (MQ=255)
                                                                          tttGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCAcgc                                                                             <  1:2331810/150‑1 (MQ=255)
                                                                              tATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGact                                                                         <  2:2707940/150‑1 (MQ=255)
                                                                                 cGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGactact                                                                      <  1:2478444/150‑1 (MQ=255)
                                                                                   gTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTgg                                                                    <  2:1945786/150‑1 (MQ=255)
                                                                                   gTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTgg                                                                    <  2:2481802/150‑1 (MQ=255)
                                                                                       ccTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATa                                                                >  1:3549544/1‑150 (MQ=255)
                                                                                        cTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATaa                                                               <  2:1905945/150‑1 (MQ=255)
                                                                                        cTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATaa                                                               >  2:4286956/1‑150 (MQ=255)
                                                                                         tATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGACTAAc                                                              >  2:401088/1‑150 (MQ=255)
                                                                                         tATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAAc                                                              >  1:3645954/1‑150 (MQ=255)
                                                                                         tATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAAc                                                              >  1:428552/1‑150 (MQ=255)
                                                                                          aTGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACg                                                             <  1:181546/150‑1 (MQ=255)
                                                                                          aTGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACg                                                             <  1:2385433/150‑1 (MQ=255)
                                                                                           tGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACgg                                                            >  1:4476315/1‑150 (MQ=255)
                                                                                             gcgcGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTa                                                          <  2:1205000/150‑1 (MQ=255)
                                                                                                cgAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTc                                                       <  1:224329/150‑1 (MQ=255)
                                                                                                cgAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTc                                                       >  1:1465890/1‑150 (MQ=255)
                                                                                                 gAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGACTAACGGTCTTCg                                                      >  2:2054715/1‑150 (MQ=255)
                                                                                                 gAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCg                                                      <  2:3769836/150‑1 (MQ=255)
                                                                                                    cTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCa                                                   <  1:283560/150‑1 (MQ=255)
                                                                                                       cGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCatat                                                <  2:731656/150‑1 (MQ=255)
                                                                                                           gAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTg                                            <  1:4176416/150‑1 (MQ=255)
                                                                                                           gAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTg                                            <  2:805563/150‑1 (MQ=255)
                                                                                                              cgcgACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGgcg                                         >  1:1242956/1‑150 (MQ=255)
                                                                                                              cgcgACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGgcg                                         <  1:2996587/150‑1 (MQ=255)
                                                                                                               gcgACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGgcgc                                        >  1:251395/1‑150 (MQ=255)
                                                                                                                 gACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTg                                      >  1:2629906/1‑150 (MQ=255)
                                                                                                                 gACCGTACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTg                                      >  2:1622294/1‑150 (MQ=255)
                                                                                                                   ccGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCg                                    >  2:2572176/1‑150 (MQ=255)
                                                                                                                     gTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTg                                  <  2:3615556/150‑1 (MQ=255)
                                                                                                                     gTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTg                                  <  2:1075639/150‑1 (MQ=255)
                                                                                                                      tACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTgg                                 >  1:3697121/1‑150 (MQ=255)
                                                                                                                      tACGCACAGCATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTgg                                 >  1:4368624/1‑150 (MQ=255)
                                                                                                                       aCGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGc                                >  2:921226/1‑150 (MQ=255)
                                                                                                                             aGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTg                          >  1:444681/1‑150 (MQ=255)
                                                                                                                               cATTATTAAAGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTgcc                        >  2:1669254/1‑150 (MQ=255)
                                                                                                                                   attaAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGccgccg                    >  2:2954292/1‑150 (MQ=255)
                                                                                                                                    ttaAAGGCATTAAAGATAGCACTGGTCAGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCgg                   >  1:3884585/1‑150 (MQ=255)
                                                                                                                                      aaaGGCATTAAAGATAGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGa                 >  2:4477856/1‑150 (MQ=255)
                                                                                                                                       aaGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGaa                <  2:2231690/150‑1 (MQ=255)
                                                                                                                                       aaGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGaa                <  2:3324433/150‑1 (MQ=255)
                                                                                                                                        aGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGAAg               <  1:2194295/150‑1 (MQ=255)
                                                                                                                                        aGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGAAg               >  1:1710191/1‑150 (MQ=255)
                                                                                                                                        aGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGAAg               >  2:877463/1‑150 (MQ=255)
                                                                                                                                         ggCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGAAGt              <  2:1323154/150‑1 (MQ=255)
                                                                                                                                                aaGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCTCCATTGCACGCGACTACTGGAATAACGGTATTCGTCATATCGTGGCGCTGCGTGGCGATCTGCCGCCGGGAAGTGGTTAGc       <  2:1557522/150‑1 (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: