Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,605,290 A→T 11.1% L24Q (CTG→CAG)  yjjW ← putative pyruvate formate lyase activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,605,2900AT11.1% 843.7 / 77.9 260L24Q (CTG→CAG) yjjWputative pyruvate formate lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base A (106/125);  new base T (16/13);  total (122/138)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 268 total aligned reads displayed.

CAGCGTTCCACACCACTTTGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCG  >  CP009273/4605141‑4605439
                                                                                                                                                     |                                                                                                                                                     
 aGCGTTCCACACCACTTTGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAg                                                                                                                                                      >  2:2983312/1‑150 (MQ=255)
 aGCGTTCCACACCACTTTGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAg                                                                                                                                                      >  1:3150677/1‑150 (MQ=255)
      tCCACACCACTTTGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGAGTGCAGCCCTGCAAAAACAGAGCCAgacgac                                                                                                                                                 >  2:3041991/1‑150 (MQ=255)
        cccaccacTTTGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCCCAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTg                                                                                                                                               <  1:1668914/148‑1 (MQ=255)
            ccacTTTGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTg                                                                                                                                           <  1:4073112/150‑1 (MQ=255)
            ccacTTTGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTg                                                                                                                                           <  1:2034423/150‑1 (MQ=255)
              acTTTGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGc                                                                                                                                         <  2:4514872/150‑1 (MQ=255)
               cTTTGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGcc                                                                                                                                        <  1:3856022/150‑1 (MQ=255)
               cTTTGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGcc                                                                                                                                        <  2:651396/150‑1 (MQ=255)
                tttGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCTGACGACTGCCTGGccc                                                                                                                                       >  1:5121312/1‑150 (MQ=255)
                 ttGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCg                                                                                                                                      >  1:4963560/1‑150 (MQ=255)
                 ttGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCg                                                                                                                                      >  2:4712004/1‑150 (MQ=255)
                   gCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTc                                                                                                                                    >  1:3360666/1‑150 (MQ=255)
                      gTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCCCAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAAc                                                                                                                                 <  2:711547/150‑1 (MQ=255)
                      gTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAAc                                                                                                                                 >  2:1651156/1‑150 (MQ=255)
                           caaTCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAgg                                                                                                                            >  2:873362/1‑150 (MQ=255)
                            aaTCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGga                                                                                                                           <  1:2504802/150‑1 (MQ=255)
                               cTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGaa                                                                                                                        <  1:794807/150‑1 (MQ=255)
                                tGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAg                                                                                                                       <  2:2869215/150‑1 (MQ=255)
                                tGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAg                                                                                                                       <  2:2538193/150‑1 (MQ=255)
                                  cAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAggg                                                                                                                     <  1:3301117/150‑1 (MQ=255)
                                   aaCGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAgggg                                                                                                                    >  2:3995890/1‑150 (MQ=255)
                                   aaCGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAgggg                                                                                                                    <  1:3510524/150‑1 (MQ=255)
                                   aaCGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAgggg                                                                                                                    >  1:4594011/1‑150 (MQ=255)
                                    aCGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGa                                                                                                                   <  2:3161621/150‑1 (MQ=255)
                                      gCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATa                                                                                                                 <  1:2582296/150‑1 (MQ=255)
                                      gCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATa                                                                                                                 >  1:4211664/1‑150 (MQ=255)
                                       ccTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCTGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATaa                                                                                                                <  1:1080107/150‑1 (MQ=255)
                                       ccTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATaa                                                                                                                <  1:4033186/150‑1 (MQ=255)
                                        cTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAAt                                                                                                               <  2:5162890/150‑1 (MQ=255)
                                         tGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATc                                                                                                              >  2:4617333/1‑150 (MQ=255)
                                           aTGCGGACACTGTGGCACGCACTCCCCACCGTCATAGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCtt                                                                                                            <  2:5095892/150‑1 (MQ=255)
                                               ggACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTg                                                                                                        <  1:2958091/150‑1 (MQ=255)
                                                  cacTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACt                                                                                                     <  2:2691874/150‑1 (MQ=255)
                                                  cacTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACt                                                                                                     <  2:344953/150‑1 (MQ=255)
                                                    cTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCTGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTaa                                                                                                   <  1:3463529/150‑1 (MQ=255)
                                                      gtgGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAg                                                                                                 <  1:2668291/150‑1 (MQ=255)
                                                       tgGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAgc                                                                                                >  1:2403504/1‑150 (MQ=255)
                                                       tgGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAgc                                                                                                <  2:1601438/150‑1 (MQ=255)
                                                        ggCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAgcg                                                                                               >  2:43256/1‑150 (MQ=255)
                                                           acgcacTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCAt                                                                                            <  1:1970182/150‑1 (MQ=255)
                                                            cgcacTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATc                                                                                           >  2:3346612/1‑150 (MQ=255)
                                                             gcacTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCTGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATct                                                                                          >  1:3789219/1‑150 (MQ=255)
                                                             gcacTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCTGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATct                                                                                          >  1:1192053/1‑150 (MQ=255)
                                                               acTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCTGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATctgc                                                                                        >  1:4966557/1‑150 (MQ=255)
                                                               acTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCTGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCATACTGACTAAAGCGCATctgc                                                                                        >  2:2888563/1‑150 (MQ=255)
                                                                cTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATctgct                                                                                       >  2:1892496/1‑150 (MQ=255)
                                                                cTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATctgct                                                                                       >  2:4342606/1‑150 (MQ=255)
                                                                 tCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATctgctg                                                                                      <  1:2795272/150‑1 (MQ=255)
                                                                   cccACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGtt                                                                                    >  1:2930207/1‑150 (MQ=255)
                                                                       caGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATg                                                                                >  1:3683716/1‑150 (MQ=255)
                                                                        aGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATgg                                                                               <  2:2887000/150‑1 (MQ=255)
                                                                        aGTCAGTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATgg                                                                               <  2:3536752/150‑1 (MQ=255)
                                                                           cATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCtt                                                                            <  2:3468252/150‑1 (MQ=255)
                                                                           cATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCtt                                                                            <  2:2326164/150‑1 (MQ=255)
                                                                             ttGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTAt                                                                          <  1:552085/150‑1 (MQ=255)
                                                                              tGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATc                                                                         >  1:1139063/1‑150 (MQ=255)
                                                                               gCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCa                                                                        <  1:1006129/150‑1 (MQ=255)
                                                                               gCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCa                                                                        >  2:3832093/1‑150 (MQ=255)
                                                                                cAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCAc                                                                       <  2:1447784/150‑1 (MQ=255)
                                                                                  aCGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCAcgc                                                                     >  2:5089612/1‑150 (MQ=255)
                                                                                   cGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCAcgcg                                                                    <  2:3166370/150‑1 (MQ=255)
                                                                                     tCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCTGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCAcgcgcg                                                                  <  1:891358/150‑1 (MQ=255)
                                                                                     tCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCAcgcgcg                                                                  >  1:4384138/1‑150 (MQ=255)
                                                                                         aTCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTg                                                              >  2:1955101/1‑150 (MQ=255)
                                                                                           cGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGc                                                            <  2:98087/150‑1 (MQ=255)
                                                                                                 cGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCa                                                      <  2:171039/150‑1 (MQ=255)
                                                                                                     ttGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAAt                                                  <  1:3582903/150‑1 (MQ=255)
                                                                                                        tgACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACg                                               <  2:273001/150‑1 (MQ=255)
                                                                                                         gACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCAGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGa                                              >  1:4322707/1‑150 (MQ=255)
                                                                                                          aCAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAg                                             >  2:4128313/1‑150 (MQ=255)
                                                                                                              tttttGCAGCGCAGAGTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGtt                                         <  1:3521908/150‑1 (MQ=255)
                                                                                                               ttttGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTtg                                        >  1:4303357/1‑150 (MQ=255)
                                                                                                                 ttGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTtgcg                                      <  1:4746088/150‑1 (MQ=255)
                                                                                                                     aGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTtgcgtgcg                                  <  2:1250285/150‑1 (MQ=255)
                                                                                                                       cgcAGATTGCAGCCCTGCAAAAACAGAGCCTGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGc                                <  1:3113200/150‑1 (MQ=255)
                                                                                                                        gcAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGct                               >  2:1135229/1‑150 (MQ=255)
                                                                                                                         cAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGctt                              <  1:3728017/150‑1 (MQ=255)
                                                                                                                           gATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGcttct                            <  2:4146795/150‑1 (MQ=255)
                                                                                                                            aTTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGcttctt                           <  1:3780483/150‑1 (MQ=255)
                                                                                                                            aTTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGcttctt                           <  1:1183482/150‑1 (MQ=255)
                                                                                                                             ttGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCttcttc                          <  1:3970991/150‑1 (MQ=255)
                                                                                                                              tGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCg                         >  1:1444370/1‑150 (MQ=255)
                                                                                                                                 aGCCCTGCAAAAACAGAGCCTGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGccc                      <  2:4966536/150‑1 (MQ=255)
                                                                                                                                  gCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCcca                     >  1:4653024/1‑150 (MQ=255)
                                                                                                                                  gCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCcca                     <  2:1106796/150‑1 (MQ=255)
                                                                                                                                    ccTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCccagc                   >  2:3902206/1‑150 (MQ=255)
                                                                                                                                     cTGCAAAAACAGAGCCTGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCccagcc                  >  1:5057300/1‑150 (MQ=255)
                                                                                                                                     cTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCccagcc                  >  2:1167155/1‑150 (MQ=255)
                                                                                                                                     cTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCccagcc                  >  2:171095/1‑150 (MQ=255)
                                                                                                                                       gCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCccagccag                >  1:1920641/1‑150 (MQ=255)
                                                                                                                                       gCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCccagccag                >  1:2186848/1‑150 (MQ=255)
                                                                                                                                       gCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCccagccag                <  2:2065960/150‑1 (MQ=255)
                                                                                                                                           aaaCAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCCCAGCCAggtgg            <  2:2393826/150‑1 (MQ=255)
                                                                                                                                            aaCAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCCCAGCCAggtggt           <  2:1500877/150‑1 (MQ=255)
                                                                                                                                            aaCAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCCCAGCCAggtggt           <  1:5179589/150‑1 (MQ=255)
                                                                                                                                               agagCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCCCAGCCAGGTGGTGtt        >  1:2234336/1‑150 (MQ=255)
                                                                                                                                                gagCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCCCAGCCAGGTGGTGTTg       >  1:166960/1‑150 (MQ=255)
                                                                                                                                                     aGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCg  >  2:4189700/1‑150 (MQ=255)
                                                                                                                                                     aGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCg  >  2:1145959/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCGTTCCACACCACTTTGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTGCGGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCG  >  CP009273/4605141‑4605439

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: