Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,138,702 T→A 25.7% S652R (AGT→AGA yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,138,7020TA25.7% 578.2 / 215.1 310S652R (AGT→AGAyegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (121/109);  new base A (44/36);  total (165/145)
Fisher's exact test for biased strand distribution p-value = 7.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 314 total aligned reads displayed.

GTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTAT  >  CP009273/2138553‑2138851
                                                                                                                                                     |                                                                                                                                                     
atgtTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGt                                                                                                                                                       >  1:4488014/2‑150 (MQ=255)
   ttGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTAct                                                                                                                                                    >  2:1699034/1‑150 (MQ=255)
    tGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTActc                                                                                                                                                   <  2:565632/150‑1 (MQ=255)
     gCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAActct                                                                                                                                                  <  2:2094558/150‑1 (MQ=255)
         atTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGAc                                                                                                                                              <  1:2079022/150‑1 (MQ=255)
            aCTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGgc                                                                                                                                           >  2:2411714/1‑150 (MQ=255)
                ttGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGgcagca                                                                                                                                       <  2:4983124/150‑1 (MQ=255)
                ttGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGgcagca                                                                                                                                       <  1:1784770/150‑1 (MQ=255)
                     gACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTCAGTACTCTGGACGGCAGCAATAtt                                                                                                                                  >  2:343769/1‑150 (MQ=255)
                        aaCGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTggt                                                                                                                               <  1:5255593/150‑1 (MQ=255)
                        aaCGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTggt                                                                                                                               >  2:1085707/1‑150 (MQ=255)
                         aCGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTggtt                                                                                                                              <  2:1438282/150‑1 (MQ=255)
                          cGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTggttc                                                                                                                             >  1:107956/1‑150 (MQ=255)
                          cGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTggttc                                                                                                                             <  2:3566925/150‑1 (MQ=255)
                           ggACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTggttcg                                                                                                                            >  1:5318295/1‑150 (MQ=255)
                           ggACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTggttcg                                                                                                                            >  1:5125018/1‑150 (MQ=255)
                           ggACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTggttcg                                                                                                                            <  1:4308454/150‑1 (MQ=255)
                              ccATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTggttcggtt                                                                                                                         <  2:3682356/150‑1 (MQ=255)
                                aTTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCt                                                                                                                       >  2:2165929/1‑150 (MQ=255)
                                    aTGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTg                                                                                                                   >  2:470136/1‑150 (MQ=255)
                                    aTGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTg                                                                                                                   <  1:2974127/150‑1 (MQ=255)
                                    aTGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTg                                                                                                                   >  1:3080785/1‑150 (MQ=255)
                                    aTGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTg                                                                                                                   >  1:4049372/1‑150 (MQ=255)
                                        agaACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTc                                                                                                               <  2:3298730/150‑1 (MQ=255)
                                          aaCATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAg                                                                                                             <  1:814090/150‑1 (MQ=255)
                                           aCATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAgg                                                                                                            >  2:4607961/1‑150 (MQ=255)
                                           aCATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAgg                                                                                                            <  1:3310511/150‑1 (MQ=255)
                                            cATTTACAGTGCCGGCACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGa                                                                                                           <  1:4650146/150‑1 (MQ=255)
                                             aTTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGAc                                                                                                          <  1:3258322/150‑1 (MQ=255)
                                               ttACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGt                                                                                                        >  1:4966433/1‑150 (MQ=255)
                                               ttACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGt                                                                                                        >  1:1254621/1‑150 (MQ=255)
                                                tACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGtt                                                                                                       <  2:2640515/150‑1 (MQ=255)
                                                     tGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGa                                                                                                  <  2:1785673/150‑1 (MQ=255)
                                                         gACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCa                                                                                              <  1:2789726/150‑1 (MQ=255)
                                                            acCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGc                                                                                           >  1:2432144/1‑150 (MQ=255)
                                                              cTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCaa                                                                                         <  1:3865727/150‑1 (MQ=255)
                                                              cTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCaa                                                                                         >  2:70601/1‑150 (MQ=255)
                                                               tCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCaaa                                                                                        >  1:2379731/1‑150 (MQ=255)
                                                               tCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCaaa                                                                                        >  2:4017082/1‑150 (MQ=255)
                                                                    ttCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAAtgct                                                                                   <  2:3348530/150‑1 (MQ=255)
                                                                      ccGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAAtgctgc                                                                                 >  2:2503412/1‑150 (MQ=255)
                                                                       cgcgATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTgcg                                                                                >  1:3446000/1‑150 (MQ=255)
                                                                             cGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCt                                                                          >  1:588497/1‑150 (MQ=255)
                                                                              gAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTg                                                                         >  2:2763777/1‑150 (MQ=255)
                                                                              gAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTg                                                                         >  2:4890106/1‑150 (MQ=255)
                                                                              gAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTg                                                                         >  1:1753014/1‑150 (MQ=255)
                                                                               aGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGa                                                                        >  2:3383904/1‑150 (MQ=255)
                                                                               aGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGa                                                                        <  2:309173/150‑1 (MQ=255)
                                                                                      tgtgGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAg                                                                 >  1:3223950/1‑150 (MQ=255)
                                                                                       gtggtgTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAgc                                                                >  2:3904115/1‑150 (MQ=255)
                                                                                               gCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCa                                                        <  1:2624024/150‑1 (MQ=255)
                                                                                                cACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCat                                                       >  1:2696603/1‑150 (MQ=255)
                                                                                                 aCTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCatg                                                      >  1:5172535/1‑150 (MQ=255)
                                                                                                 aCTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCatg                                                      <  1:2674563/150‑1 (MQ=255)
                                                                                                    gTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCatgatg                                                   <  2:1002086/150‑1 (MQ=255)
                                                                                                     tCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGc                                                  <  2:5318295/150‑1 (MQ=255)
                                                                                                      cGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCa                                                 >  2:1049647/1‑150 (MQ=255)
                                                                                                       ggAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACTGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCAc                                                >  1:1914369/1‑150 (MQ=255)
                                                                                                       ggAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCAc                                                >  1:1543161/1‑150 (MQ=255)
                                                                                                        gAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACt                                               <  1:4007765/150‑1 (MQ=255)
                                                                                                            ggcggcAGTTACGACGTGCATTACAGTATTACGCCGTTATGTACTCTGGACGGCTGCATTATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCTGCTGAGCTACAGCGCCTCCCATGATGCACTGACg                                           <  1:2261556/150‑1 (MQ=255)
                                                                                                            ggcggcAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACg                                           <  1:4995130/150‑1 (MQ=255)
                                                                                                               ggcAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCAt                                        >  1:1520831/1‑150 (MQ=255)
                                                                                                                gcAGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCAtc                                       <  2:914067/150‑1 (MQ=255)
                                                                                                                  aGTTACGACGTGCATTACAGTCTTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCAtctc                                     <  2:5287746/150‑1 (MQ=255)
                                                                                                                  aGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCAtctc                                     >  1:2346100/1‑150 (MQ=255)
                                                                                                                  aGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCAtctc                                     <  2:412598/150‑1 (MQ=255)
                                                                                                                  aGTTACGACGTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCAtctc                                     >  2:1561630/1‑150 (MQ=255)
                                                                                                                   gTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCg                                    >  2:4158198/1‑150 (MQ=255)
                                                                                                                   gTTACGACGTGCATTACAGTATTACGCCGTTAAGAACACTGGACGGCAGCAATATTTGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCg                                    >  1:1948880/1‑150 (MQ=255)
                                                                                                                     tACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAATTGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGcc                                  <  1:1450928/150‑1 (MQ=255)
                                                                                                                     tACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGcc                                  <  1:1434227/150‑1 (MQ=255)
                                                                                                                     tACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGcc                                  >  1:3624420/1‑150 (MQ=255)
                                                                                                                       cgacgTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACATCGCCTCCCATGATGCACTGACGCATCTCGCCaa                                >  2:609471/1‑150 (MQ=255)
                                                                                                                         acgTGCATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAAcc                              >  2:2117224/1‑150 (MQ=255)
                                                                                                                             gCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCgcgc                          >  1:922105/1‑150 (MQ=255)
                                                                                                                              cATTACAGTATTACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGcc                         >  2:1330035/1‑150 (MQ=255)
                                                                                                                               aTTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCa                        <  2:43714/150‑1 (MQ=255)
                                                                                                                               aTTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCa                        >  1:3388561/1‑150 (MQ=255)
                                                                                                                                  aCAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACAGACGCATCTCGCCAACCGCGCCAGtt                     <  1:660790/150‑1 (MQ=255)
                                                                                                                                   cAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGttt                    >  2:4445230/1‑150 (MQ=255)
                                                                                                                                   cAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGttt                    <  2:4345605/150‑1 (MQ=255)
                                                                                                                                   cAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGttt                    <  1:2531420/150‑1 (MQ=255)
                                                                                                                                    aGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGtttt                   >  2:4499556/1‑150 (MQ=255)
                                                                                                                                    aGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGtttt                   >  1:1710120/1‑150 (MQ=255)
                                                                                                                                    aGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGtttt                   >  2:2037572/1‑150 (MQ=255)
                                                                                                                                      tATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTga                 >  2:3069465/1‑150 (MQ=255)
                                                                                                                                      tATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTga                 >  2:2999923/1‑150 (MQ=255)
                                                                                                                                       aTTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTgag                >  2:5059262/1‑150 (MQ=255)
                                                                                                                                        ttACGCCGTTAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTgaga               <  2:5136342/150‑1 (MQ=255)
                                                                                                                                          aCGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTAACGAATCACGCAAAATCCTGCGCCAGCTGAGCTACAGAGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGaaa             >  2:5195142/1‑150 (MQ=255)
                                                                                                                                              cGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAaacaac         <  1:2153941/150‑1 (MQ=255)
                                                                                                                                               gTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACt        >  2:2323663/1‑150 (MQ=255)
                                                                                                                                                 tatgTACTCTGGAGGGCAGCATTTTTTGTTCGGTTCTTGTTATTCAGGACGTTACCGAATCTCGAAAAATGCTGCGCCAGCTGAGCTACAGCTCCTCCCATGATGCACTGACTCATCTCGCCAACCGCGCCAGTTTTGTGAAACAACTGc      <  2:3734136/147‑1 (MQ=255)
                                                                                                                                                 tAAGAACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGc      <  1:4669262/150‑1 (MQ=255)
                                                                                                                                                  aaGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCg     >  1:1009089/1‑150 (MQ=255)
                                                                                                                                                  aaGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCg     >  1:3756621/1‑150 (MQ=255)
                                                                                                                                                  aaGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGATTTGAGAAACAACTGCg     >  2:2139614/1‑150 (MQ=255)
                                                                                                                                                    gTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTa   <  1:1317254/150‑1 (MQ=255)
                                                                                                                                                     tACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTAt  <  1:4222674/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTAT  >  CP009273/2138553‑2138851

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: