Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,186,181 G→C 79.9% W205C (TGG→TGC)  yqiJ → DUF1449 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,186,1810GC79.9% 590.3 / 153.3 214W205C (TGG→TGC) yqiJDUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (18/25);  new base C (95/76);  total (113/101)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 214 total aligned reads displayed.

GCGAAGGAAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCG  >  CP009273/3186036‑3186328
                                                                                                                                                 |                                                                                                                                                   
gCGAAGGAAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCcccc                                                                                                                                                 >  2:164053/1‑150 (MQ=255)
 cgaacgaaCACTCACCGATCACTATGTCCAAATACATTATGCACGACGGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCa                                                                                                                                                <  2:2949021/150‑1 (MQ=255)
 cGAAGGAAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCa                                                                                                                                                >  1:4797464/1‑150 (MQ=255)
 cGAAGGAAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCa                                                                                                                                                <  1:298804/150‑1 (MQ=255)
  gAAGGAAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGGCCCCaa                                                                                                                                               >  1:2206209/1‑150 (MQ=255)
  gAAGGAAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCaa                                                                                                                                               >  1:342821/1‑150 (MQ=255)
  gAAGGAAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCACTGCCCCCaa                                                                                                                                               >  1:487089/1‑150 (MQ=255)
       aaaaCTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATAtt                                                                                                                                          >  1:2995114/1‑150 (MQ=255)
          aCTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGGGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTa                                                                                                                                       >  1:91995/1‑150 (MQ=255)
          aCTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTa                                                                                                                                       >  2:4906083/1‑150 (MQ=255)
           cTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTaa                                                                                                                                      <  1:1131941/150‑1 (MQ=255)
                  aTCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATttatta                                                                                                                               >  2:2545063/1‑150 (MQ=255)
                    cAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAAt                                                                                                                             >  1:4026129/1‑150 (MQ=255)
                     aaTTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATc                                                                                                                            >  2:2472498/1‑150 (MQ=255)
                      aTTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCa                                                                                                                           >  1:3323482/1‑150 (MQ=255)
                         tGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCaaaa                                                                                                                        >  1:964241/1‑150 (MQ=255)
                          ggCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAg                                                                                                                       <  1:1179983/150‑1 (MQ=255)
                            ccAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGa                                                                                                                     <  2:2787910/150‑1 (MQ=255)
                                   cATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATGTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTtaa                                                                                                              <  2:1870143/150‑1 (MQ=255)
                                      tattTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTtaataa                                                                                                           >  1:1799081/1‑150 (MQ=255)
                                       attTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTaataat                                                                                                          >  2:543669/1‑150 (MQ=255)
                                       attTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTaataat                                                                                                          >  2:1646740/1‑150 (MQ=255)
                                        tttACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTTCCGACTATCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATAGTAATTTATTAATCAAAAGGAAATATAATAATg                                                                                                         >  2:3620916/1‑150 (MQ=255)
                                        tttACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATg                                                                                                         >  2:2393724/1‑150 (MQ=255)
                                        tttACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATg                                                                                                         >  2:4666449/1‑150 (MQ=255)
                                          tactacTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGga                                                                                                       <  1:3896659/150‑1 (MQ=255)
                                               ctGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGata                                                                                                  >  2:5061770/1‑150 (MQ=255)
                                               ctGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGata                                                                                                  >  1:2262694/1‑150 (MQ=255)
                                                 ggAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATAtt                                                                                                <  1:183953/150‑1 (MQ=255)
                                                  gAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTg                                                                                               <  2:1469875/150‑1 (MQ=255)
                                                    gCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGtt                                                                                             <  2:3610436/150‑1 (MQ=255)
                                                     ccTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTa                                                                                            >  1:470555/1‑150 (MQ=255)
                                                      cTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTaa                                                                                           >  1:3622508/1‑150 (MQ=255)
                                                      cTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTaa                                                                                           >  2:2505480/1‑150 (MQ=255)
                                                       tGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAAt                                                                                          >  1:4044405/1‑150 (MQ=255)
                                                         aagaagGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTc                                                                                        >  2:2127785/1‑150 (MQ=255)
                                                          agaagGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCt                                                                                       >  2:2226146/1‑150 (MQ=255)
                                                          agaagGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCt                                                                                       >  1:252061/1‑150 (MQ=255)
                                                            aagGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGAGATTGTTAATTCtgt                                                                                     >  2:3940092/1‑150 (MQ=255)
                                                            aagGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCtgt                                                                                     <  2:3155669/150‑1 (MQ=255)
                                                            aagGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCtgt                                                                                     <  2:1064798/150‑1 (MQ=255)
                                                             agGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCtgtg                                                                                    >  2:2521636/1‑150 (MQ=255)
                                                                 aaaaTTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGccca                                                                                <  1:4288928/150‑2 (MQ=255)
                                                                 aaaaTTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCcct                                                                                <  2:320085/150‑1 (MQ=255)
                                                                     ttttCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTg                                                                            <  1:2381467/150‑1 (MQ=255)
                                                                         cACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATg                                                                        >  2:3358955/1‑150 (MQ=255)
                                                                           cTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATAATCTAATTTATTAATCAAAATGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGtt                                                                      >  2:2826251/1‑150 (MQ=255)
                                                                            tAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGttt                                                                     >  2:3540164/1‑150 (MQ=255)
                                                                             aaaGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTa                                                                    >  1:4980395/1‑150 (MQ=255)
                                                                               aGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTAcc                                                                  >  1:4444505/1‑150 (MQ=255)
                                                                                ggAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCg                                                                 <  1:1015074/150‑1 (MQ=255)
                                                                                ggAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGTATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCg                                                                 <  2:1794214/150‑1 (MQ=255)
                                                                                  agaCAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCa                                                               >  1:4122664/1‑150 (MQ=255)
                                                                                      aaaGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAatta                                                           >  1:3167190/1‑150 (MQ=255)
                                                                                       aaGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAattat                                                          <  2:3538899/150‑1 (MQ=255)
                                                                                        aGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAattatt                                                         <  2:47484/150‑1 (MQ=255)
                                                                                        aGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAattatt                                                         <  1:4908668/150‑1 (MQ=255)
                                                                                         gtgtTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTg                                                        <  2:40998/150‑1 (MQ=255)
                                                                                           gtTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGcc                                                      <  2:2383112/150‑1 (MQ=255)
                                                                                            ttAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCg                                                     >  2:4165972/1‑150 (MQ=255)
                                                                                            ttAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCg                                                     <  2:2649983/150‑1 (MQ=255)
                                                                                              aaTTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTa                                                   <  2:2172452/150‑1 (MQ=255)
                                                                                              aaTTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTa                                                   <  2:2291218/150‑1 (MQ=255)
                                                                                              aaTTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTa                                                   >  1:3949535/1‑150 (MQ=255)
                                                                                              aaTTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTa                                                   <  2:4393790/150‑1 (MQ=255)
                                                                                              aaTTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTa                                                   >  1:1897776/1‑150 (MQ=255)
                                                                                                     tGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTc                                            >  1:910272/1‑150 (MQ=255)
                                                                                                        cGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGt                                         >  2:3731034/1‑150 (MQ=255)
                                                                                                        cGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGt                                         >  2:3087649/1‑150 (MQ=255)
                                                                                                         gACTTTCTGCAACGCGATATTTGTCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCGGGCTGTTTACCGCCATTCTTGCCGTATGCATTCTGtt                                        >  2:3397994/1‑150 (MQ=255)
                                                                                                         gACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGtt                                        <  2:4132330/150‑1 (MQ=255)
                                                                                                             ttCTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGTTGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTttatt                                    <  1:3758261/150‑1 (MQ=255)
                                                                                                              tCTGCAACGCGATATTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTttatta                                   <  1:3325795/150‑1 (MQ=255)
                                                                                                               cTGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTtattat                                  >  1:1378168/1‑150 (MQ=255)
                                                                                                                tGCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTattatt                                 >  1:1473845/1‑150 (MQ=255)
                                                                                                                 gCAACGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATCATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTg                                >  2:4375204/1‑150 (MQ=255)
                                                                                                                   aaCGCGATATTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGa                              >  2:4234178/1‑150 (MQ=255)
                                                                                                                     cgcgATATTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAat                            >  1:1873718/1‑150 (MQ=255)
                                                                                                                         atatTTGGCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAattatt                        >  2:3258235/1‑150 (MQ=255)
                                                                                                                          tatTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTAttt                       <  2:3100348/150‑1 (MQ=255)
                                                                                                                           atTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTAtttt                      <  1:1474178/150‑1 (MQ=255)
                                                                                                                                gCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCa                 <  1:2127785/150‑1 (MQ=255)
                                                                                                                                gCGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCAGTCTGTTTATTATTGGAATTATTTTCGCCa                 >  1:4752648/1‑150 (MQ=255)
                                                                                                                                 cGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAg                <  1:2541027/150‑1 (MQ=255)
                                                                                                                                 cGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAg                <  2:1052271/150‑1 (MQ=255)
                                                                                                                                 cGGAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAg                <  1:1513781/150‑1 (MQ=255)
                                                                                                                                   gAAAATAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAACGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCCTTCTGTTTATTATTGGACTTATTTTCGCCAGGc              >  2:1697061/1‑150 (MQ=255)
                                                                                                                                    aaaaTAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGct             <  1:4678233/150‑1 (MQ=255)
                                                                                                                                    aaaaTAATCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGct             <  2:342821/150‑1 (MQ=255)
                                                                                                                                        taatCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTAt         <  2:1791075/150‑1 (MQ=255)
                                                                                                                                        taatCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATCCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTAt         <  1:2672207/150‑1 (MQ=255)
                                                                                                                                        taatCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTAt         >  2:4568205/1‑150 (MQ=255)
                                                                                                                                         aatCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTAtc        <  1:164066/150‑1 (MQ=255)
                                                                                                                                         aatCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTAtc        <  1:2621823/150‑1 (MQ=255)
                                                                                                                                          atCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCAGTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTAtcg       >  2:4404974/1‑150 (MQ=255)
                                                                                                                                          atCCCTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTAtcg       <  2:1809535/150‑1 (MQ=255)
                                                                                                                                            cccTGGCCCCAAATCTTGTAATTTATTAATCAAAAGGAAATTTAATAATGGAGGATATTGTTAATTCTGTGCCCTCCTGGATGTGTACCGCAATTATTGCCGGATGCATTCTGTTTAATATTGGAATTATGTTCGCCAGGCTCTAtcgtc     >  1:815157/1‑150 (MQ=255)
                                                                                                                                             ccTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTAtcgtcg    <  1:3801812/150‑1 (MQ=255)
                                                                                                                                              cTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTcgc   >  2:4512718/1‑150 (MQ=255)
                                                                                                                                              cTGCCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTcgc   >  2:2166703/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCGAAGGAAAACTCACCGATCAATTTGGCCAAATACATTATTTACTACTGGAGCCTGAAGAAGGAAAAATTTTCACTAAAGGAGACAAAGTGTTAATTATTTGCCGACTTTCTGCAACGCGATATTTGGCGGAAAATAATCCCTGGCCCCAAATATTGTAATTTATTAATCAAAAGGAAATTTAATAATGGATGATATTGTTAATTCTGTGCCCTCCTGGATGTTTACCGCAATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCG  >  CP009273/3186036‑3186328

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: ‑