Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 231,378 C→A 100% R26S (CGC→AGC)  yafS → putative S‑adenosyl‑L‑methionine‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273231,3780CA98.3% 962.1 / 5.5 299R26S (CGC→AGC) yafSputative S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (5/0);  new base A (159/135);  total (164/135)
Fisher's exact test for biased strand distribution p-value = 6.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 300 total aligned reads displayed.

TAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCC  >  CP009273/231232‑231523
                                                                                                                                                  |                                                                                                                                                 
tAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGc                                                                                                                                                 >  2:3608539/1‑149 (MQ=255)
 aaTTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCa                                                                                                                                                >  2:1116315/1‑149 (MQ=255)
 aaTTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCa                                                                                                                                                >  2:1895414/1‑149 (MQ=255)
 aaTTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCa                                                                                                                                                >  2:1116423/1‑149 (MQ=255)
  aTTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCaa                                                                                                                                               <  1:2606239/149‑1 (MQ=255)
     gTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGgc                                                                                                                                            >  1:2994883/1‑149 (MQ=255)
      tcatcaAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGgcg                                                                                                                                           >  2:1526174/1‑149 (MQ=255)
       catcaAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGgcgc                                                                                                                                          <  2:1074237/149‑1 (MQ=255)
        atcaAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCt                                                                                                                                         <  1:1151313/149‑1 (MQ=255)
         tcaAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTg                                                                                                                                        >  1:911569/1‑149 (MQ=255)
           aaaGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGa                                                                                                                                      <  2:708000/149‑1 (MQ=255)
             aGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAgc                                                                                                                                    <  2:1242215/149‑1 (MQ=255)
                cGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGcc                                                                                                                                 <  1:2069217/149‑1 (MQ=255)
                cGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGcc                                                                                                                                 <  1:2312898/149‑1 (MQ=255)
                cGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGcc                                                                                                                                 <  1:238524/149‑1 (MQ=255)
                 gggAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCa                                                                                                                                >  1:3191323/1‑149 (MQ=255)
                   gAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGc                                                                                                                              <  2:2164790/149‑1 (MQ=255)
                   gAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGc                                                                                                                              >  2:2248829/1‑149 (MQ=255)
                     aTACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTc                                                                                                                            >  1:454229/1‑149 (MQ=255)
                      tACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCa                                                                                                                           >  2:1182102/1‑149 (MQ=255)
                      tACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCa                                                                                                                           >  2:2686050/1‑149 (MQ=255)
                       aCTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCaa                                                                                                                          <  2:395818/149‑1 (MQ=255)
                        cTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAAc                                                                                                                         >  2:94206/1‑149 (MQ=255)
                            tAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGt                                                                                                                     <  1:3110295/149‑1 (MQ=255)
                            tAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGt                                                                                                                     <  1:3418182/149‑1 (MQ=255)
                             aaGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTg                                                                                                                    >  1:1568661/1‑149 (MQ=255)
                             aaGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTg                                                                                                                    >  1:2545152/1‑149 (MQ=255)
                             aaGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTg                                                                                                                    >  1:642070/1‑149 (MQ=255)
                               gATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGt                                                                                                                  >  1:18430/1‑149 (MQ=255)
                               gATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGt                                                                                                                  <  1:3811138/149‑1 (MQ=255)
                               gATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGt                                                                                                                  <  1:3628766/149‑1 (MQ=255)
                               gATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGt                                                                                                                  <  1:3760572/149‑1 (MQ=255)
                                aTTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGtt                                                                                                                 <  2:1853633/149‑1 (MQ=255)
                                   cATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCAc                                                                                                              <  1:2811185/149‑1 (MQ=255)
                                    aTAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACt                                                                                                             >  1:1087702/1‑149 (MQ=255)
                                    aTAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACt                                                                                                             >  1:4109026/1‑149 (MQ=255)
                                    aTAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACt                                                                                                             <  2:2653785/149‑1 (MQ=255)
                                      aGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTaa                                                                                                           >  2:3103015/1‑149 (MQ=255)
                                         ttACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAAt                                                                                                        >  2:476360/1‑149 (MQ=255)
                                         ttACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAACGCGCTGTAGCGCCAGCTCAACCCGTGGTTAACTAAAAt                                                                                                        >  1:530100/1‑149 (MQ=255)
                                           aCCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGt                                                                                                      >  2:881883/1‑149 (MQ=255)
                                            ccTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTa                                                                                                     >  2:866043/1‑149 (MQ=255)
                                                tcAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGt                                                                                                 >  1:2331382/1‑149 (MQ=255)
                                                tcAGTGAGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGt                                                                                                 <  2:245079/149‑1 (MQ=255)
                                                   gtgtGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGtttt                                                                                              >  2:3417927/1‑149 (MQ=255)
                                                     gtgAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCa                                                                                            <  2:1994241/149‑1 (MQ=255)
                                                            ggAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCtt                                                                                     >  1:2283379/1‑149 (MQ=255)
                                                             gAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTa                                                                                    >  1:1284283/1‑149 (MQ=255)
                                                              aaGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTaa                                                                                   >  1:1138285/1‑149 (MQ=255)
                                                               aGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAg                                                                                  >  1:3231806/1‑149 (MQ=255)
                                                                  ttGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGAtt                                                                               <  2:3748016/149‑1 (MQ=255)
                                                                    gtgATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTgg                                                                             <  2:2134777/149‑1 (MQ=255)
                                                                    gtgATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTgg                                                                             <  1:778389/149‑1 (MQ=255)
                                                                      gatgaAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCa                                                                           >  1:4082732/1‑149 (MQ=255)
                                                                      gatgaAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCa                                                                           >  2:2820273/1‑149 (MQ=255)
                                                                        tgaAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAt                                                                         >  1:980427/1‑149 (MQ=255)
                                                                         gaAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAtt                                                                        >  2:3851660/1‑149 (MQ=255)
                                                                         gaAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAtt                                                                        <  1:3926233/149‑1 (MQ=255)
                                                                         gaAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAtt                                                                        <  2:915036/149‑1 (MQ=255)
                                                                          aaaCCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAttt                                                                       >  1:1280066/1‑149 (MQ=255)
                                                                            aCCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTaa                                                                     >  1:3769485/1‑149 (MQ=255)
                                                                             ccGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGGTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAg                                                                    <  2:174441/149‑1 (MQ=255)
                                                                                   agagTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGa                                                              >  2:137568/1‑149 (MQ=255)
                                                                                      gTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAAt                                                           >  1:13818/1‑149 (MQ=255)
                                                                                          cTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAAt                                                       <  2:1340388/149‑1 (MQ=255)
                                                                                           tCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAAtt                                                      <  2:1280066/149‑1 (MQ=255)
                                                                                              aaCTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCg                                                   >  2:1418745/1‑149 (MQ=255)
                                                                                               aCTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGa                                                  <  1:1989656/149‑1 (MQ=255)
                                                                                                  gTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGc                                               <  2:3937866/149‑1 (MQ=255)
                                                                                                   tagaggCACCTGAATGCAGGCGCGATATGCCCTGGTGAAAGCTATATAGCAAGTCGATGGAGCTCCAGATCAACCCGTGGTTCACTAAAATGTATGTTTTTCATCTGCTAAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCg                                              >  2:1372809/5‑149 (MQ=255)
                                                                                                    cGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGt                                             <  2:1949427/149‑1 (MQ=255)
                                                                                                      tGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTgc                                           <  2:1685508/149‑1 (MQ=255)
                                                                                                      tGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTgc                                           <  1:1484519/149‑1 (MQ=255)
                                                                                                       ggCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTgcg                                          >  1:2681669/1‑149 (MQ=255)
                                                                                                       ggCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTgcg                                          >  1:1289611/1‑149 (MQ=255)
                                                                                                          tCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCgg                                       <  1:837081/149‑1 (MQ=255)
                                                                                                          tCCTGAGTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTAGAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCgg                                       <  2:3651137/149‑1 (MQ=255)
                                                                                                           ccTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                                      >  2:2573431/1‑149 (MQ=255)
                                                                                                           ccTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                                      <  1:1638995/149‑1 (MQ=255)
                                                                                                           ccTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                                      <  2:18430/149‑1 (MQ=255)
                                                                                                           ccTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                                      >  2:506924/1‑149 (MQ=255)
                                                                                                             tGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGttt                                    >  2:1485936/1‑149 (MQ=255)
                                                                                                                   cTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCtcatc                              <  2:3498820/149‑1 (MQ=255)
                                                                                                                     ggggCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCaa                            <  2:2892359/149‑1 (MQ=255)
                                                                                                                      gggCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAg                           >  2:3352688/1‑149 (MQ=255)
                                                                                                                      gggCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAg                           >  1:597277/1‑149 (MQ=255)
                                                                                                                          gATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGaa                       >  2:589321/1‑149 (MQ=255)
                                                                                                                           aTTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAAt                      <  2:3206485/149‑1 (MQ=255)
                                                                                                                           aTTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAAt                      <  2:2994883/149‑1 (MQ=255)
                                                                                                                           aTTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAAt                      >  1:3237047/1‑149 (MQ=255)
                                                                                                                                 ccTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCt                <  2:3137092/149‑1 (MQ=255)
                                                                                                                                   tGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTgc              <  1:3589760/149‑1 (MQ=255)
                                                                                                                                   tGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTGTCTCGTCAAGTGAATGTTTCTgc              >  1:1206045/1‑149 (MQ=255)
                                                                                                                                       gAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCaa          >  1:337537/1‑149 (MQ=255)
                                                                                                                                       gAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCaa          >  1:2724195/1‑149 (MQ=255)
                                                                                                                                          aGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGa       >  1:1773079/1‑149 (MQ=255)
                                                                                                                                          aGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGa       <  1:1887668/149‑1 (MQ=255)
                                                                                                                                           gCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGaa      >  2:2966032/1‑149 (MQ=255)
                                                                                                                                           gCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGaa      >  1:2837143/1‑149 (MQ=255)
                                                                                                                                              ttATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGc   <  2:2724195/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCC  >  CP009273/231232‑231523

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: ‑