Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,346,500 G→T 11.2% N645K (AAC→AAA cadA ← lysine decarboxylase, acid‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,346,5000GT11.2% 805.3 / 51.2 223N645K (AAC→AAAcadAlysine decarboxylase, acid‑inducible
Reads supporting (aligned to +/- strand):  ref base G (97/101);  new base T (13/12);  total (110/113)
Fisher's exact test for biased strand distribution p-value = 8.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 224 total aligned reads displayed.

CGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATTGTGGTGAACAATCAG  >  CP009273/4346352‑4346649
                                                                                                                                                    |                                                                                                                                                     
cGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGt                                                                                                                                                      <  1:3569189/150‑1 (MQ=255)
cGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGt                                                                                                                                                      <  1:2370817/150‑1 (MQ=255)
    aaTATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAAt                                                                                                                                                  <  2:3278621/150‑1 (MQ=255)
     atatCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATa                                                                                                                                                 <  1:3173108/150‑1 (MQ=255)
      tatCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATAc                                                                                                                                                <  2:3556788/150‑1 (MQ=255)
      tatCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATAc                                                                                                                                                >  2:2249071/1‑150 (MQ=255)
       atCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACg                                                                                                                                               <  1:513221/150‑1 (MQ=255)
           gTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCt                                                                                                                                           <  2:3243030/150‑1 (MQ=255)
             ttCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTAc                                                                                                                                         <  2:3309036/150‑1 (MQ=255)
              tCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTAcc                                                                                                                                        <  1:3389944/150‑1 (MQ=255)
              tCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTAcc                                                                                                                                        <  1:2171255/150‑1 (MQ=255)
                aaaGCCCGGATGGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCAt                                                                                                                                      >  2:3839668/1‑150 (MQ=255)
                aaaGCCCGGATGGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCAt                                                                                                                                      >  2:213121/1‑150 (MQ=255)
                aaaGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCAt                                                                                                                                      >  2:3644511/1‑150 (MQ=255)
                 aaGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCAtt                                                                                                                                     >  1:4499023/1‑150 (MQ=255)
                  aGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCTTTAATACGACCTACCAttt                                                                                                                                    >  1:1073552/1‑150 (MQ=255)
                    cccGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTtcg                                                                                                                                  <  1:1176836/150‑1 (MQ=255)
                     ccGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTtcgt                                                                                                                                 >  2:3662899/1‑150 (MQ=255)
                      cGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTtcgtc                                                                                                                                <  2:3449779/150‑1 (MQ=255)
                       ggATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCTTTAATACGACCTACCATTtcgtcg                                                                                                                               <  1:4071322/150‑1 (MQ=255)
                           aGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCTTTAATACGACCTACCATTTCGTCGAGGt                                                                                                                           <  2:4384889/150‑1 (MQ=255)
                               aGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAAc                                                                                                                       >  1:2299559/1‑150 (MQ=255)
                                gcgcCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAAct                                                                                                                      <  1:2512107/150‑1 (MQ=255)
                                  gcCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAActtc                                                                                                                    >  1:440650/1‑150 (MQ=255)
                                  gcCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAActtc                                                                                                                    >  1:3968917/1‑150 (MQ=255)
                                  gcCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAActtc                                                                                                                    <  2:1088720/150‑1 (MQ=255)
                                   ccGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAActtct                                                                                                                   <  1:954641/150‑1 (MQ=255)
                                    cGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCTTTAATACGACCTACCATTTCGTCGAGGTAAActtctt                                                                                                                  <  1:748809/150‑1 (MQ=255)
                                       tttCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCgg                                                                                                               >  1:4051314/1‑150 (MQ=255)
                                       tttCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCgg                                                                                                               >  1:2192308/1‑150 (MQ=255)
                                        ttCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGt                                                                                                              <  2:4006322/150‑1 (MQ=255)
                                        ttCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGt                                                                                                              >  2:519899/1‑150 (MQ=255)
                                         tCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTc                                                                                                             <  2:1986243/150‑1 (MQ=255)
                                          cacacaGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCa                                                                                                            <  1:1009181/150‑1 (MQ=255)
                                           acacaGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCAt                                                                                                           >  1:4600723/1‑150 (MQ=255)
                                             acaGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATAc                                                                                                         <  2:1051585/150‑1 (MQ=255)
                                              caGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATAcc                                                                                                        <  1:977225/150‑1 (MQ=255)
                                                    cTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAg                                                                                                  >  2:2781220/1‑150 (MQ=255)
                                                    cTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAg                                                                                                  >  2:1120754/1‑150 (MQ=255)
                                                    cTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAg                                                                                                  >  1:1809842/1‑150 (MQ=255)
                                                    cTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAg                                                                                                  >  2:2000962/1‑150 (MQ=255)
                                                         ggAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCtt                                                                                             >  1:2601757/1‑150 (MQ=255)
                                                              tCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTgg                                                                                        <  1:2327562/150‑1 (MQ=255)
                                                                   aaCAGGACGGCTTTATTCGGTGATCATTTCAACCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTTCCATTTCGTCGAGGTTAACTTCTTCGGTCATACCGTGCAGATCTTTCTGGAATgc                                                                                   <  2:1310959/150‑1 (MQ=255)
                                                                    aCCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGACTgca                                                                                  >  1:2406133/1‑150 (MQ=255)
                                                                    aCCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATgca                                                                                  >  1:1306568/1‑150 (MQ=255)
                                                                     ccGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATgcag                                                                                 <  2:4488660/150‑1 (MQ=255)
                                                                      cGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATgcagc                                                                                >  1:1925913/1‑150 (MQ=255)
                                                                       ggACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATgcagca                                                                               >  2:1454786/1‑150 (MQ=255)
                                                                       ggACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATgcagca                                                                               >  2:3745276/1‑150 (MQ=255)
                                                                        gACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCAt                                                                              <  1:997544/150‑1 (MQ=255)
                                                                          cGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATAc                                                                            <  1:4145487/150‑1 (MQ=255)
                                                                             cTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGa                                                                         >  1:3113765/1‑150 (MQ=255)
                                                                             cTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGa                                                                         <  2:1380960/150‑1 (MQ=255)
                                                                               ttcttcGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTACTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGt                                                                       <  1:4261502/150‑1 (MQ=255)
                                                                                 cttcGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCa                                                                     >  1:1580289/1‑150 (MQ=255)
                                                                                    cGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCTTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTa                                                                  >  2:916680/1‑150 (MQ=255)
                                                                                     ggTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTAc                                                                 <  1:1574543/150‑1 (MQ=255)
                                                                                      gTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTATCTGGAATGCAGCATACGGAGTCATTAcc                                                                >  1:1884545/1‑150 (MQ=255)
                                                                                         aTCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCAtc                                                             <  1:3051442/150‑1 (MQ=255)
                                                                                         aTCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCAtc                                                             <  2:785545/150‑1 (MQ=255)
                                                                                          tCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCAtcg                                                            <  1:923172/150‑1 (MQ=255)
                                                                                          tCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCAtcg                                                            >  2:712995/1‑150 (MQ=255)
                                                                                               tCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGgc                                                       <  2:3391772/150‑1 (MQ=255)
                                                                                                   ccGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCTTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGgcagca                                                   >  2:735181/1‑150 (MQ=255)
                                                                                                      gCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACtt                                                <  2:4490435/150‑1 (MQ=255)
                                                                                                      gCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATAAGGAGTCATTACCATCGTCGGCAGCACtt                                                <  1:3682334/150‑1 (MQ=255)
                                                                                                          tACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCaaa                                            >  2:3111165/1‑150 (MQ=255)
                                                                                                           aCCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAAt                                           <  2:498317/150‑1 (MQ=255)
                                                                                                              agagGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATgcg                                        >  2:1671964/1‑150 (MQ=255)
                                                                                                                 ggAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGa                                     >  1:2152880/1‑150 (MQ=255)
                                                                                                                 ggAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGa                                     >  1:2152918/1‑150 (MQ=255)
                                                                                                                  gAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATATCGTCGAGGTAAACTTCTTCGGTCATACCGTGCATCTCTTTCTGGAATGCAGCATACGGAGTCATTAACATCGTCTGCAGCACTTCAAATGCGCGAt                                    <  2:1267865/150‑1 (MQ=255)
                                                                                                                   aaCTCCCGGCGGGTACGGAAGGATCATATTGGCTTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATa                                   <  2:330895/150‑1 (MQ=255)
                                                                                                                             gggTACGGAAGGATCATATTGGCTTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATcc                         >  1:123196/1‑150 (MQ=255)
                                                                                                                              ggTACGGAAGGATCATATTGGCTTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCg                        <  2:159973/150‑1 (MQ=255)
                                                                                                                               gTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCgg                       >  1:4578801/1‑150 (MQ=255)
                                                                                                                                tACGGAAGGATCATATTGGCTTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGc                      <  2:3846826/150‑1 (MQ=255)
                                                                                                                                tACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGc                      <  2:2586856/150‑1 (MQ=255)
                                                                                                                                  cGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAg                    <  1:3801772/150‑1 (MQ=255)
                                                                                                                                   ggAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGa                   >  1:1881085/1‑150 (MQ=255)
                                                                                                                                   ggAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGa                   >  1:4230582/1‑150 (MQ=255)
                                                                                                                                   ggAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGa                   >  1:256974/1‑150 (MQ=255)
                                                                                                                                     aaGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGAtt                 <  2:1359452/150‑1 (MQ=255)
                                                                                                                                       ggATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATtgt               <  2:804953/150‑1 (MQ=255)
                                                                                                                                       ggATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATtgt               <  2:4146789/150‑1 (MQ=255)
                                                                                                                                        gATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATtgtg              >  1:457900/1‑150 (MQ=255)
                                                                                                                                         aTCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATTgtgg             >  1:935636/1‑150 (MQ=255)
                                                                                                                                          tCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATTgtggt            <  1:32873/150‑1 (MQ=255)
                                                                                                                                           cATATTGGCTTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATTgtggtg           >  2:2828089/1‑150 (MQ=255)
                                                                                                                                           cATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATTgtggtg           <  2:1248387/150‑1 (MQ=255)
                                                                                                                                           cATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATTgtggtg           <  1:2844640/150‑1 (MQ=255)
                                                                                                                                            atatTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATTGTGGTGa          >  1:1379379/1‑150 (MQ=255)
                                                                                                                                                tGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATTGTGGTGAACaa      >  2:4568029/1‑150 (MQ=255)
                                                                                                                                                 ggCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATTGTGGTGAACAAt     <  1:1674619/150‑1 (MQ=255)
                                                                                                                                                 ggCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATTGTGGTGAACAAt     >  2:1834300/1‑150 (MQ=255)
                                                                                                                                                 ggCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATTGTGGTGAACAAt     >  2:1834415/1‑150 (MQ=255)
                                                                                                                                                  gCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATTGTGGTGAACAATc    <  1:888807/150‑1 (MQ=255)
                                                                                                                                                  gCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATTGTGGTGAACAATc    <  2:921965/150‑1 (MQ=255)
                                                                                                                                                   cGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATTGTGGTGAACAATCa   <  2:3631154/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGACCTACCATTTCGTCGAGGTAAACTTCTTCGGTCATACCGTGCAGCTCTTTCTGGAATGCAGCATACGGAGTCATTACCATCGTCGGCAGCACTTCAAATGCGCGATACATCAGATCCGGCAGATTGTGGTGAACAATCAG  >  CP009273/4346352‑4346649

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: