Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,461,726 C→A 6.3% D163Y (GAC→TAC)  pyrB ← aspartate carbamoyltransferase, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,461,7260CA 6.3% 885.3 / 30.3 252D163Y (GAC→TAC) pyrBaspartate carbamoyltransferase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (124/112);  new base A (11/5);  total (135/117)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 255 total aligned reads displayed.

CAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCG  >  CP009273/4461577‑4461875
                                                                                                                                                     |                                                                                                                                                     
    cTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCacca                                                                                                                                                   >  1:9503/1‑150 (MQ=255)
      ccATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCaccaac                                                                                                                                                 >  1:1387312/1‑150 (MQ=255)
           cGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTg                                                                                                                                            >  2:5037664/1‑150 (MQ=255)
           cGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTg                                                                                                                                            >  2:1672048/1‑150 (MQ=255)
            gATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGc                                                                                                                                           >  2:2242089/1‑150 (MQ=255)
             aTCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCg                                                                                                                                          <  1:4698690/150‑1 (MQ=255)
                ccTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACg                                                                                                                                       >  2:3745949/1‑150 (MQ=255)
                ccTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACg                                                                                                                                       >  1:174834/1‑150 (MQ=255)
                 cTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGt                                                                                                                                      <  2:2098104/150‑1 (MQ=255)
                 cTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGt                                                                                                                                      >  2:187388/1‑150 (MQ=255)
                   ttttCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTgg                                                                                                                                    >  2:127224/1‑150 (MQ=255)
                     ttCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGgag                                                                                                                                  <  1:3081850/150‑1 (MQ=255)
                     ttCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGgag                                                                                                                                  <  1:3082038/150‑1 (MQ=255)
                     ttCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGgag                                                                                                                                  <  1:4320430/150‑1 (MQ=255)
                      tCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTAACCAACCATTGCGACGTGgaga                                                                                                                                 >  2:410966/1‑150 (MQ=255)
                         tCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTg                                                                                                                              <  1:1580843/150‑1 (MQ=255)
                             gCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCACCCATTGCGACGTGGAGATTGTCCa                                                                                                                          <  2:457803/150‑1 (MQ=255)
                             gCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCa                                                                                                                          <  2:3388477/150‑1 (MQ=255)
                              cATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAg                                                                                                                         <  1:541189/150‑1 (MQ=255)
                              cATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAg                                                                                                                         >  1:1235185/1‑150 (MQ=255)
                               atatCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGa                                                                                                                        >  1:2059275/1‑150 (MQ=255)
                                  tCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGc                                                                                                                     >  2:2470836/1‑150 (MQ=255)
                                   ccAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGcc                                                                                                                    >  1:2247067/1‑150 (MQ=255)
                                     aGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGcccc                                                                                                                  >  1:3112172/1‑150 (MQ=255)
                                        aTGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTgg                                                                                                               <  1:3983319/150‑1 (MQ=255)
                                        aTGTAATGCTGCATTGCCAGCGCGTCCCCCCCCATGAAGTAAAAACGGTTGCCGTCGAACTTCGATAACGCCTCACTCACGGAGTGAACGGTGCGGCCAAAAATCAGGTCACCAACCATAGCGACGACGAGCTTGTCCAGACTCCCCttg                                                                                                               >  2:2230478/1‑148 (MQ=255)
                                             ttGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTc                                                                                                          >  2:3999188/1‑150 (MQ=255)
                                               gCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCt                                                                                                        <  2:3673293/150‑1 (MQ=255)
                                                   cATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAAt                                                                                                    <  2:2300462/150‑1 (MQ=255)
                                                   cATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAAt                                                                                                    <  2:2701138/150‑1 (MQ=255)
                                                    aTTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATa                                                                                                   <  1:100212/150‑1 (MQ=255)
                                                    aTTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATa                                                                                                   >  2:3442405/1‑150 (MQ=255)
                                                    aTTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATa                                                                                                   <  2:4534194/150‑1 (MQ=255)
                                                    aTTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATa                                                                                                   >  1:3519536/1‑150 (MQ=255)
                                                    aTTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATa                                                                                                   <  1:100152/150‑1 (MQ=255)
                                                     gtGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAg                                                                                                  <  1:3403604/149‑1 (MQ=255)
                                                       gCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTg                                                                                                >  2:3148092/1‑150 (MQ=255)
                                                         cAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGaa                                                                                              <  1:3135997/150‑1 (MQ=255)
                                                         cAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGaa                                                                                              <  1:4462386/150‑1 (MQ=255)
                                                         cAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGaa                                                                                              <  1:2886032/150‑1 (MQ=255)
                                                          aGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAAt                                                                                             <  1:3162571/150‑1 (MQ=255)
                                                           gcgcGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGACTCGTGAATa                                                                                            >  2:5078701/1‑150 (MQ=255)
                                                            cgcgTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTAACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATaa                                                                                           <  1:1972744/150‑1 (MQ=255)
                                                                  cGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTAACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAg                                                                                     >  2:3813992/1‑150 (MQ=255)
                                                                   ggCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGc                                                                                    >  1:1445011/1‑150 (MQ=255)
                                                                     cgcgATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCaa                                                                                  >  1:2223276/1‑150 (MQ=255)
                                                                      gcgATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAg                                                                                 >  2:2478435/1‑150 (MQ=255)
                                                                       cgATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAgg                                                                                >  1:4344439/1‑150 (MQ=255)
                                                                           gAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTg                                                                            <  2:2976471/150‑1 (MQ=255)
                                                                             aGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCg                                                                          >  1:1870528/1‑150 (MQ=255)
                                                                              gTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGt                                                                         >  1:5037313/1‑150 (MQ=255)
                                                                              gTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGt                                                                         <  2:4402858/150‑1 (MQ=255)
                                                                              gTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCAATGCGACGTGTAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGt                                                                         >  1:5036746/1‑150 (MQ=255)
                                                                                aaaaaCGGGTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCg                                                                       <  1:751725/150‑1 (MQ=255)
                                                                                     cGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATgt                                                                  >  2:1491598/1‑150 (MQ=255)
                                                                                         tGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATgttggt                                                              <  1:1333605/150‑1 (MQ=255)
                                                                                           ccGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATgttggttg                                                            >  2:4026295/1‑150 (MQ=255)
                                                                                           ccGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATgttggttg                                                            >  1:3417482/1‑150 (MQ=255)
                                                                                            cGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGTTTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGttggttgg                                                           <  1:1315713/150‑1 (MQ=255)
                                                                                            cGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGttggttgg                                                           <  2:1869061/150‑1 (MQ=255)
                                                                                             gTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGa                                                          <  2:1388329/150‑1 (MQ=255)
                                                                                              tCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAg                                                         >  1:354645/1‑150 (MQ=255)
                                                                                                 aaCTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCa                                                      <  2:3875494/150‑1 (MQ=255)
                                                                                                 aaCTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCa                                                      >  1:2650496/1‑150 (MQ=255)
                                                                                                 aaCTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCa                                                      <  1:4946998/150‑1 (MQ=255)
                                                                                                     tCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCAc                                                  >  2:3898715/1‑150 (MQ=255)
                                                                                                       gCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCg                                                <  1:2837545/150‑1 (MQ=255)
                                                                                                        cTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCgg                                               >  1:2132087/1‑150 (MQ=255)
                                                                                                        cTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCgg                                               >  1:658081/1‑150 (MQ=255)
                                                                                                        cTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCgg                                               <  2:5181124/150‑1 (MQ=255)
                                                                                                           aCGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCAt                                            <  2:3315575/150‑1 (MQ=255)
                                                                                                           aCGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCAt                                            >  1:3584787/1‑150 (MQ=255)
                                                                                                            cGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCAtt                                           <  1:2460998/150‑1 (MQ=255)
                                                                                                              ccTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCa                                         <  2:256729/150‑1 (MQ=255)
                                                                                                                      aGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGt                                 >  1:1670217/1‑150 (MQ=255)
                                                                                                                      aGGGAGTGAACGGTGCGGCCATATTTCAGGTAACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGt                                 >  2:688479/1‑150 (MQ=255)
                                                                                                                         gAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACa                              >  2:3380363/1‑150 (MQ=255)
                                                                                                                         gAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACa                              >  2:4810642/1‑150 (MQ=255)
                                                                                                                          aGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACAt                             <  1:4517671/150‑1 (MQ=255)
                                                                                                                          aGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACAt                             <  2:4716887/150‑1 (MQ=255)
                                                                                                                           gTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACAtt                            <  2:3972398/150‑1 (MQ=255)
                                                                                                                             gAACGGTGCGGCCATATTTCAGGTAACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGc                          >  2:819904/1‑150 (MQ=255)
                                                                                                                                cGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCgg                       <  2:1978352/150‑1 (MQ=255)
                                                                                                                                     cGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAAc                  >  1:3815791/1‑150 (MQ=255)
                                                                                                                                      ggCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACt                 <  2:1522233/150‑1 (MQ=255)
                                                                                                                                          atatTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGt             <  1:788356/150‑1 (MQ=255)
                                                                                                                                            atTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTgg           >  1:5169974/1‑150 (MQ=255)
                                                                                                                                            atTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTgg           >  1:4892584/1‑150 (MQ=255)
                                                                                                                                             tttCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGc          <  2:2328152/150‑1 (MQ=255)
                                                                                                                                             tttCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGc          <  2:3515673/150‑1 (MQ=255)
                                                                                                                                              ttCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGcc         >  1:1684518/1‑150 (MQ=255)
                                                                                                                                               tCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCa        <  1:2809902/150‑1 (MQ=255)
                                                                                                                                               tCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCa        >  2:18593/1‑150 (MQ=255)
                                                                                                                                               tCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCa        <  2:1657275/150‑1 (MQ=255)
                                                                                                                                                cAGGTAACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAg       <  2:3203109/150‑1 (MQ=255)
                                                                                                                                                cAGGTAACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAg       <  2:596446/150‑1 (MQ=255)
                                                                                                                                                 aGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAgg      <  1:5088185/150‑1 (MQ=255)
                                                                                                                                                  ggTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGgc     >  2:1907953/1‑150 (MQ=255)
                                                                                                                                                    tCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGgcgc   >  1:3340787/1‑150 (MQ=255)
                                                                                                                                                     cACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGgcgcg  >  2:5113440/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGACTCCATGCGATCCCTTTTTCATCGAGCATATCCAGAATGTATTGCGGCATTGCCAGCGCGTCCGGCGCGATGAAGTAAAAACGGTTGCCGTCGAACTTCGCTAACGCCTGAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCG  >  CP009273/4461577‑4461875

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: