Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,703,866 C→A 7.1% intergenic (‑31/+293) eptB ← / ← yhjX KDO phosphoethanolamine transferase, Ca(2+)‑inducible/pyruvate‑inducible inner membrane protein, putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,703,8660CA 7.1% 1360.8 / 50.3 366intergenic (‑31/+293)eptB/yhjXKDO phosphoethanolamine transferase, Ca(2+)‑inducible/pyruvate‑inducible inner membrane protein, putative transporter
Reads supporting (aligned to +/- strand):  ref base C (176/164);  new base A (12/14);  total (188/178)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 369 total aligned reads displayed.

AAATCGTGCGCATAGCTGCCGAAGCGGCGGTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAAT  >  CP009273/3703717‑3704012
                                                                                                                                                     |                                                                                                                                                  
aaaTCGTGCGCATAGCTGCCGAAGCGGCGGTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAac                                                                                                                                                    >  1:2522576/1‑150 (MQ=255)
     gTGCGCATAGCTGCCGAAGCGGCGGTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCt                                                                                                                                               >  2:2091606/1‑150 (MQ=255)
       gcgcATAGCTGCCGAAGCGGCGGTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGc                                                                                                                                             <  1:2428043/150‑1 (MQ=255)
          cATAGCTGCCGAAGCGGCGGTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTa                                                                                                                                          <  2:4367351/150‑1 (MQ=255)
          cATAGCTGCCGAAGCGGCGGTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTa                                                                                                                                          >  1:2600459/1‑150 (MQ=255)
           aTAGCTGCCGAAGCGGCGGTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAAAACCCTGCTTAg                                                                                                                                         <  2:3961915/150‑1 (MQ=255)
                   cGAAGCGGCGGTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCtttt                                                                                                                                 >  1:2971539/1‑150 (MQ=255)
                      aGCGGCGGTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAg                                                                                                                              >  2:4455547/1‑150 (MQ=255)
                      aGCGGCGGTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAg                                                                                                                              >  1:1958890/1‑150 (MQ=255)
                        cggcggTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGc                                                                                                                            >  2:1965874/1‑150 (MQ=255)
                        cggcggTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGc                                                                                                                            >  2:729638/1‑150 (MQ=255)
                            ggTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAAAACCCTGCTTAGTGGCTTTTCAGGCCCgg                                                                                                                        >  1:1289929/1‑150 (MQ=255)
                             gTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGc                                                                                                                       >  1:3523407/1‑150 (MQ=255)
                              tAAAAAACCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCa                                                                                                                      <  2:1158210/150‑1 (MQ=255)
                               aaaaaaCCGCGCCATTCATAAAAAGGCCAATATAGATTGCCAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAg                                                                                                                     <  2:2584096/150‑1 (MQ=255)
                               aaaaaaCCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAg                                                                                                                     <  1:2177237/150‑1 (MQ=255)
                               aaaaaaCCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAg                                                                                                                     >  1:2289286/1‑150 (MQ=255)
                                aaaaaCCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAgg                                                                                                                    >  1:1659264/1‑150 (MQ=255)
                                 aaaaCCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAggg                                                                                                                   <  2:5065503/150‑1 (MQ=255)
                                            aTTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCCTGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCgg                                                                                                        <  1:2471140/150‑1 (MQ=255)
                                             ttCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCTGa                                                                                                       >  2:3786563/1‑150 (MQ=255)
                                             ttCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGa                                                                                                       >  1:3354857/1‑150 (MQ=255)
                                               cATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAgg                                                                                                     <  1:2751856/150‑1 (MQ=255)
                                               cATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAgg                                                                                                     >  2:3527846/1‑150 (MQ=255)
                                                  aaaaaGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGa                                                                                                       <  1:5148251/145‑1 (MQ=255)
                                                  aaaaaGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAg                                                                                                  >  2:2035602/1‑150 (MQ=255)
                                                   aaaaGGCCAATATAGATTGCAAGCAAAATGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGt                                                                                                 >  2:1397298/1‑150 (MQ=255)
                                                     aacgCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCa                                                                                               <  2:4976697/147‑1 (MQ=255)
                                                     aaGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTCGTCa                                                                                               >  1:4927517/1‑150 (MQ=255)
                                                     aaGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTCGTCa                                                                                               >  1:4927427/1‑150 (MQ=255)
                                                      aGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAg                                                                                              >  2:4275300/1‑150 (MQ=255)
                                                       ggCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAAAAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGa                                                                                             >  1:1203285/1‑150 (MQ=255)
                                                         ccAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATa                                                                                           >  2:3373371/1‑150 (MQ=255)
                                                           aaTATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAg                                                                                         >  2:1206211/1‑150 (MQ=255)
                                                           aaTATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAg                                                                                         <  2:2731552/150‑1 (MQ=255)
                                                              ataGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTaa                                                                                      <  2:1571841/150‑1 (MQ=255)
                                                               taGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTaaa                                                                                     >  1:3893513/1‑150 (MQ=255)
                                                               taGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTaaa                                                                                     >  2:2616634/1‑150 (MQ=255)
                                                                 gATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAAt                                                                                   <  2:886694/150‑1 (MQ=255)
                                                                     gcaagcaaAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTac                                                                               >  2:2552617/1‑150 (MQ=255)
                                                                      caagcaaATAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAAAACCCTGCTTAGTGGCTTTTCTGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTaca                                                                              <  1:2155464/150‑1 (MQ=255)
                                                                      caagcaaAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTaca                                                                              <  1:1615788/150‑1 (MQ=255)
                                                                       aagcaaAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAAAACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTacac                                                                             >  2:1918624/1‑150 (MQ=255)
                                                                        agcaaAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACAcg                                                                            >  1:4488726/1‑150 (MQ=255)
                                                                          caaAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACAcgcg                                                                          >  2:2943828/1‑150 (MQ=255)
                                                                          caaAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACAcgcg                                                                          >  1:2904141/1‑150 (MQ=255)
                                                                           aaaaaGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGACAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACAcgcgc                                                                         >  2:2377200/1‑150 (MQ=255)
                                                                            aaaaGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGAGAAGTAAAATCTACAcgcgcg                                                                        >  2:3086214/1‑150 (MQ=255)
                                                                              aaGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAAAACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGaa                                                                      >  2:1743397/1‑150 (MQ=255)
                                                                               aGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAAt                                                                     <  2:79923/150‑1 (MQ=255)
                                                                               aGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAAt                                                                     <  2:952868/150‑1 (MQ=255)
                                                                                gCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATa                                                                    >  2:5236731/1‑150 (MQ=255)
                                                                                 cTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATaa                                                                   <  2:5165500/150‑1 (MQ=255)
                                                                                  tCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCCTGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAg                                                                  <  1:1112559/150‑1 (MQ=255)
                                                                                  tCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAg                                                                  <  1:1112848/150‑1 (MQ=255)
                                                                                  tCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAg                                                                  <  1:3421561/150‑1 (MQ=255)
                                                                                    aGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGACTAAGTa                                                                >  2:3254142/1‑150 (MQ=255)
                                                                                    aGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTa                                                                >  2:618654/1‑150 (MQ=255)
                                                                                    aGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAAAACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTa                                                                >  2:1692542/1‑150 (MQ=255)
                                                                                     gCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAg                                                               <  1:824573/150‑1 (MQ=255)
                                                                                        tCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGTTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGAcc                                                            <  2:207719/150‑1 (MQ=255)
                                                                                        tCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGAcc                                                            <  2:207911/150‑1 (MQ=255)
                                                                                            ctgTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAg                                                        >  1:549025/1‑150 (MQ=255)
                                                                                             tgtgTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGc                                                       >  1:1922089/1‑150 (MQ=255)
                                                                                                 tAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAg                                                   <  1:980568/150‑1 (MQ=255)
                                                                                                 tAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAg                                                   >  2:3870184/1‑150 (MQ=255)
                                                                                                  aaTCTATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACAACCTGCGTATTGGCTTTTCAGTCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTTAAATCTACACGCGCGAATAAGTAGACCACAGCAAAgg                                                  <  2:4516369/150‑1 (MQ=255)
                                                                                                  aaTCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTCCACGCGCGAATAAGTAGACCACAGCAAAgg                                                  >  1:3639986/1‑150 (MQ=255)
                                                                                                    tCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAAAACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGa                                                <  1:5203937/150‑1 (MQ=255)
                                                                                                      gATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGaaa                                              <  1:806177/150‑1 (MQ=255)
                                                                                                        tttGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGaaaaa                                            <  2:4828766/150‑1 (MQ=255)
                                                                                                           gATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAgtg                                         >  1:136947/1‑150 (MQ=255)
                                                                                                            aTGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGACAAAgtgt                                        >  2:4629724/1‑150 (MQ=255)
                                                                                                            aTGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAgtgt                                        <  1:4375962/150‑1 (MQ=255)
                                                                                                              gTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAAAACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTgg                                      >  1:1821571/1‑150 (MQ=255)
                                                                                                                  ctcATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGATTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGaa                                  <  2:4033137/150‑1 (MQ=255)
                                                                                                                   tcATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGaaa                                 >  2:3137400/1‑150 (MQ=255)
                                                                                                                    cATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTCCACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGaatg                                >  1:2659706/1‑148 (MQ=255)
                                                                                                                     aTGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAaga                               >  1:4761493/1‑150 (MQ=255)
                                                                                                                       gCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAgt                             <  1:5004548/150‑1 (MQ=255)
                                                                                                                        caaacaaaCCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAgtg                            <  1:5065159/150‑1 (MQ=255)
                                                                                                                        caaacaaaCCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAgtg                            <  1:3065462/150‑1 (MQ=255)
                                                                                                                         aaacaaacCCAGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAgtgt                           <  2:865605/150‑1 (MQ=255)
                                                                                                                          aacaaacCCTGGACAAACAGGTGATAAAACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGtgtg                          <  1:2087264/150‑1 (MQ=255)
                                                                                                                            caaacCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCa                        <  1:672950/150‑1 (MQ=255)
                                                                                                                            caaacCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCa                        <  1:2810030/150‑1 (MQ=255)
                                                                                                                             aaacCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCAt                       <  1:4155498/150‑1 (MQ=255)
                                                                                                                                cccTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGaa                    >  2:3589127/1‑150 (MQ=255)
                                                                                                                                  cTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAATAGTGTGGGGAAAGAGTGTGCATGAAGc                  >  2:4122158/1‑150 (MQ=255)
                                                                                                                                  cTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGc                  >  1:910578/1‑150 (MQ=255)
                                                                                                                                     gACAAACCGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGcaca               <  2:2838569/150‑1 (MQ=255)
                                                                                                                                     gACAAACAGGTGATAAAACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGcaca               <  1:3257438/150‑1 (MQ=255)
                                                                                                                                       cAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACaaa             >  1:3681926/1‑150 (MQ=255)
                                                                                                                                         aaCAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGACAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAAtt           >  1:3229164/1‑150 (MQ=255)
                                                                                                                                          aCAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAAttt          <  2:3348798/150‑1 (MQ=255)
                                                                                                                                            aGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACCCGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTAc        >  1:4644667/1‑150 (MQ=255)
                                                                                                                                             ggTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACa       <  1:4064992/150‑1 (MQ=255)
                                                                                                                                              gTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCACGAAGCACAAATTTACaa      <  1:1323496/150‑1 (MQ=255)
                                                                                                                                                 aTAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACaaaaa   <  2:2684443/150‑1 (MQ=255)
                                                                                                                                                  tAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAAt  >  2:4232687/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AAATCGTGCGCATAGCTGCCGAAGCGGCGGTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAGATTGCAAGCAAAAAGCTCAGCTTCTGCTGTGTAATCGATTTGATGTATCTCATGCAAACAAACCCTGGACAAACAGGTGATAACACCCTGCTTAGTGGCTTTTCAGGCCCGGCAGGGAAAAAGTGCGGAGGTAGTCAGATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAAT  >  CP009273/3703717‑3704012

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: