Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,329,774 C→A 85.6% D189Y (GAC→TAC)  adiA ← arginine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,329,7740CA85.6% 386.8 / 58.3 184D189Y (GAC→TAC) adiAarginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (16/10);  new base A (72/84);  total (88/96)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

ATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTCGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGT  >  CP009273/4329629‑4329916
                                                                                                                                                 |                                                                                                                                              
      cGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGa                                                                                                                                        <  1:810196/148‑1 (MQ=255)
          tAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAg                                                                                                                                    <  2:3188978/148‑1 (MQ=255)
          tAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAg                                                                                                                                    <  1:1171661/148‑1 (MQ=255)
           agagCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGa                                                                                                                                   <  1:1252444/148‑1 (MQ=255)
            gagCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGAt                                                                                                                                  <  1:2403544/148‑1 (MQ=255)
             agCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTCGGTGCGGAACAGAtt                                                                                                                                 >  2:783579/1‑148 (MQ=255)
             agCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTCGGTGCGGAACAGAtt                                                                                                                                 >  2:784022/1‑148 (MQ=255)
                cGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTc                                                                                                                              <  2:657482/148‑1 (MQ=255)
                cGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTc                                                                                                                              <  2:881132/148‑1 (MQ=255)
                 ggAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCa                                                                                                                             >  1:545134/1‑148 (MQ=255)
                 ggAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCa                                                                                                                             >  1:3201303/1‑148 (MQ=255)
                 ggAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCa                                                                                                                             <  2:823121/148‑1 (MQ=255)
                  gAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTCGGTGCGGAACAGATTTTCAc                                                                                                                            >  1:3142538/1‑148 (MQ=255)
                   aaGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCAcc                                                                                                                           <  2:2583183/148‑1 (MQ=255)
                   aaGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCAcc                                                                                                                           >  1:2648423/1‑148 (MQ=255)
                     gTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCAt                                                                                                                         >  2:108121/1‑148 (MQ=255)
                      tACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCAta                                                                                                                        >  1:1283814/1‑148 (MQ=255)
                        ccGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCAtagt                                                                                                                      <  1:493419/148‑1 (MQ=255)
                        ccGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCAtagt                                                                                                                      <  2:2722812/148‑1 (MQ=255)
                         cgacgaCTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCAtagta                                                                                                                     <  1:2989021/148‑1 (MQ=255)
                           acgacTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATagtagt                                                                                                                   <  1:1403483/148‑1 (MQ=255)
                            cgacTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTc                                                                                                                  <  1:2283747/148‑1 (MQ=255)
                              acTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCAt                                                                                                                >  1:517041/1‑148 (MQ=255)
                                tACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATgg                                                                                                              >  1:706911/1‑148 (MQ=255)
                                 aCCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGt                                                                                                             >  1:1248173/1‑148 (MQ=255)
                                  ccGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTa                                                                                                            <  2:196809/148‑1 (MQ=255)
                                      ccAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGaaa                                                                                                        <  2:2076417/148‑1 (MQ=255)
                                       cAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAAc                                                                                                       <  1:2737888/148‑1 (MQ=255)
                                          gAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTc                                                                                                    >  1:1684041/1‑148 (MQ=255)
                                          gAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGAAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTc                                                                                                    <  1:1233216/148‑1 (MQ=255)
                                          gAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTAAAGAAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTc                                                                                                    <  1:1233516/148‑1 (MQ=255)
                                           aGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTcc                                                                                                   <  2:462926/148‑1 (MQ=255)
                                                  ggCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGgt                                                                                            >  1:1769407/1‑148 (MQ=255)
                                                      ccAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGttt                                                                                        <  2:1529758/148‑1 (MQ=255)
                                                       cAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCCTGGTAGAAACGTCCGGCGGGTGtttt                                                                                       >  2:2637755/1‑148 (MQ=255)
                                                       cAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGtttt                                                                                       >  1:2893327/1‑148 (MQ=255)
                                                       cAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGtttt                                                                                       >  2:1733273/1‑148 (MQ=255)
                                                       cAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCAAGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGtttt                                                                                       >  1:2832478/1‑148 (MQ=255)
                                                         aaTACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTgg                                                                                     >  1:189708/1‑148 (MQ=255)
                                                         aaTACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTgg                                                                                     <  1:2474608/148‑1 (MQ=255)
                                                         aaTACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTgg                                                                                     >  1:1498557/1‑148 (MQ=255)
                                                          aTACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGt                                                                                    >  2:2827741/1‑148 (MQ=255)
                                                           tACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTCGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTa                                                                                   >  2:331134/1‑148 (MQ=255)
                                                            aCGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTaa                                                                                  >  2:1726349/1‑148 (MQ=255)
                                                               cgTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAAc                                                                               >  2:2378288/1‑148 (MQ=255)
                                                                 tGCGGCATATTTTTCGCTTTCTCCAAATGCGCCATTATGGTCAAGCAATGATCCGAGGGATTTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTAATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAATCCa                                                                             <  2:961698/148‑1 (MQ=255)
                                                                 tGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCa                                                                             >  1:2694814/1‑148 (MQ=255)
                                                                     gCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAAcgc                                                                         >  2:1910599/1‑148 (MQ=255)
                                                                      cATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAAcgcc                                                                        >  1:663840/1‑148 (MQ=255)
                                                                         atTTTTCGCTTTCGCCAAATGCGCCATTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCAATGTCGGTGCGTAACATATTTTCACAATAGTAGTCATGGTATAAACGTCCGGCGGGTGTTTTGGTAAAACCAACgccgcc                                                                     <  1:212316/148‑1 (MQ=255)
                                                                         atTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACgccgcc                                                                     <  2:536164/148‑1 (MQ=255)
                                                                         atTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACgccgcc                                                                     >  1:3095765/1‑148 (MQ=255)
                                                                          tttttCGCTTTCGCCAAATGCGCCAGTATGATCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGccc                                                                    <  2:909351/148‑1 (MQ=255)
                                                                           ttttCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCt                                                                   >  2:1346876/1‑148 (MQ=255)
                                                                            tttcgctttcgcCAATTGCGCCAGTTGGGTCAAGAAAATAACCGTGGGAAGTTGGTTCGATGCCCTTGTGGGTGCGGAACAGTTTTTCACCATAGTTTTCATGTTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCtg                                                                  <  1:433038/148‑1 (MQ=255)
                                                                               cgctttcgcCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCtggtg                                                               >  2:1455595/1‑148 (MQ=255)
                                                                                 ctttcgcCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGc                                                             <  2:1613491/148‑1 (MQ=255)
                                                                                      gcCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTCGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGgc                                                        >  2:1204677/1‑148 (MQ=255)
                                                                                       ccAAATGCGCCATTATGGTCAAGCAAACAACCGCGGGAAGTTCGTTCGATGCCCTTGTAGGTGCGTAACATATTTTCACCATAGTATTCTTGTTATAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGgcg                                                       <  1:1499767/148‑1 (MQ=255)
                                                                                       ccAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGgcg                                                       <  2:1684041/148‑1 (MQ=255)
                                                                                       caaaaTGCGCAAGTATGGACAAGCAAAGAACCGTGGGATGTTCGTTCTTTGCCATTTTCGGTGGGGAACAGATTTTCACCATTGTTGTCTTGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGgcg                                                       <  1:699092/146‑1 (MQ=255)
                                                                                          aaTGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTg                                                    <  2:3191245/148‑1 (MQ=255)
                                                                                           aTGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGc                                                   >  1:2667977/1‑148 (MQ=255)
                                                                                                   gTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTCGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAat                                           <  2:3142538/148‑1 (MQ=255)
                                                                                                    tATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAata                                          >  1:2950243/1‑148 (MQ=255)
                                                                                                    tATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAata                                          >  2:2053765/1‑148 (MQ=255)
                                                                                                    tATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAata                                          >  1:2325967/1‑148 (MQ=255)
                                                                                                     aTGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAatat                                         <  2:849942/148‑1 (MQ=255)
                                                                                                       ggTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTc                                       <  2:252690/148‑1 (MQ=255)
                                                                                                       ggTCAACAAAAGACCCGTGGGACGGTCGTTCGATGCCCTTGTCGGTGCGGAACAGATTTACACATTATTAGTCATGGTATAAACGTCCGGCGGGTGTTTTGTTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTc                                       <  1:1432762/148‑1 (MQ=255)
                                                                                                         tCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCAt                                     >  2:3025533/1‑148 (MQ=255)
                                                                                                          cAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGACTATTCATg                                    >  2:2968697/1‑148 (MQ=255)
                                                                                                           aaGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATgg                                   <  2:2130773/148‑1 (MQ=255)
                                                                                                            aGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGa                                  >  2:2127371/1‑148 (MQ=255)
                                                                                                             gCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGACCAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGAt                                 <  2:1009724/148‑1 (MQ=255)
                                                                                                             gCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGAt                                 >  2:2030710/1‑148 (MQ=255)
                                                                                                              cAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGACTATTCATGGATg                                >  2:1308245/1‑148 (MQ=255)
                                                                                                              cAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATg                                <  2:3295753/148‑1 (MQ=255)
                                                                                                              cAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATg                                <  2:3295562/148‑1 (MQ=255)
                                                                                                              cAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATg                                <  2:164396/148‑1 (MQ=255)
                                                                                                               aaaGAACCGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGt                               <  2:2781470/148‑1 (MQ=255)
                                                                                                                  gAACCGAGGGAAGTTCGTTCGATGCCCATGTCGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCAc                            >  2:2616910/1‑148 (MQ=255)
                                                                                                                     ccGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACtat                         <  2:3090909/148‑1 (MQ=255)
                                                                                                                      cGAGGGAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACtata                        <  1:741804/148‑1 (MQ=255)
                                                                                                                         gggAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATAttt                     >  1:1977192/1‑148 (MQ=255)
                                                                                                                         gggAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATAttt                     >  1:101269/1‑148 (MQ=255)
                                                                                                                          ggAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTtc                    <  1:2935750/148‑1 (MQ=255)
                                                                                                                          ggAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTtc                    <  1:2930532/148‑1 (MQ=255)
                                                                                                                           gAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTtca                   >  1:794011/1‑148 (MQ=255)
                                                                                                                           gAAGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTtca                   >  1:2710314/1‑148 (MQ=255)
                                                                                                                            atgttcgttcGATGCCCATGTAGGTGCGGAACAGATTTTCACCATCGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTtcat                  <  1:1317677/146‑1 (MQ=255)
                                                                                                                            aaGTTCGTTCGATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTtcat                  <  1:1067677/148‑1 (MQ=255)
                                                                                                                              gttcgttcGATGCCCATGTCGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTtcatca                >  2:1896277/1‑148 (MQ=255)
                                                                                                                              gttcgttcGATGCCCATGTCGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTtcatca                >  2:2093252/1‑148 (MQ=255)
                                                                                                                               ttcgttcgATGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAg               >  1:425421/1‑148 (MQ=255)
                                                                                                                                        tGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGaa      <  2:706631/148‑1 (MQ=255)
                                                                                                                                        tGCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGaa      <  1:108121/148‑1 (MQ=255)
                                                                                                                                         gCCCATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAAc     >  2:2048171/1‑148 (MQ=255)
                                                                                                                                           ccATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAg   <  1:2928324/148‑1 (MQ=255)
                                                                                                                                           ccATGTAGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAg   <  1:1336260/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
ATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCATGTCGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGT  >  CP009273/4329629‑4329916

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: