Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 602,195 C→A 7.4% T215T (ACG→ACT ybdK ← weak gamma‑glutamyl:cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273602,1950CA 7.4% 895.4 / 30.6 237T215T (ACG→ACTybdKweak gamma‑glutamyl:cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (124/95);  new base A (7/11);  total (131/106)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

GTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGG  >  CP009273/602047‑602343
                                                                                                                                                    |                                                                                                                                                    
gTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGAACATCg                                                                                                                                                     <  2:2224605/150‑1 (MQ=255)
 tAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCgt                                                                                                                                                    >  1:1797752/1‑150 (MQ=255)
  aaCCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCgtg                                                                                                                                                   >  2:3731920/1‑150 (MQ=255)
    ccAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCgtggt                                                                                                                                                 <  2:2947221/150‑1 (MQ=255)
     cAGTGGGTGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCgtggtg                                                                                                                                                >  2:4216007/1‑150 (MQ=255)
     cAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCgtggtg                                                                                                                                                <  2:3444954/150‑1 (MQ=255)
      aGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGgtgt                                                                                                                                               <  2:1024447/150‑1 (MQ=255)
         gggCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGTTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAAc                                                                                                                                            <  2:1748920/150‑1 (MQ=255)
          ggCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACt                                                                                                                                           <  2:1085607/150‑1 (MQ=255)
            cGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCa                                                                                                                                         <  2:3127013/150‑1 (MQ=255)
            cGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCa                                                                                                                                         <  2:3127006/150‑1 (MQ=255)
             ggTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATAGTGGTGTAACTCAg                                                                                                                                        >  1:2863113/1‑150 (MQ=255)
                aGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACa                                                                                                                                     >  1:4302625/1‑150 (MQ=255)
                 gCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAg                                                                                                                                    <  1:209011/150‑1 (MQ=255)
                     gAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGaa                                                                                                                                >  1:2696842/1‑150 (MQ=255)
                       aTTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAAc                                                                                                                              >  2:3360209/1‑150 (MQ=255)
                        ttAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACa                                                                                                                             >  2:4559612/1‑150 (MQ=255)
                         tAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAg                                                                                                                            <  1:3485355/150‑1 (MQ=255)
                          aaTCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAgg                                                                                                                           >  2:3633019/1‑150 (MQ=255)
                          aaTCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAAAAgg                                                                                                                           >  2:2670278/1‑150 (MQ=255)
                           aTCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAggg                                                                                                                          <  2:1085152/150‑1 (MQ=255)
                            tCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGc                                                                                                                         >  1:302160/1‑150 (MQ=255)
                            tCCCGCCATATTTACTGCGTGGATAAGGGTTAACGGGGTATCCATCACCCGTTCCTCCACCGTGGCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGc                                                                                                                         <  2:3333192/150‑1 (MQ=255)
                              ccGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCtt                                                                                                                       <  1:1147948/150‑1 (MQ=255)
                                gCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCa                                                                                                                     >  2:1329607/1‑150 (MQ=255)
                                 ccATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCaa                                                                                                                    >  2:783999/1‑150 (MQ=255)
                                  cATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCaaa                                                                                                                   >  2:394684/1‑150 (MQ=255)
                                  cATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCaaa                                                                                                                   <  1:1516621/150‑1 (MQ=255)
                                    tatTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAAtt                                                                                                                 <  2:1236454/150‑1 (MQ=255)
                                    tatTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAAtt                                                                                                                 >  1:1007127/1‑150 (MQ=255)
                                    tatTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAAtt                                                                                                                 <  1:392524/150‑1 (MQ=255)
                                     atTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAAttg                                                                                                                >  1:1550612/1‑150 (MQ=255)
                                      tttACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAAttgt                                                                                                               <  2:2572948/150‑1 (MQ=255)
                                       ttACTGCGTGGCTAAGGGTTAACGGGGTATNCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAAttgtt                                                                                                              >  1:4275163/1‑150 (MQ=255)
                                         aCTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGc                                                                                                            <  2:2345970/150‑1 (MQ=255)
                                                    aaGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGAc                                                                                                 >  1:4060950/1‑150 (MQ=255)
                                                    aaGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATAGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGAc                                                                                                 >  1:3119911/1‑150 (MQ=255)
                                                       ggTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCa                                                                                              >  1:3294777/1‑150 (MQ=255)
                                                       ggTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCa                                                                                              >  1:1874286/1‑150 (MQ=255)
                                                         ttAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACg                                                                                            <  2:2426543/150‑1 (MQ=255)
                                                         ttAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACg                                                                                            >  2:443885/1‑150 (MQ=255)
                                                              ggggTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATc                                                                                       <  1:4428914/150‑1 (MQ=255)
                                                              ggggTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATc                                                                                       <  1:3110236/150‑1 (MQ=255)
                                                              ggggTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATc                                                                                       <  1:3070210/150‑1 (MQ=255)
                                                                ggTATCCATCACCCGACCCTCCACCGTGCCACAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCgg                                                                                     <  1:821771/150‑1 (MQ=255)
                                                                ggTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCgg                                                                                     >  1:2503554/1‑150 (MQ=255)
                                                                 gTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCggg                                                                                    >  1:2608352/1‑150 (MQ=255)
                                                                     ccATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCAt                                                                                <  1:2764086/150‑1 (MQ=255)
                                                                       aTCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTa                                                                              >  1:1809161/1‑150 (MQ=255)
                                                                       aTCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTa                                                                              >  2:896771/1‑150 (MQ=255)
                                                                        tCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTAt                                                                             >  2:3419860/1‑150 (MQ=255)
                                                                        tCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTAt                                                                             <  1:409615/150‑1 (MQ=255)
                                                                         cACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATc                                                                            >  1:2312089/1‑150 (MQ=255)
                                                                            ccGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAgg                                                                         <  1:925052/150‑1 (MQ=255)
                                                                             cGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGa                                                                        >  1:3542374/1‑150 (MQ=255)
                                                                              gAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGaa                                                                       >  1:464963/1‑150 (MQ=255)
                                                                              gAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGaa                                                                       >  2:2644833/1‑150 (MQ=255)
                                                                               aaCCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGaaa                                                                      >  2:1288177/1‑150 (MQ=255)
                                                                               aaCCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGaaa                                                                      >  1:2563830/1‑150 (MQ=255)
                                                                                aCCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAg                                                                     >  1:3642779/1‑150 (MQ=255)
                                                                                aCCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAg                                                                     >  2:1534101/1‑150 (MQ=255)
                                                                                 ccTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAgg                                                                    >  2:4126389/1‑150 (MQ=255)
                                                                                  cTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGc                                                                   <  1:2657364/150‑1 (MQ=255)
                                                                                   tCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCg                                                                  >  1:3806942/1‑150 (MQ=255)
                                                                                     cACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGa                                                                <  2:2648252/150‑1 (MQ=255)
                                                                                     cACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGa                                                                >  2:3889317/1‑150 (MQ=255)
                                                                                     cACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGa                                                                <  1:2410732/150‑1 (MQ=255)
                                                                                          tGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGaaaaaa                                                           >  2:219578/1‑150 (MQ=255)
                                                                                              aaaaTGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATAtt                                                       <  2:3551666/150‑1 (MQ=255)
                                                                                                aaTGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCg                                                     >  1:539645/1‑150 (MQ=255)
                                                                                                aaTGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATAGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCg                                                     >  2:1857347/1‑150 (MQ=255)
                                                                                                  tGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATTCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGt                                                   <  2:701686/150‑1 (MQ=255)
                                                                                                   gAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTc                                                  <  2:3402183/150‑1 (MQ=255)
                                                                                                       aCTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTga                                              >  1:744973/1‑150 (MQ=255)
                                                                                                       aCTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTga                                              <  2:1961301/150‑1 (MQ=255)
                                                                                                       aCTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAAAATTCGGTCGTga                                              >  1:3234702/1‑150 (MQ=255)
                                                                                                        cTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTgag                                             <  2:3060270/150‑1 (MQ=255)
                                                                                                           gggCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTgaggag                                          <  2:366778/150‑1 (MQ=255)
                                                                                                            ggCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGaggagg                                         <  1:4147791/150‑1 (MQ=255)
                                                                                                             gCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGc                                        <  1:981926/150‑1 (MQ=255)
                                                                                                               gAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCaa                                      >  2:1220198/1‑150 (MQ=255)
                                                                                                                aaTATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCaaa                                     <  1:4578138/150‑1 (MQ=255)
                                                                                                                        aGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTAt                             >  2:1439122/1‑150 (MQ=255)
                                                                                                                          tGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATcc                           >  2:1817053/1‑150 (MQ=255)
                                                                                                                          tGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATcc                           <  2:3813416/150‑1 (MQ=255)
                                                                                                                            cAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGt                         <  2:2792702/150‑1 (MQ=255)
                                                                                                                            cAGATCTTTAATGCTGTCGATCATAGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGt                         >  1:4216228/1‑150 (MQ=255)
                                                                                                                             aGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGACACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGtt                        >  1:275564/1‑150 (MQ=255)
                                                                                                                             aGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGtt                        <  1:3159272/150‑1 (MQ=255)
                                                                                                                                 cTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCt                    <  2:4507914/150‑1 (MQ=255)
                                                                                                                                         cTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCatatat            <  2:3739058/150‑1 (MQ=255)
                                                                                                                                          tGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATg           <  1:860644/150‑1 (MQ=255)
                                                                                                                                           gTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATgg          <  1:116781/150‑1 (MQ=255)
                                                                                                                                            tCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACAGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGc         >  2:4278278/1‑150 (MQ=255)
                                                                                                                                             cGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCg        >  2:4174899/1‑150 (MQ=255)
                                                                                                                                              gATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGa       >  2:882863/1‑150 (MQ=255)
                                                                                                                                               atcatcGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGAc      <  2:4603039/150‑1 (MQ=255)
                                                                                                                                                 catcGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGAcgc    >  1:666091/1‑150 (MQ=255)
                                                                                                                                                  atcGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGAcgcg   <  2:1133220/150‑1 (MQ=255)
                                                                                                                                                   tcGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACgcgg  <  1:3514277/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGG  >  CP009273/602047‑602343

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: