Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,994,324 A→G V31V (GTT→GTC yigE ← phosphodiester glycosidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,994,3240AG88.6% 217.7 / 20.2 79V31V (GTT→GTCyigEphosphodiester glycosidase family protein
Reads supporting (aligned to +/- strand):  ref base A (4/5);  new base G (40/30);  total (44/35)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

AGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACC  >  NZ_CP009273/3994190‑3994458
                                                                                                                                      |                                                                                                                                      
agcagAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAg                                                                                                                                    <  2:563278/139‑1 (MQ=255)
     agCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCaaa                                                                                                                                >  1:519806/1‑138 (MQ=255)
       cATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAAc                                                                                                                               >  1:168999/1‑137 (MQ=255)
       cATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAAc                                                                                                                               >  2:689697/1‑137 (MQ=255)
       cATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACaa                                                                                                                             <  2:500576/139‑1 (MQ=255)
         tGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAgg                                                                                                                           <  1:734287/139‑1 (MQ=255)
                 tCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAg                                                                                                                   >  2:66382/1‑139 (MQ=255)
                   cccACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGt                                                                                                                 >  1:837362/1‑139 (MQ=255)
                    ccACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTg                                                                                                                <  1:238519/139‑1 (MQ=255)
                    ccACGCTTCACCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTg                                                                                                                <  2:822655/139‑1 (MQ=255)
                     cACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGt                                                                                                                 >  2:403261/1‑137 (MQ=255)
                        gCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAggg                                                                                                            <  1:555364/139‑1 (MQ=255)
                                  tGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGCGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGaaaaac                                                                                                  >  2:257123/1‑138 (MQ=255)
                                         ttGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGttt                                                                                            <  1:1210979/138‑1 (MQ=255)
                                              aGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGAtt                                                                                      >  1:237537/1‑139 (MQ=255)
                                               gTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTc                                                                                     <  2:804233/139‑1 (MQ=255)
                                                  cATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAg                                                                                  >  1:329082/1‑139 (MQ=255)
                                                   aTTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAgg                                                                                 <  2:329082/139‑1 (MQ=255)
                                                     tttCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACaa                                                                                                                             >  1:87150/1‑93 (MQ=255)
                                                     tttCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACaa                                                                                                                             <  2:87150/93‑1 (MQ=255)
                                                     tttCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGt                                                                                >  1:1261633/1‑138 (MQ=255)
                                                     tttCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGt                                                                                >  2:1186290/1‑138 (MQ=255)
                                                        cACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGAt                                                                             <  2:519806/138‑1 (MQ=255)
                                                                 tGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCtt                                                                     <  1:638695/137‑1 (MQ=255)
                                                                          aTTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGt                                                          <  2:837362/139‑1 (MQ=255)
                                                                             aaCGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAAt                                                              >  2:561007/1‑132 (MQ=255)
                                                                             aaCGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAAt                                                              <  1:561007/132‑1 (MQ=255)
                                                                             aaCGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGc                                                       >  1:1167614/1‑139 (MQ=255)
                                                                             aaCGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGCGATCATTCCTTTACCAATGAGTAGc                                                       >  2:841372/1‑139 (MQ=255)
                                                                                  tataCGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATg                                                  <  1:877552/139‑1 (MQ=255)
                                                                                  tataCGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATg                                                  <  2:2996/139‑1 (MQ=255)
                                                                                            tACGGTCAGCGTCGTATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATta                                          >  1:1157700/1‑137 (MQ=255)
                                                                                            tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTGCCAATGAGTAGCTGATGCGCCATta                                          >  2:1108074/1‑137 (MQ=255)
                                                                                            tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGGTGCGCCATta                                          >  2:235529/1‑137 (MQ=255)
                                                                                            tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATgcg                                                >  2:447672/1‑131 (MQ=255)
                                                                                            tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATgcg                                                <  1:447672/131‑1 (MQ=255)
                                                                                            tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATgc                                                 >  1:748298/1‑130 (MQ=255)
                                                                                            tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATgc                                                 <  2:748298/130‑1 (MQ=255)
                                                                                            tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATtaca                                        <  1:219650/139‑3 (MQ=255)
                                                                                            tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATtaca                                        <  2:794860/139‑3 (MQ=255)
                                                                                            tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATta                                          >  2:219650/1‑137 (MQ=255)
                                                                                            tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATta                                          >  2:19679/1‑137 (MQ=255)
                                                                                            tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATta                                          >  2:104852/1‑137 (MQ=255)
                                                                                            tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATta                                          >  1:1231202/1‑137 (MQ=255)
                                                                                            tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATta                                          >  1:180808/1‑137 (MQ=255)
                                                                                            tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATta                                          >  1:794860/1‑137 (MQ=255)
                                                                                            tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATta                                          >  1:851572/1‑137 (MQ=255)
                                                                                            tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATta                                          >  1:1273148/1‑137 (MQ=255)
                                                                                                gTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTc                                    <  2:1167614/139‑1 (MQ=255)
                                                                                                 tCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCa                                                                                    <  2:1183889/90‑1 (MQ=255)
                                                                                                 tCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCa                                                                                    >  1:1183889/1‑90 (MQ=255)
                                                                                                 tCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTcc                                   >  1:1206863/1‑139 (MQ=255)
                                                                                                   aGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTg                                 >  1:876527/1‑139 (MQ=255)
                                                                                                           aTCTGAGAGTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTGGATGTgg                          >  2:440586/1‑138 (MQ=255)
                                                                                                                agagTGCGCAATCATCAGCGTCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGATAAAAACCCGTTTGATATTCAAGGTGATCATTCCTTTACCAATGAGTACCTGATGCGCCATTACCGCTCCTGGATGTTGGa                        <  2:599080/135‑1 (MQ=255)
                                                                                                                agagTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTGGATGTGGGa                        <  2:876527/135‑1 (MQ=255)
                                                                                                                agagTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTGGATGTGGGAtttt                    >  1:132911/1‑139 (MQ=255)
                                                                                                                agagTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTGGATGTGGGAtt                      <  1:66382/137‑1 (MQ=255)
                                                                                                                agagTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGCGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTGGATGTGGGa                        <  1:841372/135‑1 (MQ=255)
                                                                                                                  agTGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTGGATGTGGGAtttttt                  <  2:132911/139‑1 (MQ=255)
                                                                                                                    tGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCa                                                                >  1:1103701/1‑91 (MQ=255)
                                                                                                                    tGCGCAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCa                                                                <  2:1103701/91‑1 (MQ=255)
                                                                                                                     gcgcAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTGGATGTGGGATTTTTTTAt               >  1:432393/1‑139 (MQ=255)
                                                                                                                     gcgcAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTGGATGTGGGATTTTTTTAt               >  2:325374/1‑139 (MQ=255)
                                                                                                                      cgcAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCtt                                                                     <  2:1184737/84‑1 (MQ=255)
                                                                                                                      cgcAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCtt                                                                     >  1:1184737/1‑84 (MQ=255)
                                                                                                                      cgcAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTGGATGTGGGATTTTTTTATc              >  1:692396/1‑139 (MQ=255)
                                                                                                                      cgcAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTGGATGTGGGATTTTTTTATc              >  2:1094046/1‑139 (MQ=255)
                                                                                                                      cgcAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTGGATGTGGGATTTTTTTATc              >  1:1260187/1‑139 (MQ=255)
                                                                                                                      cgcAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCt                                                                      <  1:412971/83‑1 (MQ=255)
                                                                                                                      cgcAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCt                                                                      >  2:412971/1‑83 (MQ=255)
                                                                                                                       gcAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTGGATGTGGGATTTTTTTATcc             <  2:598235/139‑1 (MQ=255)
                                                                                                                       gcAATCATCAGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTGGATGTGGGATTTTTTTATcc             <  1:1241307/139‑1 (MQ=255)
                                                                                                                             atcaGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAg       <  2:1157700/139‑1 (MQ=255)
                                                                                                                                aGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATtac                                         >  1:158545/1‑101 (MQ=255)
                                                                                                                                aGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATtac                                         <  2:158545/102‑2 (MQ=255)
                                                                                                                                aGCGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGa    >  2:1195380/1‑139 (MQ=255)
                                                                                                                                  cGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGAcc  >  1:378214/1‑139 (MQ=255)
                                                                                                                                  cGGCGACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAACCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTACAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGAcc  <  2:432393/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
AGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACC  >  NZ_CP009273/3994190‑3994458

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