Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,936,686 G→T 5.2% intergenic (‑182/‑10) ygdI ← / → csdA DUF903 family verified lipoprotein/cysteine sulfinate desulfinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,936,6860GT 5.2% 724.6 / 11.8 190intergenic (‑182/‑10)ygdI/csdADUF903 family verified lipoprotein/cysteine sulfinate desulfinase
Reads supporting (aligned to +/- strand):  ref base G (85/95);  new base T (7/3);  total (92/98)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 199 total aligned reads displayed.

AGTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAG  >  CP009273/2936539‑2936832
                                                                                                                                                   |                                                                                                                                                  
aGTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCg                                                                                                                                                    <  2:2344579/148‑1 (MQ=255)
 gTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCga                                                                                                                                                   <  2:2672721/148‑1 (MQ=255)
 gTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCga                                                                                                                                                   <  2:3178318/148‑1 (MQ=255)
 gTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCga                                                                                                                                                   <  2:2430318/148‑1 (MQ=255)
 gTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCga                                                                                                                                                   <  1:1948944/148‑1 (MQ=255)
 gTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCga                                                                                                                                                   >  1:384016/1‑148 (MQ=255)
  tATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCgag                                                                                                                                                  >  1:2908984/1‑148 (MQ=255)
      gATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGt                                                                                                                                              <  1:2953657/148‑1 (MQ=255)
       aTTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTa                                                                                                                                             <  1:355519/148‑1 (MQ=255)
       aTTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTa                                                                                                                                             >  1:3405828/1‑148 (MQ=255)
       aTTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTa                                                                                                                                             >  2:1037532/1‑148 (MQ=255)
       aTTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTa                                                                                                                                             <  1:355560/148‑1 (MQ=255)
          aaCAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCa                                                                                                                                          <  2:652449/148‑1 (MQ=255)
             aGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGa                                                                                                                                       <  2:2772750/148‑1 (MQ=255)
             aGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGa                                                                                                                                       <  1:3242007/148‑1 (MQ=255)
                cTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACg                                                                                                                                    <  1:1401260/148‑1 (MQ=255)
                  aaaaTTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGtt                                                                                                                                  <  1:3020853/148‑1 (MQ=255)
                  aaaaTTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGtt                                                                                                                                  >  2:349682/1‑148 (MQ=255)
                  aaaaTTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGtt                                                                                                                                  <  2:3104080/148‑1 (MQ=255)
                   aaaTTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGttt                                                                                                                                 <  2:3276095/148‑1 (MQ=255)
                         aCGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATc                                                                                                                           <  2:2301091/148‑1 (MQ=255)
                           gCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATccc                                                                                                                         <  1:1239418/148‑1 (MQ=255)
                            ccTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCcg                                                                                                                        >  2:1486035/1‑148 (MQ=255)
                              tAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCcgcg                                                                                                                      >  1:2024978/1‑148 (MQ=255)
                              tAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCcgcg                                                                                                                      >  1:2024477/1‑148 (MQ=255)
                               aaCACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCgcgc                                                                                                                     >  2:3120945/1‑148 (MQ=255)
                                acacTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTGTTTAATCCCGCGCa                                                                                                                    >  1:905688/1‑148 (MQ=255)
                                 cacTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAg                                                                                                                   >  1:3428145/1‑148 (MQ=255)
                                  acTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGt                                                                                                                  <  1:3273212/148‑1 (MQ=255)
                                  acTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGt                                                                                                                  <  1:2286076/148‑1 (MQ=255)
                                  acTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGt                                                                                                                  <  1:2673629/148‑1 (MQ=255)
                                  acTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGt                                                                                                                  <  1:3357624/148‑1 (MQ=255)
                                    tATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGttt                                                                                                                >  1:2729102/1‑148 (MQ=255)
                                      ttCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCTAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTcg                                                                                                              >  2:509284/1‑148 (MQ=255)
                                            atatGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCa                                                                                                        >  1:666376/1‑148 (MQ=255)
                                            atatGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCa                                                                                                        >  2:747409/1‑148 (MQ=255)
                                               tGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGtt                                                                                                     >  1:2194366/1‑148 (MQ=255)
                                                gTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGttt                                                                                                    <  1:782818/148‑1 (MQ=255)
                                                 ttacttacGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTc                                                                                                   >  2:470562/1‑148 (MQ=255)
                                                 ttacttacGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTc                                                                                                   <  2:2785328/148‑1 (MQ=255)
                                                     ttacGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGc                                                                                               >  2:2818997/1‑148 (MQ=255)
                                                     ttacGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGc                                                                                               >  2:2818525/1‑148 (MQ=255)
                                                      tacGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCa                                                                                              >  2:1395393/1‑148 (MQ=255)
                                                       acGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCTAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCAc                                                                                             >  1:2241222/1‑148 (MQ=255)
                                                        cgcgGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACt                                                                                            >  2:1259650/1‑148 (MQ=255)
                                                          cgGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTAc                                                                                          <  1:2869783/148‑1 (MQ=255)
                                                                cGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATg                                                                                    >  2:938568/1‑148 (MQ=255)
                                                                 gCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGc                                                                                   <  1:187384/148‑1 (MQ=255)
                                                                  ccAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCg                                                                                  >  1:533762/1‑148 (MQ=255)
                                                                    aGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCggg                                                                                >  2:3266706/1‑148 (MQ=255)
                                                                       atatCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGt                                                                             >  2:2699141/1‑148 (MQ=255)
                                                                        tatCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTc                                                                            <  1:3013300/148‑1 (MQ=255)
                                                                         atCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCt                                                                           <  1:1770322/148‑1 (MQ=255)
                                                                          tCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTa                                                                          >  2:1931815/1‑148 (MQ=255)
                                                                           ccAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTAt                                                                         >  1:1482899/1‑148 (MQ=255)
                                                                           ccAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTAt                                                                         <  1:866124/148‑1 (MQ=255)
                                                                           ccAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTAt                                                                         <  1:1016095/148‑1 (MQ=255)
                                                                           ccAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTAt                                                                         >  2:220618/1‑148 (MQ=255)
                                                                           ccAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTAt                                                                         <  2:2988532/148‑1 (MQ=255)
                                                                           ccAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAACTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTAt                                                                         >  2:1522998/1‑148 (MQ=255)
                                                                               aTAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCg                                                                     >  1:1592470/1‑148 (MQ=255)
                                                                                tAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGa                                                                    >  2:1320552/1‑148 (MQ=255)
                                                                                     cTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCTAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAgcg                                                               <  1:122342/148‑1 (MQ=255)
                                                                                      tGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAgcgc                                                              >  2:105827/1‑148 (MQ=255)
                                                                                       gATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGcgcc                                                             >  2:548883/1‑148 (MQ=255)
                                                                                            ttAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGAc                                                        >  1:3348435/1‑148 (MQ=255)
                                                                                             tAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGAcc                                                       <  2:666619/148‑1 (MQ=255)
                                                                                               gCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACcgc                                                     >  1:1121164/1‑148 (MQ=255)
                                                                                               gCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACcgc                                                     >  1:976125/1‑148 (MQ=255)
                                                                                                 aCTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCgcgc                                                   <  1:2567794/148‑1 (MQ=255)
                                                                                                 aCTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCgcgc                                                   >  1:897423/1‑148 (MQ=255)
                                                                                                  ctctCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCt                                                  >  2:1113493/1‑148 (MQ=255)
                                                                                                   tctcCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCtt                                                 >  2:1187237/1‑148 (MQ=255)
                                                                                                    ctcCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCTAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTa                                                >  1:209377/1‑148 (MQ=255)
                                                                                                        ttGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAAcc                                            <  1:1947106/148‑1 (MQ=255)
                                                                                                          gTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTg                                          <  2:3348979/148‑1 (MQ=255)
                                                                                                            aTCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGaa                                        >  1:1051233/1‑148 (MQ=255)
                                                                                                            aTCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGaa                                        >  1:3122511/1‑148 (MQ=255)
                                                                                                             tCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAg                                       <  2:2024978/148‑1 (MQ=255)
                                                                                                             tCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAg                                       <  1:1929719/148‑1 (MQ=255)
                                                                                                                aaGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCg                                    <  2:3088158/148‑1 (MQ=255)
                                                                                                                 aGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGt                                   <  1:224144/148‑1 (MQ=255)
                                                                                                                   tGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTgg                                 <  1:1548862/148‑1 (MQ=255)
                                                                                                                   tGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTgg                                 >  2:36637/1‑148 (MQ=255)
                                                                                                                       ttttGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGa                             >  2:2561694/1‑148 (MQ=255)
                                                                                                                         ttGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAg                           >  2:2247314/1‑148 (MQ=255)
                                                                                                                           gCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGcc                         <  1:349682/148‑1 (MQ=255)
                                                                                                                               cACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCAccc                     >  1:1729755/1‑148 (MQ=255)
                                                                                                                                aCGATCGGTGCATCAAGCCTAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCa                    <  1:509284/148‑1 (MQ=255)
                                                                                                                                    tCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAg                >  1:296125/1‑148 (MQ=255)
                                                                                                                                    tCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAg                <  1:2120813/148‑1 (MQ=255)
                                                                                                                                      ggTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGtt              <  1:626634/148‑1 (MQ=255)
                                                                                                                                       gTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGttt             >  1:3037789/1‑148 (MQ=255)
                                                                                                                                        tGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGtttt            >  1:2411079/1‑148 (MQ=255)
                                                                                                                                        tGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGtttt            <  1:2392335/148‑1 (MQ=255)
                                                                                                                                             cAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGt       <  1:3338031/148‑1 (MQ=255)
                                                                                                                                               aGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCt     >  2:2832872/1‑148 (MQ=255)
                                                                                                                                                gCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGTGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTg    <  1:3099467/148‑1 (MQ=255)
                                                                                                                                                  cGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAg  <  2:3428145/148‑1 (MQ=255)
                                                                                                                                                  cGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAg  <  2:1663157/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAG  >  CP009273/2936539‑2936832

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: