Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 504,372 T→C 4.5% S822G (AGT→GGT)  copA ← copper transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273504,3720TC 4.5% 309.4 / 3.6 110S822G (AGT→GGT) copAcopper transporter
Reads supporting (aligned to +/- strand):  ref base T (53/52);  new base C (3/2);  total (56/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-01

Only 100 of 112 total aligned reads displayed.

AACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTG  >  CP009273/504224‑504520
                                                                                                                                                    |                                                                                                                                                    
aaCCTGTGCCTGAGCCGTAGGTCGAATAAGGCGTTCACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCGAAAGCAGCGCATCCGCAATGATGTACTAATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACt                                                                                                                                                      <  1:633361/149‑1 (MQ=39)
  ccTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCa                                                                                                                                                    >  2:937338/1‑149 (MQ=255)
  ccTGGGCCTGACCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGCCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCa                                                                                                                                                    <  2:359554/149‑1 (MQ=255)
    tgtgCCTAAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTGGGCGTTACTCACt                                                                                                                                                  >  2:819167/1‑149 (MQ=255)
     ttgCCCGAACCGTAGGCCGGTTAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTa                                                                                                                                                 <  2:1107748/148‑1 (MQ=255)
     gtgCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGTAGCAACCGGTTGGCGTTACTCACTa                                                                                                                                                 <  1:740848/149‑1 (MQ=255)
       gCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCAGTCACGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATAATGTACTTATTCCGTCGGTTTAACCCGCAGCAACCGGTTGACGTTGCTCAcagcg                                                                                                                                               >  2:96800/1‑145 (MQ=38)
         cTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTGTTCCTTCGGTTAAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGt                                                                                                                                             >  2:309630/1‑149 (MQ=255)
         cTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGACAGTCATGCGTCGAGGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAGTGATGTAGTTATTCCTTCGGTTTACACCCCAGCCACCGGTTGGCGTTACCCCCTCCGGt                                                                                                                                             >  2:101859/1‑149 (MQ=38)
          tGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGGTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTa                                                                                                                                            >  2:648869/1‑149 (MQ=255)
          tGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGAGAAGGTTCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTa                                                                                                                                            >  1:551258/1‑149 (MQ=39)
             cccGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATc                                                                                                                                         <  2:1162472/148‑1 (MQ=255)
             aCCGTGCGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATc                                                                                                                                         <  2:108780/149‑1 (MQ=255)
             aCCGTAGGCCGGATAAGGCGTTGACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATc                                                                                                                                         <  1:367038/149‑1 (MQ=255)
              ccGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCCTGCGTCGATGCCAAATGCGCCAACCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTGGTTTTAAACCGCAGCAACCGGTTGGCGTTACCCACTCCGTGAATcg                                                                                                                                        >  2:476621/1‑149 (MQ=38)
              ccGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACCCACTACGGTAATcg                                                                                                                                        >  2:505745/1‑149 (MQ=255)
               cGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTCAACTGCAGCAACCGGTTGGCGTTACTTATTACGGTAATcga                                                                                                                                       >  1:87826/1‑149 (MQ=255)
               cGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATcga                                                                                                                                       <  1:499910/149‑1 (MQ=255)
                gTAGGCCGGATACGGCGTGTACGCCGCCTCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCTTCCGCAATGATGTACTCATTCCTTCGGTTTAAACCGCAGCAACCTGTTCGCGTTACTCACTACGGTAATcgac                                                                                                                                      <  1:926335/149‑1 (MQ=255)
                 tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATcgacg                                                                                                                                     >  2:553905/1‑149 (MQ=255)
                  aGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATcgacga                                                                                                                                    <  2:517232/149‑1 (MQ=255)
                    gCCGGATAAGGGGTTTACGCCGCACCCGGCAGTCATGCGGCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGAAATGATGTACTTATTCTTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACgaga                                                                                                                                  <  2:305358/149‑1 (MQ=255)
                    gCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCCTCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAGCCGGTTGGCGCTACTCACTACGGTAGTcgacgaaa                                                                                                                                  >  2:601055/1‑147 (MQ=255)
                    gCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACgagg                                                                                                                                  >  1:313441/1‑148 (MQ=255)
                    gCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACgaga                                                                                                                                  >  1:446962/1‑149 (MQ=255)
                    gCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCACAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACgaga                                                                                                                                  >  1:1193162/1‑149 (MQ=255)
                      cGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACCACGGTAATCGACGAGAgc                                                                                                                                >  1:1084067/1‑149 (MQ=255)
                              gCGTTTACGCCGCATCCGGCAGTCATGCGTCGAGGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGc                                                                                                                        >  1:907196/1‑149 (MQ=255)
                              gCGATTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCGTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCCCTACGGTAATCGACGCGACTGCCCTTGc                                                                                                                        >  2:498568/1‑149 (MQ=255)
                               cGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGcc                                                                                                                       >  1:447223/1‑149 (MQ=255)
                                gTTTCCGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTGAAGCAGCGCATCCGCAATGATGTTCTTTTTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCTCTCCGGTAATCGACGAGAGCGCCATTGCCg                                                                                                                      <  2:256032/149‑1 (MQ=255)
                                gTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCGGAAACCGGTTGGCGTTACTCACTACGGTAATCGACgagcgggcacttcgcg                                                                                                                      >  2:1183037/1‑136 (MQ=255)
                                   tACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTc                                                                                                                   >  1:1003885/1‑149 (MQ=255)
                                   tACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTc                                                                                                                   >  1:178953/1‑149 (MQ=255)
                                    aCGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTcc                                                                                                                  <  2:153637/149‑1 (MQ=255)
                                        cgcATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCgggc                                                                                                              >  2:432123/1‑147 (MQ=255)
                                             ccGGCTGTCTGGCGTCGGTGCCAAATGCCCCACCCTAGCGCACCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAGCCGCAGCCACCGGTTGGTGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTAc                                                                                                         <  2:1208683/149‑1 (MQ=255)
                                             ccGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTAc                                                                                                         <  1:766073/149‑1 (MQ=255)
                                               ggtaGTCATGCGTCGAAGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCTCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCg                                                                                                       <  2:831787/146‑1 (MQ=255)
                                               ggCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCg                                                                                                       <  1:855880/149‑1 (MQ=255)
                                                      aTGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGGAATCGAGGAGAGCGCCATTGCCGCTCCGGCAACACCCGGGGTAAg                                                                                                >  2:1121776/1‑149 (MQ=255)
                                                       tGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGc                                                                                               <  1:837309/149‑1 (MQ=255)
                                                       cgCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAAGCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGc                                                                                               >  2:908454/2‑149 (MQ=255)
                                                                 ccAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAg                                                                                     <  2:1003885/149‑1 (MQ=255)
                                                                 ccAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAg                                                                                     <  2:178953/149‑1 (MQ=255)
                                                                  cAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGt                                                                                    >  1:598757/1‑149 (MQ=255)
                                                                   aaaTGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTg                                                                                   >  1:388969/1‑149 (MQ=255)
                                                                     aTGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGaa                                                                                 <  1:1071843/149‑1 (MQ=255)
                                                                        cgcCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACgg                                                                              <  1:1069599/149‑1 (MQ=255)
                                                                        cgcCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACgg                                                                              <  1:438631/149‑1 (MQ=255)
                                                                        cgcCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACgg                                                                              <  1:250249/149‑1 (MQ=255)
                                                                         gcCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGc                                                                             <  2:884255/149‑1 (MQ=255)
                                                                           cACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGccc                                                                           >  1:584905/1‑148 (MQ=255)
                                                                           cACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCca                                                                           >  1:19691/1‑149 (MQ=255)
                                                                            aCCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCcac                                                                          <  2:446962/149‑1 (MQ=255)
                                                                             cccTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCccca                                                                         >  1:139751/1‑149 (MQ=255)
                                                                              ccTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACaa                                                                        >  1:152296/1‑149 (MQ=255)
                                                                              ccTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACaa                                                                        >  1:245216/1‑149 (MQ=255)
                                                                               cTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTGTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACaaa                                                                       <  2:396847/149‑1 (MQ=255)
                                                                                tAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGCCGGGACTACCTGGTCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACaaaa                                                                      >  1:595673/1‑149 (MQ=255)
                                                                                tAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCCGTGTTCCAGTGCACGGCCACaaaa                                                                      >  2:91932/1‑149 (MQ=255)
                                                                                tAAAGCAGCGCATCCGCAATGATGTACTTATGCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACaaaa                                                                      <  2:71354/149‑1 (MQ=255)
                                                                                  aaGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATa                                                                    <  2:981526/149‑1 (MQ=255)
                                                                                   agcagcGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATAc                                                                   <  2:1209242/149‑1 (MQ=255)
                                                                                      gtcgcATCCGCAATGATTTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACcgg                                                                <  2:252590/147‑1 (MQ=255)
                                                                                        cgcATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACcggcg                                                              >  2:841408/1‑149 (MQ=255)
                                                                                         gcATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACcggcgg                                                             <  2:752610/149‑1 (MQ=255)
                                                                                         gcATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACcggcgg                                                             <  2:752764/149‑1 (MQ=255)
                                                                                         gcATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACcggcgg                                                             <  1:1230054/149‑1 (MQ=255)
                                                                                         gcATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACcggcgg                                                             <  2:1048395/149‑1 (MQ=255)
                                                                                         gcACCCCCAACGATGTACTTACTCCTTCGGTTTAAACCGCCGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACcggcgg                                                             <  1:383547/149‑1 (MQ=255)
                                                                                          cATCCGCAATGATGTGCTTATTCCTTCGGTTTAAACCGCAGCGACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCggcggc                                                            <  1:726655/149‑1 (MQ=255)
                                                                                            tCCGCAATGGTGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGGTAAGCAGTGTTCCAGTGAACGGCCACAAAAAACCggcggcac                                                          >  2:1089500/1‑147 (MQ=255)
                                                                                               gCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCg                                                       >  1:623984/1‑149 (MQ=255)
                                                                                                  atgatgTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACCCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGaa                                                    <  1:505745/149‑1 (MQ=255)
                                                                                                     atgTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTGACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATAc                                                 <  1:432123/149‑1 (MQ=255)
                                                                                                       gTACTGATCCCTTCGGTGTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCg                                               <  2:233025/149‑1 (MQ=255)
                                                                                                        tACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGa                                              >  2:432637/1‑149 (MQ=255)
                                                                                                          cTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGAAACGAAATACCGGCGGCGACCGGAATACCGATa                                            <  2:882708/149‑1 (MQ=255)
                                                                                                                       tttAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAACATACCGGCGGCGACCGGAATACCGATACTGTAGTAGATaa                               >  2:1120650/1‑149 (MQ=255)
                                                                                                                        ttAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATaaa                              <  1:311891/149‑1 (MQ=255)
                                                                                                                            aCCGCAGCAACCGGTTGGCGTTACTACCTACGATAAACGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGCATACCGATACGGTTGGAGATAAACGcc                          >  2:268811/1‑148 (MQ=255)
                                                                                                                              cGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAACTACCGGCGGCGACCGGCATACCGATACTGTCGTAGATAACAGCAcc                        >  2:1017423/1‑149 (MQ=255)
                                                                                                                              cGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTGCATAAACGcccc                        >  2:187863/1‑149 (MQ=255)
                                                                                                                                   caACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGTACGGCCACAAAATACCGGCGGCGACCGGAATACTGATACTGTTGTAGATAAACGCACCGAGCa                   <  1:806357/149‑1 (MQ=255)
                                                                                                                                   caACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACGACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCa                   <  1:185231/149‑1 (MQ=255)
                                                                                                                                     aCCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAgg                 <  2:388969/149‑1 (MQ=255)
                                                                                                                                        ggTTGGCGGTACTCACTCCGGTACTCGACGAGAGCGCCATGGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTc              <  2:949626/149‑1 (MQ=255)
                                                                                                                                            gtcGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCATTGTTCCAGTGAACGGTCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCt          <  1:1117186/147‑1 (MQ=255)
                                                                                                                                            ggCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCCCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGAAACAGTTGTAGAAAAACGCACCGAGCAGGTTCTGCt          >  2:1063238/1‑149 (MQ=255)
                                                                                                                                             gCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCtt         <  1:657897/149‑1 (MQ=255)
                                                                                                                                               gTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCa       <  1:432637/149‑1 (MQ=255)
                                                                                                                                               gTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACAGGGACGGGTCTGCTTCa       >  1:300979/1‑149 (MQ=255)
                                                                                                                                               gTTACTCACTACGGTAATCGACGAGAGAGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGAAACAGTTGTAAATAAACGAACCGAGCAGGTTCTGCTCCa       >  2:386020/1‑149 (MQ=255)
                                                                                                                                                ttACTCCCTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCAt      <  2:355315/149‑1 (MQ=255)
                                                                                                                                                ttACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCAt      >  2:717598/1‑149 (MQ=255)
                                                                                                                                                 tACTCACTCCGGCATTCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAAGACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCAtg     <  2:758059/149‑1 (MQ=255)
                                                                                                                                                  aCTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCAtgt    >  1:81010/1‑149 (MQ=255)
                                                                                                                                                  aCTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCAtgt    >  2:1019279/1‑149 (MQ=255)
                                                                                                                                                  aCTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTAAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCAtgt    <  2:133741/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTG  >  CP009273/504224‑504520

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