Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 596,812 G→T 10.2% M881I (ATG→ATT)  cusA → copper/silver efflux system, membrane component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273596,8120GT10.2% 896.1 / 48.3 245M881I (ATG→ATT) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base G (128/92);  new base T (16/9);  total (144/101)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 254 total aligned reads displayed.

ACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTGGCTTTAT  >  CP009273/596664‑596961
                                                                                                                                                    |                                                                                                                                                     
     gtgtCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATc                                                                                                                                                 <  2:104970/150‑1 (MQ=255)
     gtgtCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATc                                                                                                                                                 <  1:1996673/150‑1 (MQ=255)
          ggtggtTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGt                                                                                                                                            <  2:4478251/150‑1 (MQ=255)
           gtggtTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTg                                                                                                                                           <  2:1138122/150‑1 (MQ=255)
             ggtTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCt                                                                                                                                         <  1:1364611/150‑1 (MQ=255)
               ttCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCtgt                                                                                                                                       >  2:3827217/1‑150 (MQ=255)
                 cACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCtgttg                                                                                                                                     <  1:3478798/150‑1 (MQ=255)
                    gATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTAt                                                                                                                                  >  2:4467220/1‑150 (MQ=255)
                    gATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTAt                                                                                                                                  >  1:3968694/1‑150 (MQ=255)
                     aTTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATc                                                                                                                                 >  1:4216643/1‑150 (MQ=255)
                       ttGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTACCTg                                                                                                                               <  2:1882820/150‑1 (MQ=255)
                          cAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCg                                                                                                                            >  2:2778610/1‑150 (MQ=255)
                           aaaatGCGATATCTTAAAATGTGCAGTTAAATCCGGGCACCAGCGTGGCATTCTACGGGGAGTTCGAGCTTCTGTAGCGTGCCAACCATAATCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGt                                                                                                                           <  2:3394779/150‑1 (MQ=255)
                           aaaaaGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGt                                                                                                                           >  2:2276996/1‑150 (MQ=255)
                            aaaaGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAACCTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGtt                                                                                                                          <  1:2582348/150‑1 (MQ=255)
                            aaaaGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGtt                                                                                                                          >  1:2682141/1‑150 (MQ=255)
                              aaGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTcc                                                                                                                        >  2:2648261/1‑150 (MQ=255)
                              aaGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTcc                                                                                                                        <  1:2480024/150‑1 (MQ=255)
                               aGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCg                                                                                                                       >  1:1222817/1‑150 (MQ=255)
                                   aTAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATTATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCgg                                                                                                                   >  2:2041960/1‑150 (MQ=255)
                                   aTAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCgg                                                                                                                   >  1:4618710/1‑150 (MQ=255)
                                    tAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCggg                                                                                                                  >  2:4382179/1‑150 (MQ=255)
                                     aGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGt                                                                                                                 >  1:3752405/1‑150 (MQ=255)
                                     aGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGt                                                                                                                 >  1:4196224/1‑150 (MQ=255)
                                       cTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTgg                                                                                                               <  2:674285/150‑1 (MQ=255)
                                         gAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGc                                                                                                             >  2:2413881/1‑150 (MQ=255)
                                         gAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGc                                                                                                             >  2:2541012/1‑150 (MQ=255)
                                         gAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGc                                                                                                             <  1:2275714/150‑1 (MQ=255)
                                          aaaaaGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCg                                                                                                            >  2:920749/1‑150 (MQ=255)
                                           aaaaGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATTATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGa                                                                                                           >  1:2762031/1‑150 (MQ=255)
                                            aaaGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGaa                                                                                                          <  2:4136162/150‑1 (MQ=255)
                                              aGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGc                                                                                                        >  2:3691128/1‑150 (MQ=255)
                                              aGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGc                                                                                                        >  1:1051399/1‑150 (MQ=255)
                                                 gCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGtt                                                                                                     >  2:4508506/1‑150 (MQ=255)
                                                  cAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTg                                                                                                    <  1:2758938/150‑1 (MQ=255)
                                                  cAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTAATTGAGCGCGCCAACCATATGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGATGCGTTg                                                                                                    <  1:325725/150‑1 (MQ=255)
                                                      tAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGa                                                                                                <  1:2414983/150‑1 (MQ=255)
                                                       aaaaCCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGAt                                                                                               >  2:1779734/1‑150 (MQ=255)
                                                        aaaCCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATTATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGAtt                                                                                              <  2:309698/150‑1 (MQ=255)
                                                        aaaCCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCGGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGAtt                                                                                              >  2:1926353/1‑150 (MQ=255)
                                                           ccGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATc                                                                                           <  2:4317808/150‑1 (MQ=255)
                                                           ccGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATc                                                                                           >  2:2992724/1‑150 (MQ=255)
                                                           ccGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATc                                                                                           >  1:3331682/1‑150 (MQ=255)
                                                             gggCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATcag                                                                                         >  1:611203/1‑150 (MQ=255)
                                                                 aCCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCTTCGGGTGGGCGAAGCGTTGCTGATTATCagcagc                                                                                     >  2:2870597/1‑150 (MQ=255)
                                                                 aCCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCagcagc                                                                                     <  2:2452269/150‑1 (MQ=255)
                                                                    aGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTa                                                                                  >  1:4118204/1‑150 (MQ=255)
                                                                       gTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCg                                                                               >  1:3815452/1‑150 (MQ=255)
                                                                       gTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCg                                                                               >  1:1303776/1‑150 (MQ=255)
                                                                          gCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATTATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGttt                                                                            <  1:335574/150‑1 (MQ=255)
                                                                           cATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTg                                                                           >  2:1773763/1‑150 (MQ=255)
                                                                           cATTCTCCGGGCAGTTCGAGCGACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTg                                                                           <  2:1697578/150‑1 (MQ=255)
                                                                                tcCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATTATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCtg                                                                      >  1:4235411/1‑150 (MQ=255)
                                                                                      cAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGggc                                                                <  2:19902/150‑1 (MQ=255)
                                                                                           cGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCAt                                                           >  1:2178961/1‑150 (MQ=255)
                                                                                            gAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATc                                                          >  1:567785/1‑150 (MQ=255)
                                                                                            gAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATc                                                          >  2:4394181/1‑150 (MQ=255)
                                                                                             aGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATTATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCt                                                         >  1:2298282/1‑150 (MQ=255)
                                                                                             aGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATTATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCt                                                         >  1:4136637/1‑150 (MQ=255)
                                                                                             aGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCt                                                         <  1:4382179/150‑1 (MQ=255)
                                                                                                tACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGt                                                      >  2:3604099/1‑150 (MQ=255)
                                                                                                    ggAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCt                                                  >  1:4083098/1‑150 (MQ=255)
                                                                                                      aGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCtgt                                                <  2:788884/150‑1 (MQ=255)
                                                                                                       gcgcgcCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCtgtg                                               <  2:3046962/150‑1 (MQ=255)
                                                                                                        cgcgcCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATTATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTgtgg                                              <  1:354701/150‑1 (MQ=255)
                                                                                                          cgcCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTgtggtg                                            >  2:1170513/1‑150 (MQ=255)
                                                                                                            ctaaCCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGa                                          >  1:3143362/3‑150 (MQ=255)
                                                                                                             cAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGAt                                         <  1:2921095/150‑1 (MQ=255)
                                                                                                             cAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGAt                                         >  1:4311358/1‑150 (MQ=255)
                                                                                                              aaCCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATg                                        <  2:4167143/150‑1 (MQ=255)
                                                                                                              aaCCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATg                                        <  1:1261188/150‑1 (MQ=255)
                                                                                                                ccATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATGATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATggg                                      <  2:1602981/150‑1 (MQ=255)
                                                                                                                    aaGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGTCttt                                  <  2:619836/150‑1 (MQ=255)
                                                                                                                     aGCTTAAACTCATGGTGCCGATGACGTTGATTATTATCTTCGTGCTGGTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGGGGTGGATGGGCTTTc                                 >  2:199999/1‑150 (MQ=255)
                                                                                                                     aGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTc                                 >  1:3992213/1‑150 (MQ=255)
                                                                                                                      gCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCa                                <  1:2990515/150‑1 (MQ=255)
                                                                                                                       cTTATAGTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGTGTGGGCGAAGCGTTGCTGATTATCAGCAGAGTACCGTTTGCGCGGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCAt                               <  2:2332178/150‑1 (MQ=255)
                                                                                                                        ttAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATc                              >  1:1496128/1‑150 (MQ=255)
                                                                                                                         tAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCt                             >  2:2224391/1‑150 (MQ=255)
                                                                                                                          aaaCTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCtt                            >  2:2242773/1‑150 (MQ=255)
                                                                                                                          aaaCTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCtt                            >  1:65425/1‑150 (MQ=255)
                                                                                                                           aaCTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCttt                           <  1:1220563/150‑1 (MQ=255)
                                                                                                                           aaCTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCttt                           >  1:569679/1‑150 (MQ=255)
                                                                                                                            aCTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTc                          >  2:3534588/1‑150 (MQ=255)
                                                                                                                             cTCATGGTGCCGATGAGGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTcc                         <  2:1542805/150‑1 (MQ=255)
                                                                                                                             cTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTcc                         <  2:1646383/150‑1 (MQ=255)
                                                                                                                                aTGGAGCCGATGACGTTGATTATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGAGTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTg                      <  2:1493225/150‑1 (MQ=255)
                                                                                                                                  ggTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGc                    <  2:4216530/150‑1 (MQ=255)
                                                                                                                                       cGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACgg               >  1:1623018/1‑150 (MQ=255)
                                                                                                                                          tGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCa            >  1:755242/1‑150 (MQ=255)
                                                                                                                                             cGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTg         <  2:3418917/150‑1 (MQ=255)
                                                                                                                                             cGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTg         <  1:1841733/150‑1 (MQ=255)
                                                                                                                                              gTTGATTATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTgg        >  1:3696499/1‑150 (MQ=255)
                                                                                                                                              gTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGTGACGGGCACTgg        >  1:1088142/1‑150 (MQ=255)
                                                                                                                                              gTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTgg        <  2:136466/150‑1 (MQ=255)
                                                                                                                                              gTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTgg        >  2:844953/1‑150 (MQ=255)
                                                                                                                                               ttGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTGGc       >  1:1734523/1‑150 (MQ=255)
                                                                                                                                                 gatgatTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTGGCtt     <  1:1914115/150‑1 (MQ=255)
                                                                                                                                                  atgatTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTGGCttt    >  1:1537099/1‑150 (MQ=255)
                                                                                                                                                    gatTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTGGCTTTAt  >  1:99802/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCATCTGGTTGCTGTGGTGGATGGGCTTTCATCTTTCCGTGGCGACGGGCACTGGCTTTAT  >  CP009273/596664‑596961

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