Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,020,629 C→A 26.2% I50I (ATC→ATA ygfQ → putative purine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,020,6290CA26.2% 563.1 / 192.6 262I50I (ATC→ATAygfQputative purine permease
Reads supporting (aligned to +/- strand):  ref base C (100/93);  new base A (38/31);  total (138/124)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 267 total aligned reads displayed.

GTCTGGAGACATCCTACAAACACCGGACGCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTCCTTGACGGCGTT  >  CP009273/3020482‑3020778
                                                                                                                                                   |                                                                                                                                                     
gTCTGGAGACATCCTACAAACACCGGACGCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATGACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTAtcgt                                                                                                                                                     <  1:4511205/150‑1 (MQ=255)
 tCTGGAGACATCCTACAAACACCGGACGCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTAtagtc                                                                                                                                                    <  1:287239/150‑1 (MQ=255)
          aTCCTACAAACACCGGACGCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGa                                                                                                                                           >  1:3075419/1‑150 (MQ=255)
            ccTACAAACACCGGACGCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAAt                                                                                                                                         <  1:2007649/150‑1 (MQ=255)
               aCAAACACCGGACGCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCt                                                                                                                                      >  1:4210993/1‑150 (MQ=255)
               aCAAACACCGGACGCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTAGTCCGTTATCGTCGTTCCGGGAATGCt                                                                                                                                      <  1:1902696/150‑1 (MQ=255)
                 aaaCACCGGACGCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTgg                                                                                                                                    <  2:2588419/150‑1 (MQ=255)
                       cGGACGCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAg                                                                                                                              <  2:2486498/150‑1 (MQ=255)
                        ggACGCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTGTATTCCGTTATAGTCGTTCCGGGACTGCTGGGCAAAGc                                                                                                                             >  2:4984561/1‑150 (MQ=255)
                         gACGCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCa                                                                                                                            >  2:1935941/1‑150 (MQ=255)
                          aCGCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAg                                                                                                                           <  2:1323838/150‑1 (MQ=255)
                            gCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCgggg                                                                                                                         >  1:4906956/1‑149 (MQ=255)
                            gCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGt                                                                                                                         >  1:2571263/1‑150 (MQ=255)
                            gCACCAAAGCCACAGGGCGCGCTGGATAAATATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGACAGGAAGTACTGGCTGGCGTAACGACCTTTCAGGCCATGGTTTATTCCGTTATAGTCGTTACGGGACTGCTGGGCAAAGCAGGt                                                                                                                         >  2:4123/1‑150 (MQ=255)
                             cACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGtt                                                                                                                        >  1:1533338/1‑150 (MQ=255)
                               ccAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGtttt                                                                                                                      >  2:4412007/1‑150 (MQ=255)
                                cAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTtc                                                                                                                     <  2:1237564/150‑1 (MQ=255)
                                 aaaGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTtcc                                                                                                                    >  1:779447/1‑150 (MQ=255)
                                  aaGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTtcct                                                                                                                   <  2:4322085/150‑1 (MQ=255)
                                   aGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTtcctc                                                                                                                  >  1:1122764/1‑150 (MQ=255)
                                      cacaGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTACGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCg                                                                                                               >  2:2887556/1‑150 (MQ=255)
                                         aGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAg                                                                                                            >  1:1725104/1‑150 (MQ=255)
                                         aGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAg                                                                                                            <  2:4680724/150‑1 (MQ=255)
                                              gcgcTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCtgtg                                                                                                       <  2:406739/150‑1 (MQ=255)
                                              gcgcTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCtgtg                                                                                                       <  2:2099888/150‑1 (MQ=255)
                                              gcgcTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCtgtg                                                                                                       <  2:2684099/150‑1 (MQ=255)
                                                gcTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGtgtt                                                                                                     <  1:1031782/150‑1 (MQ=255)
                                                 cTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGtgttt                                                                                                    <  2:978509/150‑1 (MQ=255)
                                                 cTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACAGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGtgttt                                                                                                    >  2:3310363/1‑150 (MQ=255)
                                                   ggATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGtgtttgt                                                                                                  >  1:2478895/1‑150 (MQ=255)
                                                    gATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGtgtttgtt                                                                                                 >  1:2527273/1‑150 (MQ=255)
                                                     aTAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTttgttg                                                                                                >  1:731605/1‑150 (MQ=255)
                                                     aTAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTttgttg                                                                                                <  2:589473/150‑1 (MQ=255)
                                                      tAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGc                                                                                               >  1:1409416/1‑150 (MQ=255)
                                                        attattTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCa                                                                                             >  1:67092/1‑150 (MQ=255)
                                                         ttattTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCAc                                                                                            <  2:879739/150‑1 (MQ=255)
                                                           attTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCt                                                                                          >  1:4061567/1‑150 (MQ=255)
                                                              ttAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTc                                                                                       >  2:1574554/1‑150 (MQ=255)
                                                               tAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCt                                                                                      <  1:2830895/150‑1 (MQ=255)
                                                               tAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCt                                                                                      <  2:4016392/150‑1 (MQ=255)
                                                                aaaaTTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTg                                                                                     <  1:1093694/150‑1 (MQ=255)
                                                                 aaaTTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTgg                                                                                    >  2:669610/1‑150 (MQ=255)
                                                                  aaTTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGt                                                                                   <  2:1027597/150‑1 (MQ=255)
                                                                    ttACCGCTCGTGGCAGTCCCGTTCGTCAGGAAGGACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTAGTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTcg                                                                                 <  1:2436481/150‑1 (MQ=255)
                                                                    ttACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTGTCTGGCCATGGTTTATTCCTTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTcg                                                                                 >  2:1488257/1‑150 (MQ=255)
                                                                     tACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTcgc                                                                                >  1:2218607/1‑150 (MQ=255)
                                                                         gCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGc                                                                            >  1:1924897/1‑150 (MQ=255)
                                                                            cGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTc                                                                         <  1:2965065/150‑1 (MQ=255)
                                                                            cGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTc                                                                         <  1:1935941/150‑1 (MQ=255)
                                                                            cGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCAGGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCGCGCAGCTGTGTTTGTTGCCACCTGGCTGGTCGCGGGCTTc                                                                         >  2:2806368/1‑150 (MQ=255)
                                                                               ggCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGc                                                                      >  2:1775992/1‑150 (MQ=255)
                                                                               ggCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGc                                                                      >  2:1116236/1‑150 (MQ=255)
                                                                                  aGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCg                                                                   >  1:3899032/1‑150 (MQ=255)
                                                                                  aGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCg                                                                   >  1:1950524/1‑150 (MQ=255)
                                                                                   gTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGt                                                                  >  2:2998137/1‑150 (MQ=255)
                                                                                   gTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCACCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGt                                                                  >  1:1447942/1‑150 (MQ=255)
                                                                                    tACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGtt                                                                 <  1:1308516/150‑1 (MQ=255)
                                                                                     aCCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCGGTCTGGTCGCGGGCTTCGGCTCGTTg                                                                >  2:3155241/1‑150 (MQ=255)
                                                                                       cGTTCGTCAGGAAGTTCTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCt                                                              <  1:913964/150‑1 (MQ=255)
                                                                                        gTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTg                                                             >  2:2425764/1‑150 (MQ=255)
                                                                                         ttCGTCAGGAAGTACTGGCTGTCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGa                                                            <  1:4483620/150‑1 (MQ=255)
                                                                                           cGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATg                                                          >  1:2069716/1‑150 (MQ=255)
                                                                                           cGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATg                                                          >  1:4674353/1‑150 (MQ=255)
                                                                                            gTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGACTGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGGTGCTGATgg                                                         >  2:3123043/1‑150 (MQ=255)
                                                                                            gTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATgg                                                         >  2:3123127/1‑150 (MQ=255)
                                                                                                  caGTCCTGGCTGGCCTCACGCCCTCTCTGGCCCGGGTTTCTTCCGTGATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTAt                                                   <  1:2689763/149‑1 (MQ=255)
                                                                                                  aaGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATAATCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTAt                                                   <  2:1764403/150‑1 (MQ=255)
                                                                                                  aaGTACTGGCTGGCTTAACGACCTTTCTGGCCAGGGTTTATTCCGTTATCGTCGGTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTAt                                                   <  2:5011025/150‑1 (MQ=255)
                                                                                                     tACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATggg                                                >  1:5192487/1‑150 (MQ=255)
                                                                                                     tACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATggg                                                >  1:516598/1‑150 (MQ=255)
                                                                                                     tACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTAGTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATggg                                                <  1:4281281/150‑1 (MQ=255)
                                                                                                      aCTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGc                                               <  2:1835050/150‑1 (MQ=255)
                                                                                                        tggctggcTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTa                                             >  1:3873656/1‑150 (MQ=255)
                                                                                                         ggctggctTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTaa                                            <  2:2061037/150‑1 (MQ=255)
                                                                                                          gctggctTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAAt                                           >  2:3475651/1‑150 (MQ=255)
                                                                                                            tggctTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAAttt                                         >  2:1368595/1‑150 (MQ=255)
                                                                                                            tggctTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAAttt                                         >  1:5065698/1‑150 (MQ=255)
                                                                                                            tggctTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAAttt                                         <  2:3847288/150‑1 (MQ=255)
                                                                                                                 tAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCaa                                    <  2:3600203/150‑1 (MQ=255)
                                                                                                                  aaCGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAAt                                   >  2:2588826/1‑150 (MQ=255)
                                                                                                                     gACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGc                                >  1:3605697/1‑150 (MQ=255)
                                                                                                                     gACCTTTCGGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGc                                <  2:5074866/150‑1 (MQ=255)
                                                                                                                        cTTTCTGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGAt                             >  1:4096032/1‑150 (MQ=255)
                                                                                                                             tGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGtt                        >  1:587777/1‑150 (MQ=255)
                                                                                                                             tGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGtt                        >  2:5191746/1‑150 (MQ=255)
                                                                                                                             tGGCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAAAGGCGATTGGtt                        >  2:3269409/1‑150 (MQ=255)
                                                                                                                              ggCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTg                       >  1:4880991/1‑150 (MQ=255)
                                                                                                                              ggCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCGTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTg                       >  2:4833252/1‑150 (MQ=255)
                                                                                                                              ggCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTg                       <  2:2231968/150‑1 (MQ=255)
                                                                                                                               gCCATGGTTTATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTgc                      <  2:738683/150‑1 (MQ=255)
                                                                                                                                ccATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTgcg                     <  1:3167499/150‑1 (MQ=255)
                                                                                                                                  aTGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGcgcg                   >  1:3758736/1‑150 (MQ=255)
                                                                                                                                      tttATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGAttt               >  2:4647827/1‑150 (MQ=255)
                                                                                                                                       ttATTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTc              >  2:4664010/1‑150 (MQ=255)
                                                                                                                                         aTTCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTCCt            >  2:1604992/1‑150 (MQ=255)
                                                                                                                                          ttCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTCCtt           <  1:4437847/150‑1 (MQ=255)
                                                                                                                                           tCCGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTCCTTg          <  1:4654203/150‑1 (MQ=255)
                                                                                                                                             cGTTATAGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTCCTTGAc        >  1:4731458/1‑150 (MQ=255)
                                                                                                                                                tATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGGTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTCCTTGACGGc     <  2:927481/150‑1 (MQ=255)
                                                                                                                                                  tcgtcgTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTGGCGCGATTTCATTGACGGCGt   >  2:4817569/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTCTGGAGACATCCTACAAACACCGGACGCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTCCTTGACGGCGTT  >  CP009273/3020482‑3020778

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: