Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,141,321 (T)5→4 coding (1706/2652 nt) ppc ← phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,141,3170T.90.3% 253.0 / 24.8 72coding (1710/2652 nt)ppcphosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base . (32/33);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-02

CCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATAC  >  NZ_CP009273/4141187‑4141449
                                                                                                                                   |                                                                                                                                    
ccTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAg                                                                                                                             >  1:465629/1‑139 (MQ=255)
    ttCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGttt                                                                                                                         >  2:1231739/1‑139 (MQ=255)
    ttCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGgtt                                                                                                                         >  1:1100406/1‑139 (MQ=255)
        gTTCCGCGCGGGCCGCTTTTCGGTCTTCCGGGCGGTTGTGACAGCAGCGCCGCATGAGGAGGTGCGCCGCCGCAACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGa                                                                                                                      <  2:862809/139‑1 (MQ=255)
             gcgcAGGCCGCTTTTCAGGCTTCCTGGGGGTTGGGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGCCCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATg                                                                                                                <  1:316016/139‑1 (MQ=255)
             gcgcAGGCCGCCTTTCGGGCTTCCTGGGGGTTGGGACAGCAGCGCCGCATGAGCGGGTGCCCCGCCGCGCCCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATg                                                                                                                <  2:465629/139‑1 (MQ=255)
              cgcAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATg                                                                                                                >  1:1093238/1‑139 (MQ=255)
              cgcAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGAttctgtc                                                                                                               >  1:283805/1‑134 (MQ=255)
              cgcAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATg                                                                                                                >  2:321103/1‑138 (MQ=255)
              cgcAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATg                                                                                                                >  2:175120/1‑138 (MQ=255)
              cgcAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGc                                                                                                               >  2:98563/1‑139 (MQ=255)
                 gggCCGCCTTTCGGGCTTCCTGGGGGTTGTGCCAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATc                                                                                                            <  2:358501/138‑1 (MQ=255)
                   gccgccTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCAGGAGCAGGTGCGCCGCCGCGACCAAGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCt                                                                                                          <  1:98563/139‑1 (MQ=255)
                   gccgccTTCCAGGCTTCCTGGGGGTTGGGACAGCAGCGCCGAAGGAGCGGGTGCGCCGCCGCGACCAATGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCt                                                                                                          <  2:124398/139‑1 (MQ=255)
                   gcccccTGTCAGGCTTCCTGGCGTTTGTGCCGGCAGCGCCGCATGAGCAGGGGCGCCGCCGCGCCCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCt                                                                                                          <  2:1278162/139‑1 (MQ=255)
                    ccgccTTTCAGGCTCCCGGGGGTTTGTGACACCCGCCCCGACTGGGGGGGTGCGCCGCCGCTACAAAGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCt                                                                                                          <  1:640400/139‑1 (MQ=255)
                       ccTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGc                                                                                                      >  1:713788/1‑139 (MQ=255)
                                  tCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATTGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGACATTacg                                                                                           >  1:1155320/1‑139 (MQ=255)
                                    cTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGcc                                                                                         >  2:1170248/1‑139 (MQ=255)
                                       ggggTTGTGACAGCGCCGCCGCATGAGGGGGGGCCCCCCCGCGACCAAGGTAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAg                                                                                      <  2:1220807/137‑1 (MQ=255)
                                        cGGTGGTGACCGCAGCGCGGCAGGGGGGGGGGCGCCCCCGCCACCAGGGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAgg                                                                                     <  2:656844/139‑1 (MQ=255)
                                              tgACAGCGGCGCCGCATGAGCAGGTGCGCCGCCGCAACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCt                                                                                <  1:904134/138‑1 (MQ=255)
                                                   ccggcgcCGCGTGGGCGGGGGCGCCCCCGCGACAATGGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATc                                                                          <  2:867842/136‑1 (MQ=255)
                                                       cgccgcATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACt                                                                       >  1:1033060/1‑139 (MQ=255)
                                                       cgccgcATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTc                                                                      >  2:640122/1‑139 (MQ=255)
                                                          cgctggAGCGGGTGCGCCCCCGCCACCAAGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTccc                                                                    <  2:1093238/134‑1 (MQ=255)
                                                           gcATGAGCAGGCGCCCCGCCCCGCCAAATGGACCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTCCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGc                                                                  <  2:111042/139‑1 (MQ=255)
                                                                    ggggcgcCGCCGCGACAAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGc                                                         <  1:71617/136‑1 (MQ=255)
                                                                    ggTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTCCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGc                                                         >  1:111042/1‑139 (MQ=255)
                                                                     gTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCt                                                        >  1:896953/1‑139 (MQ=255)
                                                                        cgcccccgcGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAg                                                     <  2:22878/139‑1 (MQ=255)
                                                                        ccccgccgaGCCCAAGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAg                                                     <  2:1155320/137‑1 (MQ=255)
                                                                         gccgccgcGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGt                                                    >  2:298471/1‑139 (MQ=255)
                                                                         gccgccgcGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGt                                                    <  2:792863/139‑1 (MQ=255)
                                                                         gccgccgcGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGt                                                    >  2:1188824/1‑139 (MQ=255)
                                                                          ccgccgcGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTc                                                   >  1:529115/1‑139 (MQ=255)
                                                                           cgccgcGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCg                                                  >  1:990422/1‑139 (MQ=255)
                                                                             ccgcGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGa                                                <  2:529115/139‑1 (MQ=255)
                                                                               gcgACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGaa                                               >  2:1083794/1‑138 (MQ=255)
                                                                                 gcccAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAg                                            <  1:1144042/137‑1 (MQ=255)
                                                                                        ggAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAatca                                     >  2:736579/1‑139 (MQ=255)
                                                                                          aaaCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAatcat                                    <  1:955793/138‑1 (MQ=255)
                                                                                            aCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAAtcatca                                  >  2:66065/1‑139 (MQ=255)
                                                                                              cgccgcGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTg                                                          <  2:390362/113‑1 (MQ=255)
                                                                                              cgccgcGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTg                                                          >  1:390362/1‑113 (MQ=255)
                                                                                                ccgcGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTccc                                                                    >  2:355685/1‑101 (MQ=255)
                                                                                                ccgcGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTccc                                                                    <  1:355685/101‑1 (MQ=255)
                                                                                                ccgcGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCAt                              >  1:534924/1‑139 (MQ=255)
                                                                                                ccgcGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCAt                              >  1:828669/1‑139 (MQ=255)
                                                                                                    gACCGTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGt                          <  1:175120/139‑1 (MQ=255)
                                                                                                        gTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGCTATTGCGCCCAGGAAGCTGCCATCACTCCCg                                                                   >  2:1137318/1‑94 (MQ=255)
                                                                                                        gTGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCg                                                                   <  1:1137318/94‑1 (MQ=255)
                                                                                                         tGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGa                                                                                                                      >  2:187787/1‑42 (MQ=255)
                                                                                                         tGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGa                                                                                                                      <  1:187787/42‑1 (MQ=255)
                                                                                                         tGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGcc                     >  1:702850/1‑139 (MQ=255)
                                                                                                         tGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAACGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGtt                         <  1:668684/135‑1 (MQ=255)
                                                                                                         tGGAACAACGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAACGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGtt                         >  2:668684/1‑135 (MQ=255)
                                                                                                               aacGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAAt               <  1:904209/139‑1 (MQ=255)
                                                                                                                acGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCt                                                                                                          <  2:768460/47‑1 (MQ=255)
                                                                                                                acGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCt                                                                                                          >  1:768460/1‑47 (MQ=255)
                                                                                                                 cGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCt                                                                                                          >  1:259008/1‑46 (MQ=255)
                                                                                                                 cGTCAGCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCt                                                                                                          <  2:259008/46‑1 (MQ=255)
                                                                                                                      gCTCAATACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCa        <  1:1168980/139‑1 (MQ=255)
                                                                                                                       cTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCa        <  1:930999/139‑1 (MQ=255)
                                                                                                                          aaTACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGa     <  2:896953/138‑1 (MQ=255)
                                                                                                                           aTACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa   <  2:534924/139‑1 (MQ=255)
                                                                                                                           aTACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa   <  2:747979/139‑1 (MQ=255)
                                                                                                                           aTACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa   <  1:736579/139‑1 (MQ=255)
                                                                                                                           aTACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa   >  1:747979/1‑139 (MQ=255)
                                                                                                                           aTACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa   <  1:1098056/139‑1 (MQ=255)
                                                                                                                            tACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa   <  1:923646/139‑1 (MQ=255)
                                                                                                                            tACCCGCTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATAc  >  1:162258/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                    
CCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATAC  >  NZ_CP009273/4141187‑4141449

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