Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 943,573 G→T 100% A297S (GCC→TCC)  ycaM → putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273943,5730GT95.1% 849.6 / 27.6 266A297S (GCC→TCC) ycaMputative transporter
Reads supporting (aligned to +/- strand):  ref base G (8/5);  new base T (140/113);  total (148/118)
Fisher's exact test for biased strand distribution p-value = 7.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 268 total aligned reads displayed.

TTTGTGCCATTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACGCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGTTAATGGCTATTTT  >  CP009273/943427‑943722
                                                                                                                                                  |                                                                                                                                                     
  tgtgCCATTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCAtt                                                                                                                                                  <  1:995666/150‑1 (MQ=255)
    tgCCATTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGc                                                                                                                                                <  1:4634072/150‑1 (MQ=255)
    tgCCATTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGc                                                                                                                                                >  2:3492177/1‑150 (MQ=255)
     gCCATTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCg                                                                                                                                               >  1:918370/1‑150 (MQ=255)
     gCCATTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCg                                                                                                                                               >  2:141778/1‑150 (MQ=255)
     gCCATTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCg                                                                                                                                               <  2:1859108/150‑1 (MQ=255)
         ttCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATa                                                                                                                                           <  1:1294622/150‑1 (MQ=255)
           cTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATAcc                                                                                                                                         >  2:1936307/1‑150 (MQ=255)
             gggCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGTTGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCt                                                                                                                                       <  1:2092695/150‑1 (MQ=255)
               gCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTgg                                                                                                                                     <  1:2162370/150‑1 (MQ=255)
                cTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTggg                                                                                                                                    <  1:4348636/150‑1 (MQ=255)
                 tCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGa                                                                                                                                   >  2:840828/1‑150 (MQ=255)
                 tCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGa                                                                                                                                   >  1:1394031/1‑150 (MQ=255)
                    cTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACaa                                                                                                                                <  2:716574/150‑1 (MQ=255)
                     tGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAg                                                                                                                               >  1:120653/1‑150 (MQ=255)
                     tGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAg                                                                                                                               >  2:4410803/1‑150 (MQ=255)
                        cGATGGGGATGATGTTTGATTCTCGTAATATCACGGATGACTAAATGACCAACGGTCAGTATAACGCCTATCAGAAGCTGGGCGAGAATTACAACATGCGTAATACTTTAATGGTGATTAACTCCAATCCTAATACACTGGGACAAGTAg                                                                                                                            >  1:2671318/1‑150 (MQ=255)
                         gATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAgc                                                                                                                           >  1:2582380/1‑150 (MQ=255)
                          aTGGGGTTGATGTTTGTTTCTCGTAATATCCCGGATGACTTAATGACCATCGGTCAGTATTACTCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTTATGGTGATTTTCTCCATTGCGAATACCCTGGGACAAGTAgcg                                                                                                                          <  1:1969415/150‑1 (MQ=255)
                             gggATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAgcggcg                                                                                                                       <  1:4389155/150‑1 (MQ=255)
                                    tgTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTAt                                                                                                                >  1:2831044/1‑150 (MQ=255)
                                     gTTTGATTCGCGTAATCTCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGTTGCGTGGCGTGTATTACAACATGGGTTATACTTTAATGGTGATTTACTCCATTGCGATTTCCCTGGGACAAGTAGCGGCGCTGGTAtt                                                                                                               <  1:3361790/150‑1 (MQ=255)
                                           ttCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGAt                                                                                                         <  2:2292927/150‑1 (MQ=255)
                                           ttCGCGTAATATCCCGGATGACTGAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGAt                                                                                                         <  2:2760799/150‑1 (MQ=255)
                                            tCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGAtt                                                                                                        >  2:303340/1‑150 (MQ=255)
                                            tCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGAtt                                                                                                        >  2:869605/1‑150 (MQ=255)
                                             cgcgTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACCTGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTg                                                                                                       <  2:2444383/150‑1 (MQ=255)
                                                 tAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGc                                                                                                   >  2:3029166/1‑150 (MQ=255)
                                                  aaTATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGcc                                                                                                  >  2:1356873/1‑150 (MQ=255)
                                                  aaTATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGcc                                                                                                  >  1:2967841/1‑150 (MQ=255)
                                                    tatCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTTATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGcccc                                                                                                <  1:4193659/150‑1 (MQ=255)
                                                        ccGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTCTTCTCGATTGATGCCCCGCtt                                                                                            >  2:3981439/1‑150 (MQ=255)
                                                        ccGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCtt                                                                                            >  2:3661892/1‑150 (MQ=255)
                                                          ggATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTaa                                                                                          <  2:349931/150‑1 (MQ=255)
                                                           gATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTaaa                                                                                         <  1:430717/150‑1 (MQ=255)
                                                            aTGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAg                                                                                        <  2:1069361/150‑1 (MQ=255)
                                                            aTGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAg                                                                                        <  2:1038912/150‑1 (MQ=255)
                                                            aTGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGAGGCCCCGCTTAAAg                                                                                        <  2:3530035/150‑1 (MQ=255)
                                                              gACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTg                                                                                      >  2:2394157/1‑150 (MQ=255)
                                                               aCTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGc                                                                                     >  2:3718004/1‑150 (MQ=255)
                                                               aCTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGc                                                                                     >  1:3810145/1‑150 (MQ=255)
                                                                cTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCt                                                                                    >  1:3500198/1‑150 (MQ=255)
                                                                   aaTGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTAtt                                                                                 >  2:3795918/1‑150 (MQ=255)
                                                                     tGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAg                                                                               <  2:1107654/150‑1 (MQ=255)
                                                                      gACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAgg                                                                              <  1:4530129/150‑1 (MQ=255)
                                                                       aCCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGt                                                                             >  1:3949347/1‑150 (MQ=255)
                                                                                cAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGAc                                                                    <  1:522757/150‑1 (MQ=255)
                                                                                 aGTATTACGCCTTTCAGAAGCTGGGCGAGTCTTACAACATGGGTAATACTTTAATGGTGATTTACGCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACa                                                                   >  2:2719516/1‑150 (MQ=255)
                                                                                 aGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACa                                                                   >  2:4463729/1‑150 (MQ=255)
                                                                                   tATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGc                                                                 <  1:1925382/150‑1 (MQ=255)
                                                                                   tATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGc                                                                 >  1:3464206/1‑150 (MQ=255)
                                                                                      tACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCaaa                                                              <  2:3317399/150‑1 (MQ=255)
                                                                                       aCGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAat                                                             <  1:152295/150‑1 (MQ=255)
                                                                                       aCGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAat                                                             <  2:2194176/150‑1 (MQ=255)
                                                                                        cGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAata                                                            <  1:2338823/150‑1 (MQ=255)
                                                                                         gCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAatat                                                           >  1:145938/1‑150 (MQ=255)
                                                                                           cTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAatatat                                                         >  2:1700700/1‑150 (MQ=255)
                                                                                           cTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAatatat                                                         <  1:2310684/150‑1 (MQ=255)
                                                                                                     cTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGt                                               >  2:3012925/1‑150 (MQ=255)
                                                                                                       gggCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGttt                                             <  2:1946757/150‑1 (MQ=255)
                                                                                                          cGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATg                                          >  1:4337481/1‑150 (MQ=255)
                                                                                                            aGTATTACAACATGGGTAATACTTTAATGGTGATTTACGCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTc                                        >  1:3670926/1‑150 (MQ=255)
                                                                                                            aGTATTACAACATGGGTAATACTTTAATGGTGATTTACGCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTc                                        >  2:2111440/1‑150 (MQ=255)
                                                                                                               aTTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTa                                     >  1:966312/1‑150 (MQ=255)
                                                                                                               aTTACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTa                                     >  1:3863708/1‑150 (MQ=255)
                                                                                                                 tACAACATGGGTAATACTTTTTTGGTGATTTACTCGATTGAGAATTCCCTTTTACAAGTAGCGGCGCTGGTATTCTAGATTGATTCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAATTATATTCAAGCCAGTTTATTTCGTaac                                   <  2:2174891/150‑3 (MQ=255)
                                                                                                                 tACAACATGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTAcc                                   >  2:1723659/1‑150 (MQ=255)
                                                                                                                  acaacaTGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCa                                  >  1:1543646/1‑150 (MQ=255)
                                                                                                                  acaacaTGGGTAATACTATAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCa                                  <  2:4159624/150‑1 (MQ=255)
                                                                                                                   caacaTGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCaa                                 >  1:1125927/1‑150 (MQ=255)
                                                                                                                   caacaTGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCaa                                 >  1:3222256/1‑150 (MQ=255)
                                                                                                                    aacaTGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAAc                                <  2:3564747/150‑1 (MQ=255)
                                                                                                                     acaTGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACg                               >  1:3826782/1‑150 (MQ=255)
                                                                                                                        tGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCtt                            <  1:204843/150‑1 (MQ=255)
                                                                                                                        tGGGTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCtt                            >  1:3004570/1‑150 (MQ=255)
                                                                                                                         gggTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTc                           <  2:4337131/150‑1 (MQ=255)
                                                                                                                         gggTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTc                           >  1:2213087/1‑150 (MQ=255)
                                                                                                                          ggTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCt                          <  1:2245699/150‑1 (MQ=255)
                                                                                                                          ggTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCt                          <  2:2844080/150‑1 (MQ=255)
                                                                                                                          ggTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCt                          >  1:4427850/1‑150 (MQ=255)
                                                                                                                           gTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTg                         >  1:2445397/1‑150 (MQ=255)
                                                                                                                           gTAATACTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTg                         <  2:437322/150‑1 (MQ=255)
                                                                                                                                aCTTTAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACg                    >  1:1393823/1‑150 (MQ=255)
                                                                                                                                  tttAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGcc                  >  2:62326/1‑150 (MQ=255)
                                                                                                                                  tttAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGcc                  >  2:61371/1‑150 (MQ=255)
                                                                                                                                   ttAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGccc                 <  1:3164748/150‑1 (MQ=255)
                                                                                                                                   ttAATGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGccc                 >  2:3933165/1‑150 (MQ=255)
                                                                                                                                   ttAATGGTGATGTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGccc                 <  2:1199105/150‑1 (MQ=255)
                                                                                                                                    tAATGGTGTTTTACTCCATTGCGAATACCCTGGGACAAGTCGCGGCGCTGGTATTCTCGATTGATTCCACGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCg                <  2:2080112/150‑1 (MQ=255)
                                                                                                                                      aTGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGtt              >  1:4090523/1‑150 (MQ=255)
                                                                                                                                      aTGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGtt              >  1:2692844/1‑150 (MQ=255)
                                                                                                                                      aTGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGtt              >  1:3137663/1‑150 (MQ=255)
                                                                                                                                      aTGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGtt              >  1:794127/1‑150 (MQ=255)
                                                                                                                                       tGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGTTa             <  1:1848129/150‑1 (MQ=255)
                                                                                                                                       tGGTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGTTa             >  2:1614847/1‑150 (MQ=255)
                                                                                                                                        ggTGATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGTTaa            <  2:4593739/150‑1 (MQ=255)
                                                                                                                                           gATTTACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGTTAATgg         >  2:3522116/1‑150 (MQ=255)
                                                                                                                                              atACTCCATAGCGAATACCCTGGGACAAGTAGCAGCGCTGGAATTCTCTATTCATGCCCCGAATAAAGTGCTATTAGGTGATGCTGACAGCAAATATAATCATGCCATTTTATGTCGTACCAAAGCTTATGGAACGCCAGTTAATGGCTa      >  2:3419665/2‑150 (MQ=255)
                                                                                                                                               tACTCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGTTAATGGCTAt     >  1:4106878/1‑150 (MQ=255)
                                                                                                                                                  gCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGTTAATGGCTAtttt  >  1:4418325/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TTTGTGCCATTCTGGGCTCGCTGGCGATGGGGATGATGTTTGATTCGCGTAATATCCCGGATGACTTAATGACCAACGGTCAGTATTACGCCTTTCAGAAGCTGGGCGAGTATTACAACATGGGTAATACTTTAATGGTGATTTACGCCATTGCGAATACCCTGGGACAAGTAGCGGCGCTGGTATTCTCGATTGATGCCCCGCTTAAAGTGCTATTAGGTGATGCTGACAGCAAATATATTCCAGCCAGTTTATGTCGTACCAACGCTTCTGGTACGCCCGTTAATGGCTATTTT  >  CP009273/943427‑943722

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: