Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 378,179 A→T 100% pseudogene (150/282 nt) yaiX ← pseudogene, interrupted by IS2A, acetyltransferase homolog; nonfunctional; interruped by IS2; putative transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273378,1790AT87.8% 832.0 / 90.4 279pseudogene (150/282 nt)yaiXpseudogene, interrupted by IS2A, acetyltransferase homolog; nonfunctional; interruped by IS2; putative transferase
Reads supporting (aligned to +/- strand):  ref base A (14/20);  new base T (103/142);  total (117/162)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 296 total aligned reads displayed.

ATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATAAGTC  >  CP009273/378032‑378325
                                                                                                                                                   |                                                                                                                                                  
 taGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAAtttt                                                                                                                                                  <  1:2814771/149‑1 (MQ=255)
  aGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATtttt                                                                                                                                                 >  2:3259349/1‑149 (MQ=255)
    gggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATtttttt                                                                                                                                               <  1:1982698/149‑1 (MQ=255)
     ggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTtc                                                                                                                                              <  1:505929/149‑1 (MQ=255)
     ggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTtc                                                                                                                                              <  2:4085021/149‑1 (MQ=255)
       cAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTtctc                                                                                                                                            <  2:103789/149‑1 (MQ=255)
       cAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTtctc                                                                                                                                            <  2:1601827/149‑1 (MQ=255)
       cAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAAGAATTACCTGGTCTGCAATCTGAATTTTTTtctc                                                                                                                                            >  2:1395666/1‑149 (MQ=255)
        aaaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTctct                                                                                                                                           <  1:2951700/149‑1 (MQ=255)
        aaaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTctct                                                                                                                                           <  1:632660/149‑1 (MQ=255)
        aaaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTctct                                                                                                                                           <  2:2138350/149‑1 (MQ=255)
         aaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCtt                                                                                                                                          <  1:4690245/149‑1 (MQ=255)
          aTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCttt                                                                                                                                         >  2:628088/1‑149 (MQ=255)
            ccAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTaa                                                                                                                                       <  2:2672000/149‑1 (MQ=255)
             cAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAAc                                                                                                                                      >  1:4124910/1‑149 (MQ=255)
                ttACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACAtt                                                                                                                                   >  2:2805716/1‑149 (MQ=255)
                 tACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTa                                                                                                                                  >  2:1955626/1‑149 (MQ=255)
                 tACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTa                                                                                                                                  <  2:4621835/149‑1 (MQ=255)
                   ccTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAg                                                                                                                                <  1:3093032/149‑1 (MQ=255)
                    cTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGc                                                                                                                               >  2:3515730/1‑149 (MQ=255)
                     tATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCt                                                                                                                              <  1:4350875/149‑1 (MQ=255)
                      aTCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCtt                                                                                                                             <  1:873290/149‑1 (MQ=255)
                      aTCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCtt                                                                                                                             >  1:411415/1‑149 (MQ=255)
                       tCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTc                                                                                                                            <  2:2470060/149‑1 (MQ=255)
                         aGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAg                                                                                                                          <  2:3117928/149‑1 (MQ=255)
                          ggCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGa                                                                                                                         <  2:4005169/149‑1 (MQ=255)
                          ggCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGa                                                                                                                         >  2:3809681/1‑149 (MQ=255)
                          ggCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGa                                                                                                                         <  1:4297444/149‑1 (MQ=255)
                          ggCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGa                                                                                                                         <  1:4003325/149‑1 (MQ=255)
                            cAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATa                                                                                                                       <  2:787565/149‑1 (MQ=255)
                             aGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAg                                                                                                                      >  2:2049804/1‑149 (MQ=255)
                               tttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGtt                                                                                                                    >  1:2351207/1‑149 (MQ=255)
                                ttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGttt                                                                                                                   <  1:3425890/149‑1 (MQ=255)
                                ttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGttt                                                                                                                   <  1:779060/149‑1 (MQ=255)
                                   cgcCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTc                                                                                                                <  1:4668669/149‑1 (MQ=255)
                                   cgcCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTc                                                                                                                <  2:1430162/149‑1 (MQ=255)
                                   cgcCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTc                                                                                                                <  2:1781400/149‑1 (MQ=255)
                                     ccAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAg                                                                                                              <  1:2575319/149‑1 (MQ=255)
                                     ccAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAg                                                                                                              >  2:1244761/1‑149 (MQ=255)
                                     ccAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGAGAGTTTTTCAg                                                                                                              >  2:4461079/1‑149 (MQ=255)
                                        aTCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTc                                                                                                           <  1:1843095/149‑1 (MQ=255)
                                         tCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTcc                                                                                                          <  1:4480510/149‑1 (MQ=255)
                                         tCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTcc                                                                                                          <  1:3010413/149‑1 (MQ=255)
                                         tCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTcc                                                                                                          <  1:3742198/149‑1 (MQ=255)
                                               ccGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGGTTCATCAATACTTACCTGGTCTGCACTCTGAATTTTTTTCTCTTTAACATTAAGCTGCAGATAGTTTTTCAGTTCCTGTGAc                                                                                                    >  2:4488343/1‑149 (MQ=255)
                                                  gACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGa                                                                                                 <  2:726540/149‑1 (MQ=255)
                                                   aCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGaa                                                                                                <  2:3059110/149‑1 (MQ=255)
                                                   aCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGaa                                                                                                <  2:262616/149‑1 (MQ=255)
                                                   aCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGaa                                                                                                >  2:2283530/1‑149 (MQ=255)
                                                           tAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGa                                                                                        >  1:1419279/1‑149 (MQ=255)
                                                             aCGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAAt                                                                                      >  2:4256393/1‑149 (MQ=255)
                                                               ggCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATaa                                                                                    <  2:4371322/149‑1 (MQ=255)
                                                               ggCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATaa                                                                                    <  2:3005115/149‑1 (MQ=255)
                                                               ggCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATaa                                                                                    >  2:3555961/1‑149 (MQ=255)
                                                                gCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAg                                                                                   <  1:1220291/149‑1 (MQ=255)
                                                                 cACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGc                                                                                  <  2:3748612/149‑1 (MQ=255)
                                                                   ccATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCgg                                                                                <  1:4612546/149‑1 (MQ=255)
                                                                    cATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGt                                                                               <  2:682166/149‑1 (MQ=255)
                                                                     aTGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGtt                                                                              >  1:2743407/1‑149 (MQ=255)
                                                                      tGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTg                                                                             <  1:3135292/149‑1 (MQ=255)
                                                                      tGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTg                                                                             <  1:2194113/149‑1 (MQ=255)
                                                                         cAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACg                                                                          >  1:2093748/1‑149 (MQ=255)
                                                                         cAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACg                                                                          >  2:2952278/1‑149 (MQ=255)
                                                                             tACGGGTATTCGCGCCGATAACAACTTCACCTGCAGGTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTc                                                                      <  1:2283530/149‑1 (MQ=255)
                                                                                   tATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTAACTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATAAAGATTCAGATAGATTTTAAGTTACTGTTACATAAGGTTGGGAATAAGCGGTTGTCTCTTCCACCAg                                                                >  2:1191626/1‑149 (MQ=255)
                                                                                     ttCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTg                                                              >  1:541543/1‑149 (MQ=255)
                                                                                        gcgcCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCgg                                                           >  1:1763631/1‑149 (MQ=255)
                                                                                          gcCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGaa                                                         >  1:4528127/1‑149 (MQ=255)
                                                                                           ccGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAAt                                                        <  2:2382347/149‑1 (MQ=255)
                                                                                            cGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATa                                                       <  2:2462510/149‑1 (MQ=255)
                                                                                             gATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAg                                                      <  2:3458411/149‑1 (MQ=255)
                                                                                              aTAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGc                                                     >  1:2614633/1‑149 (MQ=255)
                                                                                                aacaacTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgg                                                   >  1:3067097/1‑149 (MQ=255)
                                                                                                aacaacTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgg                                                   >  1:51857/1‑149 (MQ=255)
                                                                                                     cTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgggtggg                                              >  1:3699137/1‑149 (MQ=255)
                                                                                                      ttCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgggtgggt                                             >  2:1327021/1‑149 (MQ=255)
                                                                                                        cACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGttt                                           <  1:2910981/149‑1 (MQ=255)
                                                                                                        cACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGttt                                           >  1:3359560/1‑149 (MQ=255)
                                                                                                         aCCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTg                                          >  1:3607786/1‑149 (MQ=255)
                                                                                                             gCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTcgc                                      >  1:3538690/1‑149 (MQ=255)
                                                                                                               aGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCgcgc                                    <  1:99509/149‑1 (MQ=255)
                                                                                                                gTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCt                                   >  1:3159636/1‑149 (MQ=255)
                                                                                                                    catcaATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAgg                               <  1:2329415/149‑1 (MQ=255)
                                                                                                                     atcaATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGt                              >  1:3847322/1‑149 (MQ=255)
                                                                                                                     atcaATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGt                              <  1:1187410/149‑1 (MQ=255)
                                                                                                                      tcaATAATTACCGGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGtt                             <  1:1744086/149‑1 (MQ=255)
                                                                                                                             ttACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTcc                      <  2:2602336/149‑1 (MQ=255)
                                                                                                                             ttACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTcc                      <  1:1591039/149‑1 (MQ=255)
                                                                                                                             ttACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTcc                      <  2:2727398/149‑1 (MQ=255)
                                                                                                                              tACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGTAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCa                     >  1:2413276/1‑149 (MQ=255)
                                                                                                                              tACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCa                     >  1:3230074/1‑149 (MQ=255)
                                                                                                                                   ggTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTaaa                <  2:4197217/149‑1 (MQ=255)
                                                                                                                                     tCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTaaaaa              >  1:3281796/1‑149 (MQ=255)
                                                                                                                                       tGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAg            >  2:289794/1‑149 (MQ=255)
                                                                                                                                         cAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGc          <  2:1833775/149‑1 (MQ=255)
                                                                                                                                         cAATCAGAATTTTTTTCTCTTGAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGc          <  2:4433661/149‑1 (MQ=255)
                                                                                                                                           aTCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCaa        >  2:1221453/1‑149 (MQ=255)
                                                                                                                                             cTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATa      <  2:3650100/149‑1 (MQ=255)
                                                                                                                                             cTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATa      <  2:3212991/149‑1 (MQ=255)
                                                                                                                                              tGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATaa     >  2:2468309/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATAAGTC  >  CP009273/378032‑378325

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: