Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,007,147 A→C 5.4% E245D (GAA→GAC metE → 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,007,1470AC 5.4% 1096.0 / 43.6 336E245D (GAA→GACmetE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (182/135);  new base C (8/10);  total (190/146)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 345 total aligned reads displayed.

GCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATGGTCGTAA  >  CP009273/4006998‑4007296
                                                                                                                                                     |                                                                                                                                                     
  aGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAgg                                                                                                                                                     >  2:4565214/1‑150 (MQ=255)
  aGAACTGGCGAAACGCGGCATCGAGCGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAgg                                                                                                                                                     >  2:4506089/1‑150 (MQ=255)
   gAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGc                                                                                                                                                    >  1:743652/1‑150 (MQ=255)
     aCTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGt                                                                                                                                                  >  2:1481718/1‑150 (MQ=255)
        ggCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAAc                                                                                                                                               <  1:1668924/150‑1 (MQ=255)
         gCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACg                                                                                                                                              <  2:3280303/150‑1 (MQ=255)
          cGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGc                                                                                                                                             >  1:1932801/1‑150 (MQ=255)
           gAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGcc                                                                                                                                            >  2:2240459/1‑150 (MQ=255)
           gAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGcc                                                                                                                                            >  2:3560949/1‑150 (MQ=255)
            aaaCGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCa                                                                                                                                           >  2:3010019/1‑150 (MQ=255)
             aaCGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCaa                                                                                                                                          >  1:637212/1‑150 (MQ=255)
             aaCGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCaa                                                                                                                                          >  1:2763125/1‑150 (MQ=255)
                 cgGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAAtct                                                                                                                                      <  1:1058884/150‑1 (MQ=255)
                    cATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGACGGCGTAACGCCAAATCTCGa                                                                                                                                   <  2:3374073/150‑1 (MQ=255)
                        gAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACg                                                                                                                               >  2:904739/1‑150 (MQ=255)
                           tGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGAtt                                                                                                                            <  2:2000510/150‑1 (MQ=255)
                            gggTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTa                                                                                                                           <  2:5437402/150‑1 (MQ=255)
                              gTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGACGGCGTAACGCCAAATCTCGACACGATTACt                                                                                                                         <  1:2909903/150‑1 (MQ=255)
                               tACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTg                                                                                                                        <  2:863986/150‑1 (MQ=255)
                                 cAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTgcg                                                                                                                      >  1:3360260/1‑150 (MQ=255)
                                  aGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTgcgc                                                                                                                     >  2:4178698/1‑150 (MQ=255)
                                   gATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCt                                                                                                                    >  2:721684/1‑150 (MQ=255)
                                        atgaACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTg                                                                                                               >  1:3926829/1‑150 (MQ=255)
                                        atgaACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTg                                                                                                               >  1:4904602/1‑150 (MQ=255)
                                        atgaACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATATTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTg                                                                                                               >  2:4378570/1‑150 (MQ=255)
                                          gaACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGtt                                                                                                             <  1:118395/150‑1 (MQ=255)
                                            aCCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCa                                                                                                           >  2:4954285/1‑150 (MQ=255)
                                            aCCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCa                                                                                                           >  2:3838997/1‑150 (MQ=255)
                                             cccGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAg                                                                                                          >  1:4618900/1‑150 (MQ=255)
                                             cccGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAg                                                                                                          >  2:3199831/1‑150 (MQ=255)
                                             cccGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAg                                                                                                          >  2:1269860/1‑150 (MQ=255)
                                              ccGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAgg                                                                                                         <  1:3668177/150‑1 (MQ=255)
                                               cgcgCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAggg                                                                                                        <  1:4726216/150‑1 (MQ=255)
                                               cgcgCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAggg                                                                                                        >  1:419482/1‑150 (MQ=255)
                                               cgcgCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAggg                                                                                                        >  1:50637/1‑150 (MQ=255)
                                                  gcTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCt                                                                                                     <  1:1317579/150‑1 (MQ=255)
                                                  gcTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCt                                                                                                     <  1:1481710/150‑1 (MQ=255)
                                                   cTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTg                                                                                                    >  1:5132798/1‑150 (MQ=255)
                                                   cTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTg                                                                                                    >  1:531165/1‑150 (MQ=255)
                                                    tGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGc                                                                                                   >  2:5284425/1‑150 (MQ=255)
                                                       tACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATg                                                                                                >  2:584459/1‑150 (MQ=255)
                                                       tACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATg                                                                                                <  1:5578692/150‑1 (MQ=255)
                                                        aCTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGt                                                                                               <  1:3160042/150‑1 (MQ=255)
                                                         cTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGACGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGtt                                                                                              <  2:4858563/150‑1 (MQ=255)
                                                         cTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGtt                                                                                              >  1:4241932/1‑150 (MQ=255)
                                                          tGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTg                                                                                             <  1:3401795/150‑1 (MQ=255)
                                                               cTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGACGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTc                                                                                        >  2:3681893/1‑150 (MQ=255)
                                                               cTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTc                                                                                        >  1:816963/1‑150 (MQ=255)
                                                                 aCCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGt                                                                                      <  1:2725370/150‑1 (MQ=255)
                                                                 aCCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGt                                                                                      <  1:4738171/150‑1 (MQ=255)
                                                                  ccACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTa                                                                                     <  2:4226222/150‑1 (MQ=255)
                                                                    acaGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACa                                                                                   >  2:3111143/1‑150 (MQ=255)
                                                                          gTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTaa                                                                             >  2:4806440/1‑150 (MQ=255)
                                                                          gTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTaa                                                                             >  2:4252917/1‑150 (MQ=255)
                                                                          gTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTaa                                                                             >  1:2371033/1‑150 (MQ=255)
                                                                            ggCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAg                                                                           <  2:2073525/150‑1 (MQ=255)
                                                                            ggCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAg                                                                           <  2:4908389/150‑1 (MQ=255)
                                                                               tGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATg                                                                        <  2:1346918/150‑1 (MQ=255)
                                                                                ggACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGa                                                                       <  2:5135828/150‑1 (MQ=255)
                                                                                    gCATACAAACCCTCTTAGGACGCGCTCCATGGACAGGTTAAACTGCTGCTGACCACCTATTTTGAAGGCGTAAGGCCAATTCTCGTCCCGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTAGATGGTAAAGATGACGtt                                                                   <  2:4674735/150‑1 (MQ=255)
                                                                                    gCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGtt                                                                   >  2:284863/1‑150 (MQ=255)
                                                                                    gCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGtt                                                                   <  2:4062067/150‑1 (MQ=255)
                                                                                     cATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGACGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTg                                                                  >  2:2800673/1‑150 (MQ=255)
                                                                                        aCAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTg                                                               >  2:4415259/1‑150 (MQ=255)
                                                                                          aaaCCCTCTTAATACGCGCTCCATGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCATATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTTCATGTTGACCTCGTTCATGGTAAAGATGACGTTGCTGaa                                                             <  2:1725104/150‑1 (MQ=255)
                                                                                          aaaCCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGaa                                                             <  2:1752573/150‑1 (MQ=255)
                                                                                           aaCCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAAc                                                            <  2:2539544/150‑1 (MQ=255)
                                                                                            aCCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACCCGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACt                                                           >  2:2734715/1‑150 (MQ=255)
                                                                                             cccGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTg                                                          >  1:2627162/1‑150 (MQ=255)
                                                                                                    acgacgCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgc                                                   >  1:396015/1‑150 (MQ=255)
                                                                                                      gacgCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgcgc                                                 >  2:5365099/1‑150 (MQ=255)
                                                                                                             tCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctgccta                                          >  2:3031398/1‑149 (MQ=255)
                                                                                                             tCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCtt                                          <  2:1438025/150‑1 (MQ=255)
                                                                                                             tCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCtt                                          <  2:743663/150‑1 (MQ=255)
                                                                                                             tCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCtt                                          >  2:5242989/1‑150 (MQ=255)
                                                                                                              ccAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTc                                         <  1:4042773/150‑1 (MQ=255)
                                                                                                                aGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTg                                       >  1:3582337/1‑150 (MQ=255)
                                                                                                                aGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTg                                       >  1:1836435/1‑150 (MQ=255)
                                                                                                                   gACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACt                                    <  1:3199831/150‑1 (MQ=255)
                                                                                                                     cAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTgg                                  <  2:1042949/150‑1 (MQ=255)
                                                                                                                       ggTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGtt                                >  2:1809490/1‑150 (MQ=255)
                                                                                                                       ggTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGtt                                <  2:298969/150‑1 (MQ=255)
                                                                                                                        gTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTg                               >  1:4368981/1‑150 (MQ=255)
                                                                                                                           aaaCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTg                            <  1:1247340/150‑1 (MQ=255)
                                                                                                                            aaCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGt                           >  1:4478480/1‑150 (MQ=255)
                                                                                                                             aCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTc                          <  2:1901352/150‑1 (MQ=255)
                                                                                                                                gctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGc                       >  1:2284737/1‑150 (MQ=255)
                                                                                                                                    ctgACCACCTATTTTGACGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGt                   >  1:4436828/1‑150 (MQ=255)
                                                                                                                                    ctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGt                   <  2:1213893/150‑1 (MQ=255)
                                                                                                                                    ctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGt                   >  1:588418/1‑150 (MQ=255)
                                                                                                                                      gACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCt                 >  2:2662228/1‑150 (MQ=255)
                                                                                                                                      gACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCt                 >  2:5505779/1‑150 (MQ=255)
                                                                                                                                         caccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGAt              <  1:904739/150‑1 (MQ=255)
                                                                                                                                          accTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATc             >  1:2480398/1‑150 (MQ=255)
                                                                                                                                            cTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCaa           >  2:2385389/1‑150 (MQ=255)
                                                                                                                                             tATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAAt          <  1:4539561/150‑1 (MQ=255)
                                                                                                                                              aatattcaGGCGTAACGCCAAATCTCGACCCGATTACGGCGCTGCCTGTTCAGGGTCTGCAGGTTGACCTCGTACATGGTAACGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATg         <  2:2877765/143‑1 (MQ=255)
                                                                                                                                               ttttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATgg        <  1:2065349/150‑1 (MQ=255)
                                                                                                                                                 ttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATGGTc      <  2:2191383/150‑1 (MQ=255)
                                                                                                                                                  tGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATGGTCg     <  2:3613583/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACGACGCGCTCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATGGTCGTAA  >  CP009273/4006998‑4007296

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: