Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 337,720 G→T 100% E295D (GAG→GAT yahJ → putative metallo‑dependent hydrolase domain deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273337,7200GT96.1% 478.3 / 7.5 165E295D (GAG→GATyahJputative metallo‑dependent hydrolase domain deaminase
Reads supporting (aligned to +/- strand):  ref base G (1/6);  new base T (80/78);  total (81/84)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 167 total aligned reads displayed.

CGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTG  >  CP009273/337575‑337864
                                                                                                                                                 |                                                                                                                                                
cGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATaa                                                                                                                                                >  1:700922/1‑148 (MQ=255)
 gACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATaaa                                                                                                                                               <  1:2645384/148‑1 (MQ=255)
 gACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATaaa                                                                                                                                               <  1:2646458/148‑1 (MQ=255)
  aCCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATaaaa                                                                                                                                              >  1:3525782/1‑148 (MQ=255)
    cAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACg                                                                                                                                            <  2:804274/148‑1 (MQ=255)
      gtgtTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGcc                                                                                                                                          <  2:208944/148‑1 (MQ=255)
      gtgtTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGAGATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGcc                                                                                                                                          <  2:2651857/148‑1 (MQ=255)
       tgtTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCca                                                                                                                                         >  2:1636273/1‑148 (MQ=255)
          tGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAg                                                                                                                                      >  2:1407459/1‑148 (MQ=255)
            aTGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCt                                                                                                                                    <  2:1816255/148‑1 (MQ=255)
            aTGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCt                                                                                                                                    >  2:1856394/1‑148 (MQ=255)
             tGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTg                                                                                                                                   <  2:853070/148‑1 (MQ=255)
              ggCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGa                                                                                                                                  <  1:461379/148‑1 (MQ=255)
               gcgcGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAACGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGaa                                                                                                                                 <  2:1495536/148‑1 (MQ=255)
               gcgcGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACTAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGaa                                                                                                                                 <  1:2914500/148‑1 (MQ=255)
               gcgcGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGaa                                                                                                                                 <  2:1495527/148‑1 (MQ=255)
                 gcgATGTAAAAATCCCTCGACACCTTGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAgg                                                                                                                               <  1:1896055/148‑1 (MQ=255)
                 gcgATGGAAAAAGCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATAAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAgg                                                                                                                               <  1:1898349/148‑1 (MQ=255)
                    aTGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCa                                                                                                                            >  1:130448/1‑148 (MQ=255)
                      ggAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCCCCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAg                                                                                                                          <  1:1115910/148‑1 (MQ=255)
                      ggAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAg                                                                                                                          >  1:3388704/1‑148 (MQ=255)
                      ggAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAg                                                                                                                          <  2:1142896/148‑1 (MQ=255)
                          aaaTCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGa                                                                                                                      >  2:460718/1‑148 (MQ=255)
                           aaTCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGAc                                                                                                                     >  2:753913/1‑148 (MQ=255)
                                cTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCa                                                                                                                >  2:3340235/1‑148 (MQ=255)
                                cTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCa                                                                                                                >  1:2743259/1‑148 (MQ=255)
                                cTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCa                                                                                                                >  2:1571607/1‑148 (MQ=255)
                                  cGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGt                                                                                                              <  1:2154434/148‑1 (MQ=255)
                                  cGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGt                                                                                                              >  1:1758643/1‑148 (MQ=255)
                                  cGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGt                                                                                                              >  1:2957556/1‑148 (MQ=255)
                                     cacCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCAc                                                                                                           >  2:1706496/1‑148 (MQ=255)
                                      acCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATCTGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACg                                                                                                          >  2:3483658/1‑148 (MQ=255)
                                       ccATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTGGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGc                                                                                                         >  1:25537/1‑148 (MQ=255)
                                         aTGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCt                                                                                                       >  1:977989/1‑148 (MQ=255)
                                         aTGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCt                                                                                                       <  1:1891711/148‑1 (MQ=255)
                                         aTGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCt                                                                                                       <  2:1482362/148‑1 (MQ=255)
                                          tGTTCCAGATTGCGCTGGACTACGACACAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCtt                                                                                                      <  2:2377974/148‑1 (MQ=255)
                                          tGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCtt                                                                                                      <  2:1924762/148‑1 (MQ=255)
                                          tGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCtt                                                                                                      <  1:2337051/148‑1 (MQ=255)
                                            ttCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTg                                                                                                    >  1:1787843/1‑148 (MQ=255)
                                               cAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAACACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTgcgc                                                                                                 <  1:862909/148‑1 (MQ=255)
                                               cAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTgcgc                                                                                                 <  2:432275/148‑1 (MQ=255)
                                               cAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTgcgc                                                                                                 <  1:761677/148‑1 (MQ=255)
                                                        cTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGc                                                                                        >  1:303170/1‑148 (MQ=255)
                                                           gACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCa                                                                                     <  2:700922/148‑1 (MQ=255)
                                                           gACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCa                                                                                     >  2:3062275/1‑148 (MQ=255)
                                                           gACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCa                                                                                     <  1:1317846/148‑1 (MQ=255)
                                                              tACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACg                                                                                  >  2:2967221/1‑148 (MQ=255)
                                                              tACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACg                                                                                  >  2:2968465/1‑148 (MQ=255)
                                                               aCGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGa                                                                                 >  1:2236069/1‑148 (MQ=255)
                                                                 gACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGc                                                                               <  1:481033/148‑1 (MQ=255)
                                                                  aCAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCa                                                                              >  1:1615078/1‑148 (MQ=255)
                                                                     aaGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACa                                                                           >  1:3296492/1‑148 (MQ=255)
                                                                      aGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAg                                                                          <  1:2252849/148‑1 (MQ=255)
                                                                      aGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAg                                                                          <  2:2482587/148‑1 (MQ=255)
                                                                       ggCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAgg                                                                         >  2:1772832/1‑148 (MQ=255)
                                                                         cGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTa                                                                       >  2:2703618/1‑148 (MQ=255)
                                                                          gTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAg                                                                      <  1:1024584/148‑1 (MQ=255)
                                                                           tCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGa                                                                     >  2:616444/1‑148 (MQ=255)
                                                                           tCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGa                                                                     <  2:3533367/148‑1 (MQ=255)
                                                                             gATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATg                                                                   <  2:3234632/148‑1 (MQ=255)
                                                                              atatTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTTCGCTGGCAACGCTCAACGAGCAACAGGTAGATGa                                                                  <  2:57698/148‑1 (MQ=255)
                                                                                      cTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGa                                                          <  1:2026871/148‑1 (MQ=255)
                                                                                       tGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGaa                                                         <  2:224167/148‑1 (MQ=255)
                                                                                          aCGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCg                                                      >  1:435469/1‑148 (MQ=255)
                                                                                           cGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCgg                                                     >  2:238292/1‑148 (MQ=255)
                                                                                           cGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCgg                                                     >  2:2164907/1‑148 (MQ=255)
                                                                                              aaCCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtg                                                  >  2:1935987/1‑148 (MQ=255)
                                                                                                ccACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtggt                                                >  1:2867130/1‑148 (MQ=255)
                                                                                                    tCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGgtggtg                                            <  1:3419300/148‑1 (MQ=255)
                                                                                                       ggCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGcaa                                         <  1:52733/148‑1 (MQ=255)
                                                                                                         cAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGcaaca                                       >  1:192048/1‑148 (MQ=255)
                                                                                                           ggCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCGGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACaaa                                     <  2:2102513/148‑1 (MQ=255)
                                                                                                            gCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAt                                    >  2:2335493/1‑148 (MQ=255)
                                                                                                            gCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAt                                    >  1:2364827/1‑148 (MQ=255)
                                                                                                            gCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAt                                    >  1:2090177/1‑148 (MQ=255)
                                                                                                               tGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTc                                 >  2:161248/1‑148 (MQ=255)
                                                                                                                ggCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCt                                <  2:2816395/148‑1 (MQ=255)
                                                                                                                 gCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTa                               >  2:1321884/1‑148 (MQ=255)
                                                                                                                     ccATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGc                           >  2:2658633/1‑148 (MQ=255)
                                                                                                                     ccATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGc                           <  1:2173866/148‑1 (MQ=255)
                                                                                                                     ccATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGc                           >  2:3274068/1‑148 (MQ=255)
                                                                                                                     acataatttATATGGTTTAAACTGTAGATAAAACGACACAGCTGAAGGGAAATCTGACAATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGATCATCAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGc                           <  1:2173592/141‑1 (MQ=255)
                                                                                                                      cATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGcc                          >  2:1147705/1‑148 (MQ=255)
                                                                                                                      cATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGcc                          >  2:1084809/1‑148 (MQ=255)
                                                                                                                      cATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGcc                          >  2:2864342/1‑148 (MQ=255)
                                                                                                                         cAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCg                       >  1:1086177/1‑148 (MQ=255)
                                                                                                                         cAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCg                       >  1:3125996/1‑148 (MQ=255)
                                                                                                                         cAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCg                       >  1:1743001/1‑148 (MQ=255)
                                                                                                                             tataTGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACgg                   >  2:3212414/1‑148 (MQ=255)
                                                                                                                               tatGTTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTg                 <  1:1408444/148‑1 (MQ=255)
                                                                                                                               tatGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTg                 >  2:945373/1‑148 (MQ=255)
                                                                                                                               tatGGTTGAAAATGTAGATAATACGCAACAGCTGAAGGGCAAGCTGACCATCAGTCTCGCCTTTGCGCTGTCAACTCTCAAAGAGCAACAGGTAGATGAAATGGCGAACCGGATGGTTGTTCAACATATTTCTATCGCCTCGACGGTg                 <  2:2179101/148‑1 (MQ=255)
                                                                                                                                    ttGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGAt            >  1:420279/1‑148 (MQ=255)
                                                                                                                                     tGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGAtt           <  2:472525/148‑1 (MQ=255)
                                                                                                                                     tGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGAtt           >  1:3088240/1‑148 (MQ=255)
                                                                                                                                     tGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGAtt           <  1:1856241/148‑1 (MQ=255)
                                                                                                                                        aaCGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGc        >  1:1618739/1‑148 (MQ=255)
                                                                                                                                             tAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCt   >  2:3155759/1‑148 (MQ=255)
                                                                                                                                              agatAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTg  <  2:2024502/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTG  >  CP009273/337575‑337864

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: