Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,127,022 T→G 100% E188D (GAA→GAC yijF ← DUF1287 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,127,0220TG96.2% 1009.2 / 11.9 257E188D (GAA→GACyijFDUF1287 family protein
Reads supporting (aligned to +/- strand):  ref base T (6/4);  new base G (137/110);  total (143/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 262 total aligned reads displayed.

GATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTT  >  CP009273/4126875‑4127168
                                                                                                                                                   |                                                                                                                                                  
gagtgtgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACAtcgtc                                                                                                                                                  <  2:1428826/147‑1 (MQ=255)
 attgtgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACAtcttct                                                                                                                                                 >  2:4199042/1‑150 (MQ=255)
   tgtgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCgtcttc                                                                                                                                               >  2:3968353/1‑150 (MQ=255)
   tgtgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCgtcttc                                                                                                                                               >  2:1205373/1‑150 (MQ=255)
    gtgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTcc                                                                                                                                              <  2:1574203/150‑1 (MQ=255)
    gtgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTcc                                                                                                                                              >  2:4404249/1‑150 (MQ=255)
     tgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCt                                                                                                                                             <  1:463158/150‑1 (MQ=255)
              tGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTg                                                                                                                                    <  2:2716903/150‑1 (MQ=255)
               ggCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTgc                                                                                                                                   <  1:495602/150‑1 (MQ=255)
               ggCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTgc                                                                                                                                   >  2:3185885/1‑150 (MQ=255)
                     ggTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAAt                                                                                                                             >  1:1362264/1‑150 (MQ=255)
                       tGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGt                                                                                                                           <  1:939471/150‑1 (MQ=255)
                        gAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGtt                                                                                                                          <  2:663506/150‑1 (MQ=255)
                        gAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGtt                                                                                                                          <  1:711176/150‑1 (MQ=255)
                            cGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGt                                                                                                                      <  2:923006/150‑1 (MQ=255)
                             gTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTa                                                                                                                     <  2:783467/150‑1 (MQ=255)
                              tAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTAt                                                                                                                    <  2:4300417/150‑1 (MQ=255)
                               aaGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATc                                                                                                                   >  2:3499967/1‑150 (MQ=255)
                                 gAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCAc                                                                                                                 >  1:4316597/1‑150 (MQ=255)
                                 gAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCAc                                                                                                                 >  1:1479263/1‑150 (MQ=255)
                                   aaaaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCa                                                                                                               <  1:3706975/150‑1 (MQ=255)
                                    aaaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAg                                                                                                              >  1:3321064/1‑150 (MQ=255)
                                      aagaagCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAgcg                                                                                                            <  1:738608/150‑1 (MQ=255)
                                      aagaagCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAgcg                                                                                                            >  2:1029460/1‑150 (MQ=255)
                                      aagaagCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAgcg                                                                                                            <  1:2896043/150‑1 (MQ=255)
                                            ccGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTAc                                                                                                      >  2:1527532/1‑150 (MQ=255)
                                              gTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCg                                                                                                    >  1:808802/1‑150 (MQ=255)
                                              gTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCg                                                                                                    >  1:3787188/1‑150 (MQ=255)
                                                   tgtAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTcgcg                                                                                               >  2:2571637/1‑150 (MQ=255)
                                                       aacaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCg                                                                                           <  2:876982/147‑1 (MQ=255)
                                                       aaaaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCg                                                                                           <  1:3015745/150‑1 (MQ=255)
                                                        aaaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGa                                                                                          <  1:1183576/150‑1 (MQ=255)
                                                        aaaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGa                                                                                          >  2:4310782/1‑150 (MQ=255)
                                                           ggATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGc                                                                                       >  2:2801006/1‑150 (MQ=255)
                                                           ggATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGc                                                                                       >  2:4624798/1‑150 (MQ=255)
                                                              tCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTACCAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCAt                                                                                    <  1:4243644/150‑1 (MQ=255)
                                                                tGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCt                                                                                  <  2:2994490/150‑1 (MQ=255)
                                                                 gAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAAGGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTg                                                                                 >  2:2756370/1‑150 (MQ=255)
                                                                  aaTTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGa                                                                                >  2:3096879/1‑150 (MQ=255)
                                                                         tcttcATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAaccacc                                                                         >  2:41734/1‑150 (MQ=255)
                                                                         tcttcATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAaccacc                                                                         <  1:669800/150‑1 (MQ=255)
                                                                          cttcATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCAccc                                                                        >  1:861623/1‑150 (MQ=255)
                                                                             cATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCaa                                                                     >  2:137716/1‑150 (MQ=255)
                                                                             cATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCaa                                                                     >  1:3675619/1‑150 (MQ=255)
                                                                              aTGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAat                                                                    >  2:629758/1‑150 (MQ=255)
                                                                                 tgGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATAtg                                                                 <  1:2374015/150‑1 (MQ=255)
                                                                                  ggggCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATAtgt                                                                <  2:3367626/150‑1 (MQ=255)
                                                                                  ggggCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATAtgt                                                                >  2:844367/1‑150 (MQ=255)
                                                                                    ggCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGc                                                              >  1:189889/1‑150 (MQ=255)
                                                                                       gATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAg                                                           <  2:2222393/150‑1 (MQ=255)
                                                                                           tctTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCg                                                       <  1:3283625/150‑1 (MQ=255)
                                                                                            ctTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGt                                                      >  2:4502082/1‑150 (MQ=255)
                                                                                            ctTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGt                                                      <  2:4428098/150‑1 (MQ=255)
                                                                                             ttATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGtt                                                     <  2:3502937/150‑1 (MQ=255)
                                                                                             ttATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGtt                                                     <  2:1957344/150‑1 (MQ=255)
                                                                                               aTTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTAt                                                   >  2:2005615/1‑150 (MQ=255)
                                                                                                tttAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATc                                                  >  1:1038026/1‑150 (MQ=255)
                                                                                                   aaCAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAg                                               >  1:177292/1‑150 (MQ=255)
                                                                                                   aaCAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAg                                               >  1:598531/1‑150 (MQ=255)
                                                                                                    aCAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGt                                              <  1:2329972/150‑1 (MQ=255)
                                                                                                     cAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTc                                             <  1:2620413/150‑1 (MQ=255)
                                                                                                     cAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTc                                             <  1:826316/150‑1 (MQ=255)
                                                                                                      aaaaTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCg                                            >  2:3640851/1‑150 (MQ=255)
                                                                                                      aaaaTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCg                                            >  2:2898520/1‑150 (MQ=255)
                                                                                                       aaaTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGc                                           >  2:3074765/1‑150 (MQ=255)
                                                                                                           aaCGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATTTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCACAATCTGTGCCAGCCCGTTATCCAGTCGCCAgg                                       <  1:3070343/150‑1 (MQ=255)
                                                                                                             cGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgag                                     >  1:4453324/1‑150 (MQ=255)
                                                                                                             cGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgag                                     >  1:4453198/1‑150 (MQ=255)
                                                                                                             cGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgag                                     >  1:4009094/1‑150 (MQ=255)
                                                                                                              gATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgaga                                    <  1:2229353/150‑1 (MQ=255)
                                                                                                               aTAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGAc                                   >  2:3417149/1‑150 (MQ=255)
                                                                                                                tAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACa                                  >  2:1012770/1‑150 (MQ=255)
                                                                                                                tAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACa                                  >  2:1181691/1‑150 (MQ=255)
                                                                                                                 aaTGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACaa                                 <  2:1114667/150‑1 (MQ=255)
                                                                                                                  aTGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAat                                >  1:1003109/1‑150 (MQ=255)
                                                                                                                    gCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAatat                              >  1:166811/1‑150 (MQ=255)
                                                                                                                        caccaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGcc                          >  2:538823/1‑150 (MQ=255)
                                                                                                                          ccaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCg                        <  1:1565274/150‑1 (MQ=255)
                                                                                                                           caTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGc                       >  2:1659252/1‑150 (MQ=255)
                                                                                                                            aTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCt                      >  2:1341205/1‑150 (MQ=255)
                                                                                                                            aTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCt                      <  1:2512608/150‑1 (MQ=255)
                                                                                                                             tCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCtt                     <  1:74839/150‑1 (MQ=255)
                                                                                                                              ccGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTg                    <  1:2286677/150‑1 (MQ=255)
                                                                                                                               cGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGa                   >  2:75855/1‑150 (MQ=255)
                                                                                                                                gCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGAt                  >  1:3858634/1‑150 (MQ=255)
                                                                                                                                gCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGAt                  <  2:1881132/150‑1 (MQ=255)
                                                                                                                                gCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGAt                  <  2:1329982/150‑1 (MQ=255)
                                                                                                                                 ccAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATa                 <  1:29694/150‑1 (MQ=255)
                                                                                                                                    gTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTc              >  1:2292658/1‑150 (MQ=255)
                                                                                                                                     ttAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCa             >  1:3213849/1‑150 (MQ=255)
                                                                                                                                      tAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCAc            >  1:3928432/1‑150 (MQ=255)
                                                                                                                                      tAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCAc            <  2:2399631/150‑1 (MQ=255)
                                                                                                                                          cAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTggg        >  2:1795222/1‑150 (MQ=255)
                                                                                                                                            gCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCCGGAGACACTATCGCCCGCTTGATAGTCACTGGGGt      >  2:4608013/1‑150 (MQ=255)
                                                                                                                                             cacaTCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGtt     <  2:717313/150‑1 (MQ=255)
                                                                                                                                              acaTCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTc    >  1:1893575/1‑150 (MQ=255)
                                                                                                                                               caTCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCt   >  1:3401942/1‑150 (MQ=255)
                                                                                                                                                aTCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCtt  <  2:3265458/150‑1 (MQ=255)
                                                                                                                                                aTCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCtt  <  1:2678763/150‑1 (MQ=255)
                                                                                                                                                aTCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCtt  >  2:374006/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTT  >  CP009273/4126875‑4127168

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: