Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,936,672 C→A 15.4% intergenic (‑168/‑24) ygdI ← / → csdA DUF903 family verified lipoprotein/cysteine sulfinate desulfinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,936,6720CA15.4% 897.5 / 111.8 277intergenic (‑168/‑24)ygdI/csdADUF903 family verified lipoprotein/cysteine sulfinate desulfinase
Reads supporting (aligned to +/- strand):  ref base C (146/88);  new base A (19/24);  total (165/112)
Fisher's exact test for biased strand distribution p-value = 2.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 288 total aligned reads displayed.

TTAATAAATCTTAATAGTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTT  >  CP009273/2936524‑2936820
                                                                                                                                                    |                                                                                                                                                    
ttAATAAATCTTAATAGTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCCGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATc                                                                                                                                                      >  2:1327907/1‑149 (MQ=255)
ttAATAAATCTTAATAGTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATc                                                                                                                                                      <  2:726477/149‑1 (MQ=255)
ttAATAAATCTTAATAGTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATc                                                                                                                                                      <  1:518167/149‑1 (MQ=255)
  aataaATCTTAATAGTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCgg                                                                                                                                                    >  2:455513/1‑149 (MQ=255)
      aaTCTTAATAGTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCa                                                                                                                                                >  2:1006780/1‑149 (MQ=255)
         cTTAATAGTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCa                                                                                                                                             >  2:351546/1‑149 (MQ=255)
            aaTAGTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGc                                                                                                                                          <  1:1101804/149‑1 (MQ=255)
                gTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCgag                                                                                                                                      >  2:20546/1‑149 (MQ=255)
                 tATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCgagg                                                                                                                                     <  2:1871045/149‑1 (MQ=255)
                   ttGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCgaggag                                                                                                                                   <  2:1379404/149‑1 (MQ=255)
                   ttGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCgaggag                                                                                                                                   >  2:4019573/1‑149 (MQ=255)
                   ttGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCgaggag                                                                                                                                   >  1:279518/1‑149 (MQ=255)
                   ttGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCgaggag                                                                                                                                   >  1:280441/1‑149 (MQ=255)
                     gATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGTGTa                                                                                                                                 >  2:1721890/1‑149 (MQ=255)
                     gATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTa                                                                                                                                 >  1:1422249/1‑149 (MQ=255)
                     gATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTa                                                                                                                                 <  1:4677593/149‑1 (MQ=255)
                       ttAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTAcc                                                                                                                               >  2:734165/1‑149 (MQ=255)
                        tAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCa                                                                                                                              <  2:3203616/149‑1 (MQ=255)
                           cAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGa                                                                                                                           >  2:5215630/1‑149 (MQ=255)
                            aGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATAGGTGCATCAAGCCGAGGAGTACCATGaa                                                                                                                          <  1:4139332/149‑1 (MQ=255)
                            aGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATAGGTGCATCAAGCCGAGGAGTACCATGaa                                                                                                                          <  1:4139197/149‑1 (MQ=255)
                             ggCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAAc                                                                                                                         <  1:820403/149‑1 (MQ=255)
                                 aaaaTTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGttt                                                                                                                     <  1:4915831/149‑1 (MQ=255)
                                  aaaTTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGtttt                                                                                                                    >  2:3311557/1‑149 (MQ=255)
                                     ttAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTaa                                                                                                                 >  1:4508034/1‑149 (MQ=255)
                                        aCGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATAGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATcc                                                                                                              <  1:215012/149‑1 (MQ=255)
                                         cGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATccc                                                                                                             <  2:4157960/149‑1 (MQ=255)
                                          gCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCcg                                                                                                            >  1:4256066/1‑149 (MQ=255)
                                                  cTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGttt                                                                                                    <  1:4355750/149‑1 (MQ=255)
                                                    aTTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTcg                                                                                                  >  1:4632942/1‑149 (MQ=255)
                                                       cAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCgcgc                                                                                               <  1:4724301/149‑1 (MQ=255)
                                                       cAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCgcgc                                                                                               <  1:2111070/149‑1 (MQ=255)
                                                       cAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATAGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCgcgc                                                                                               >  2:492464/1‑149 (MQ=255)
                                                        atcATCTGTGACTTACGCGGCTTCGCCAGGATATCCAGAAAATTCGGATGGTTAGCACGCTCCTTGTATCAAAGTGAATTTTGCTTCTCGATCGGGGCATCAAGCCGTGGAGTACCATTACCGTTTTTAATCCCGCGCAGTTTCGCGcc                                                                                              <  1:224484/147‑1 (MQ=255)
                                                        aGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGcc                                                                                              >  2:1811598/1‑149 (MQ=255)
                                                         gCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGccc                                                                                             >  1:2892663/1‑149 (MQ=255)
                                                         gCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATAGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGccc                                                                                             <  1:484609/149‑1 (MQ=255)
                                                           atatGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAg                                                                                           >  2:111931/1‑149 (MQ=255)
                                                           atatGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATAGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAg                                                                                           >  2:3650581/1‑149 (MQ=255)
                                                            tatGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGt                                                                                          >  1:3297818/1‑149 (MQ=255)
                                                            tatGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGt                                                                                          >  1:3112058/1‑149 (MQ=255)
                                                             atGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGtt                                                                                         >  1:2074843/1‑149 (MQ=255)
                                                              tGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGttt                                                                                        >  1:3422043/1‑149 (MQ=255)
                                                              tGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGttt                                                                                        >  1:1693962/1‑149 (MQ=255)
                                                                ttacttacGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTcc                                                                                      <  1:2440241/149‑1 (MQ=255)
                                                                ttacttacGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTcc                                                                                      <  2:2371680/149‑1 (MQ=255)
                                                                ttacttacGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTcc                                                                                      >  2:472327/1‑149 (MQ=255)
                                                                ttacttacGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATAGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTcc                                                                                      <  2:3529139/149‑1 (MQ=255)
                                                                 tacttacGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTccc                                                                                     <  1:3610989/149‑1 (MQ=255)
                                                                   cttacGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGc                                                                                   >  2:1774181/1‑149 (MQ=255)
                                                                    ttacGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCa                                                                                  >  2:3135882/1‑149 (MQ=255)
                                                                      acGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACt                                                                                <  1:2313542/149‑1 (MQ=255)
                                                                        gcgGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTAc                                                                              >  1:926476/1‑149 (MQ=255)
                                                                           gCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAgg                                                                           >  1:670053/1‑149 (MQ=255)
                                                                            cTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGa                                                                          <  2:2790333/149‑1 (MQ=255)
                                                                            cTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGa                                                                          <  2:827209/149‑1 (MQ=255)
                                                                                gCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCg                                                                      >  1:1294236/1‑149 (MQ=255)
                                                                                 ccAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCgg                                                                     >  1:1282602/1‑149 (MQ=255)
                                                                                 ccAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCgg                                                                     >  1:4326191/1‑149 (MQ=255)
                                                                                  cAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATAGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCggg                                                                    <  1:2334945/149‑1 (MQ=255)
                                                                                     gATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGt                                                                 >  1:3171720/1‑149 (MQ=255)
                                                                                        atCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTa                                                              >  2:548454/1‑149 (MQ=255)
                                                                                         tCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTAt                                                             >  2:5115915/1‑149 (MQ=255)
                                                                                           cAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTAtct                                                           <  2:76451/149‑1 (MQ=255)
                                                                                           cAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTAtct                                                           <  2:626083/149‑1 (MQ=255)
                                                                                            aGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTAtctc                                                          >  2:4131366/1‑149 (MQ=255)
                                                                                              aTAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGa                                                        <  1:261002/149‑1 (MQ=255)
                                                                                              aTAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGa                                                        >  1:3991955/1‑149 (MQ=255)
                                                                                                    cTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAgcgc                                                  <  2:2217223/149‑1 (MQ=255)
                                                                                                         ggTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGa                                             <  1:1328138/149‑1 (MQ=255)
                                                                                                          gTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGAc                                            <  1:4440090/149‑1 (MQ=255)
                                                                                                          gTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGAc                                            <  1:3599076/149‑1 (MQ=255)
                                                                                                          gTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGAc                                            >  1:2856606/1‑149 (MQ=255)
                                                                                                           ttAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATAGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGAcc                                           >  1:611985/1‑149 (MQ=255)
                                                                                                            tAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATAGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACcg                                          >  1:3010063/1‑149 (MQ=255)
                                                                                                              gCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACcgcg                                        >  1:2224716/1‑149 (MQ=255)
                                                                                                                aCTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCt                                      <  1:1006759/149‑1 (MQ=255)
                                                                                                                 ctctCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCtt                                     <  2:3839010/149‑1 (MQ=255)
                                                                                                                 ctctCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCtt                                     >  1:5082013/1‑149 (MQ=255)
                                                                                                                  tctcCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTa                                    >  1:1847899/1‑149 (MQ=255)
                                                                                                                      cTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAAcc                                <  1:603180/149‑1 (MQ=255)
                                                                                                                        tGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTg                              <  1:3887646/149‑1 (MQ=255)
                                                                                                                         gTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGa                             <  2:273761/149‑1 (MQ=255)
                                                                                                                           aTCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAg                           >  2:3383104/1‑149 (MQ=255)
                                                                                                                              aaaGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCg                        >  2:4366925/1‑149 (MQ=255)
                                                                                                                              aaaGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCg                        >  1:3742988/1‑149 (MQ=255)
                                                                                                                               aaGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGt                       >  1:2528381/1‑149 (MQ=255)
                                                                                                                                aGTGAATTTTGCGTCACGATAGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTg                      >  1:270769/1‑149 (MQ=255)
                                                                                                                                 gTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTgg                     <  2:4666252/149‑1 (MQ=255)
                                                                                                                                 gTGAATTTTGCGTCACGATAGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTgg                     <  1:2994652/149‑1 (MQ=255)
                                                                                                                                  tGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGt                    >  2:2296741/1‑149 (MQ=255)
                                                                                                                                   gAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGtt                   <  1:758072/149‑1 (MQ=255)
                                                                                                                                   gAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGtt                   <  2:303723/149‑1 (MQ=255)
                                                                                                                                       tttGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAg               <  1:455513/149‑1 (MQ=255)
                                                                                                                                        ttGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGc              <  2:1584123/149‑1 (MQ=255)
                                                                                                                                              cACGATAGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCa        >  1:653068/1‑149 (MQ=255)
                                                                                                                                               aCGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCaa       >  2:2650610/1‑149 (MQ=255)
                                                                                                                                                  aTCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAg    <  2:2489996/149‑1 (MQ=255)
                                                                                                                                                    cGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGtt  <  1:3265728/149‑1 (MQ=255)
                                                                                                                                                    cGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAAACTGAATCCGTGGTTGAAGCCACCCAACAGtt  >  2:4066657/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTAATAAATCTTAATAGTATTGATTAACAGGCTAAAATTAACGCCTAACACTATTCAGCATATGTTACTTACGCGGCTTCGCCAGGATATCCAGATAATTCTGATGGTTAGCACTCTCCTTGTATCAAAGTGAATTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTT  >  CP009273/2936524‑2936820

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: