Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,256,589 G→T 6.5% E294* (GAA→TAA)  alr → alanine racemase, biosynthetic, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,256,5890GT 6.5% 1475.7 / 45.0 378E294* (GAA→TAA) alralanine racemase, biosynthetic, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base G (172/181);  new base T (11/14);  total (183/195)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 391 total aligned reads displayed.

CGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAGCGCTTAC  >  CP009273/4256442‑4256738
                                                                                                                                                   |                                                                                                                                                     
   gAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTTGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTa                                                                                                                                                  <  2:2613107/150‑1 (MQ=255)
     gCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTAcc                                                                                                                                                <  2:4862610/150‑1 (MQ=255)
       aTAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGa                                                                                                                                              >  2:1146901/1‑150 (MQ=255)
         aaaGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGAtt                                                                                                                                            >  2:3000708/1‑150 (MQ=255)
            gCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTc                                                                                                                                         >  2:3306690/1‑150 (MQ=255)
              cGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCgg                                                                                                                                       >  2:1198391/1‑150 (MQ=255)
                     ccTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCTAAGTACCGATTGTCGGGCGCGTg                                                                                                                                <  2:3982215/150‑1 (MQ=255)
                      cTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTgg                                                                                                                               <  2:3637020/150‑1 (MQ=255)
                      cTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTgg                                                                                                                               <  2:3590845/150‑1 (MQ=255)
                        gTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCg                                                                                                                             <  1:781448/150‑1 (MQ=255)
                         ttGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGa                                                                                                                            >  1:695565/1‑150 (MQ=255)
                            gTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATgg                                                                                                                         >  2:2513136/1‑150 (MQ=255)
                               aTGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGata                                                                                                                      >  2:5125945/1‑150 (MQ=255)
                               aTGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGata                                                                                                                      <  1:2317148/150‑1 (MQ=255)
                                 ggtggAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCTAAGTACCGATTGTCGGGCGCGTGGCGATGGATATg                                                                                                                    >  1:1197856/1‑150 (MQ=255)
                                   tggAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGAt                                                                                                                  >  2:2882968/1‑150 (MQ=255)
                                    ggAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATc                                                                                                                 <  2:1604234/150‑1 (MQ=255)
                                       aCCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGc                                                                                                              >  1:1385903/1‑150 (MQ=255)
                                          tGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTa                                                                                                           <  1:2539835/150‑1 (MQ=255)
                                          tGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTa                                                                                                           >  1:1741327/1‑150 (MQ=255)
                                           gggTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAg                                                                                                          >  2:3076260/1‑150 (MQ=255)
                                           gggTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAg                                                                                                          <  1:119068/150‑1 (MQ=255)
                                           gggTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAg                                                                                                          >  1:530064/1‑150 (MQ=255)
                                              tAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACt                                                                                                       >  1:2493451/1‑150 (MQ=255)
                                                aGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTa                                                                                                     <  2:3737538/150‑1 (MQ=255)
                                                 gCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAg                                                                                                    >  2:1540853/1‑150 (MQ=255)
                                                  cGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAgg                                                                                                   <  2:4710326/150‑1 (MQ=255)
                                                   gAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGt                                                                                                  <  1:1662461/150‑1 (MQ=255)
                                                     aCGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTcc                                                                                                <  2:1634577/150‑1 (MQ=255)
                                                        tGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCcaca                                                                                             >  2:2350703/1‑150 (MQ=255)
                                                        tGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCcaca                                                                                             <  1:2781509/150‑1 (MQ=255)
                                                         gATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAg                                                                                            <  2:2076713/150‑1 (MQ=255)
                                                          aTACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAgg                                                                                           <  2:1567642/150‑1 (MQ=255)
                                                            aCCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGgcg                                                                                         >  1:3398004/1‑150 (MQ=255)
                                                              ccGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCTAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCa                                                                                       <  1:3606634/150‑1 (MQ=255)
                                                                gTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCTAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAgg                                                                                     <  1:2885851/150‑1 (MQ=255)
                                                                gTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAgg                                                                                     >  2:1978777/1‑150 (MQ=255)
                                                                  cTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCTAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGAc                                                                                   <  2:4758508/150‑1 (MQ=255)
                                                                  cTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGAc                                                                                   <  1:4601617/150‑1 (MQ=255)
                                                                   ttGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCTAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACa                                                                                  >  1:2906990/1‑150 (MQ=255)
                                                                   ttGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACa                                                                                  >  2:1800666/1‑150 (MQ=255)
                                                                   ttGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACa                                                                                  <  1:4896866/150‑1 (MQ=255)
                                                                    tGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACaa                                                                                 <  1:3237964/150‑1 (MQ=255)
                                                                    tGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACaa                                                                                 >  1:861671/1‑150 (MQ=255)
                                                                     ggCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACaag                                                                                <  2:2144174/150‑2 (MQ=255)
                                                                       cGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGc                                                                              <  2:2682934/150‑1 (MQ=255)
                                                                        gTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGcc                                                                             >  2:3879536/1‑150 (MQ=255)
                                                                          aGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCgg                                                                           >  1:4813912/1‑150 (MQ=255)
                                                                          aGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCgg                                                                           >  2:473327/1‑150 (MQ=255)
                                                                             cgcgATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACCAAGCCGGGGa                                                                        >  1:3024264/1‑150 (MQ=255)
                                                                              gcgATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGACCGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGAt                                                                       >  1:4785474/1‑150 (MQ=255)
                                                                               cgATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATc                                                                      >  1:4296902/1‑150 (MQ=255)
                                                                                gATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCTAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATcc                                                                     <  1:2228637/150‑1 (MQ=255)
                                                                                gATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATcc                                                                     <  2:4212823/150‑1 (MQ=255)
                                                                                 aTGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCg                                                                    <  1:4733529/150‑1 (MQ=255)
                                                                                 aTGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCg                                                                    >  2:1667759/1‑150 (MQ=255)
                                                                                 aTGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCg                                                                    <  2:1413874/150‑1 (MQ=255)
                                                                                  tGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCgg                                                                   >  1:3350452/1‑150 (MQ=255)
                                                                                   gggCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGt                                                                  <  1:881854/150‑1 (MQ=255)
                                                                                    ggCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTc                                                                 >  1:4239356/1‑150 (MQ=255)
                                                                                     gCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCa                                                                <  1:611023/150‑1 (MQ=255)
                                                                                         tGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCAtttt                                                            <  2:4477253/150‑1 (MQ=255)
                                                                                              aTGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATggg                                                       <  1:5081572/150‑1 (MQ=255)
                                                                                               tGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCTAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATgggg                                                      <  1:17851/150‑1 (MQ=255)
                                                                                               tGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATgggg                                                      <  1:3559617/150‑1 (MQ=255)
                                                                                                ggTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGc                                                     >  2:701386/1‑150 (MQ=255)
                                                                                                 gTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCg                                                    >  1:1624871/1‑150 (MQ=255)
                                                                                                 gTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCg                                                    <  1:3829131/150‑1 (MQ=255)
                                                                                                    aTCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAg                                                 <  1:2925907/150‑1 (MQ=255)
                                                                                                      ccGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGt                                               <  1:1017148/150‑1 (MQ=255)
                                                                                                       cgcgcgCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGtt                                              <  2:492464/150‑1 (MQ=255)
                                                                                                        gcgcgcCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGttt                                             <  2:1307650/150‑1 (MQ=255)
                                                                                                        gcgcgcCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGttt                                             <  1:4446431/150‑1 (MQ=255)
                                                                                                        gcgcgcCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGttt                                             <  2:4586215/150‑1 (MQ=255)
                                                                                                           cgccgcgccgTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGcc                                          <  2:4998583/150‑1 (MQ=255)
                                                                                                            gccgcgccgTCCGGTCCGCCAGTGCTGGTGAACGGTCGCTAAGGACCGATTGTCGGGCGCGTGGCGATGGCTATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGccc                                         <  1:3735155/150‑1 (MQ=255)
                                                                                                               gcgccgTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTa                                      <  1:398550/150‑1 (MQ=255)
                                                                                                               gcgccgTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTa                                      <  1:3192782/150‑1 (MQ=255)
                                                                                                                 gccgTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGa                                    <  2:4145722/150‑1 (MQ=255)
                                                                                                                  ccgTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGaa                                   <  2:4549402/150‑1 (MQ=255)
                                                                                                                  ccgTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGaa                                   <  2:2587278/150‑1 (MQ=255)
                                                                                                                        ggTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATc                             <  1:3088970/150‑1 (MQ=255)
                                                                                                                         gTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCg                            >  1:4165877/1‑150 (MQ=255)
                                                                                                                         gTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCg                            >  1:2593271/1‑150 (MQ=255)
                                                                                                                         gTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCg                            <  1:2483271/150‑1 (MQ=255)
                                                                                                                            cGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTg                         >  2:4410196/1‑150 (MQ=255)
                                                                                                                            cGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTg                         >  2:444273/1‑150 (MQ=255)
                                                                                                                              ccAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGaa                       >  2:554666/1‑150 (MQ=255)
                                                                                                                              ccAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGaa                       >  1:1912076/1‑150 (MQ=255)
                                                                                                                                aGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAAt                     >  1:4547609/1‑150 (MQ=255)
                                                                                                                                    cTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACg                 <  2:192953/150‑1 (MQ=255)
                                                                                                                                    cTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACg                 >  1:903791/1‑150 (MQ=255)
                                                                                                                                     tGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGa                >  2:3500261/1‑150 (MQ=255)
                                                                                                                                           aCGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTaa          >  1:2075368/1‑150 (MQ=255)
                                                                                                                                              gTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAgcg       <  2:2374146/150‑1 (MQ=255)
                                                                                                                                              gTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAgcg       <  2:2194901/150‑1 (MQ=255)
                                                                                                                                              gTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAgcg       <  2:3812019/150‑1 (MQ=255)
                                                                                                                                                  cgAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAGCGCTTa   <  2:317144/150‑1 (MQ=255)
                                                                                                                                                  cgAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAGCGCTTa   >  2:4247247/1‑150 (MQ=255)
                                                                                                                                                   gAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGTCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAGCTCTTAc  >  1:4287470/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAGTGCTGGTGAACGGTCGCGAAGTACCGATTGTCGGGCGCGTGGCGATGGATATGATCTGCGTAGACTTAGGTCCACAGGCGCAGGACAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAGCGCTTAC  >  CP009273/4256442‑4256738

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: