Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273351,8220TC63.9% 74.9 / 94.5 83I397T (ATC→ACC) codBcytosine permease
Reads supporting (aligned to +/- strand):  ref base T (16/14);  new base C (28/25);  total (44/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTACCAGGCGAAGAGGGGCTGTGGCAGATTCATCTGCAG  >  NZ_CP009273/351690‑351959
                                                                                                                                    |                                                                                                                                         
ttGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTaa                                                                                                                                     >  2:453681/1‑139 (MQ=255)
         cgcgTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGAtttt                                                                                                                            >  1:567484/1‑139 (MQ=255)
         cgcgTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGAtttt                                                                                                                            >  1:531004/1‑139 (MQ=255)
         cgcgTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGGTACAGGGAATTGTTCCGGTCAAGGCGGTATTAGGTGGCGCGCTGAGGTATCTGATCCTTAAgcctgtctc                                                                                                                            >  2:235632/1‑133 (MQ=255)
             tATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAAt                                                                                                                        >  1:525790/1‑139 (MQ=255)
                    aGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTaaaa                                                                                                                 >  1:385282/1‑139 (MQ=255)
                    aGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTaaaa                                                                                                                 >  2:97758/1‑139 (MQ=255)
                    aGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTaaaa                                                                                                                 >  2:556988/1‑139 (MQ=255)
                    aGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTaaaa                                                                                                                 >  1:137451/1‑139 (MQ=255)
                        tCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGa                                                                                                              >  1:114620/1‑138 (MQ=255)
                          aaTTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGAcag                                                                                                           <  1:556988/139‑1 (MQ=255)
                           aTTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGAcag                                                                                                           >  2:350150/1‑138 (MQ=255)
                           aTTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGAcag                                                                                                           >  1:121722/1‑138 (MQ=255)
                           aTTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGAcag                                                                                                           >  2:251510/1‑138 (MQ=255)
                               ggTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAgcagca                                                                                                      <  2:567484/139‑1 (MQ=255)
                               ggTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAgcagca                                                                                                      <  2:121722/139‑1 (MQ=255)
                                gTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCaa                                                                                                     >  1:17765/1‑139 (MQ=255)
                                 tGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAAt                                                                                                    <  1:350150/139‑1 (MQ=255)
                                   gCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGa                                                                                                  <  1:534724/139‑1 (MQ=255)
                                    cGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGAc                                                                                                 <  1:251510/139‑1 (MQ=255)
                                     gATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACg                                                                                                <  1:159142/139‑1 (MQ=255)
                                            gCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCAtgtg                                                                                          <  1:302761/138‑1 (MQ=255)
                                            gCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCAtgtg                                                                                          <  1:449156/138‑1 (MQ=255)
                                            gCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCAtgtg                                                                                          <  2:17765/138‑1 (MQ=255)
                                            gCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCAtgtg                                                                                          <  2:144883/138‑1 (MQ=255)
                                            gCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCAtgtg                                                                                          <  2:137451/138‑1 (MQ=255)
                                            gCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCAtgtg                                                                                          <  2:66425/138‑1 (MQ=255)
                                                   ccTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTaa                                                                                  <  2:130020/139‑1 (MQ=255)
                                                   ccTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTaa                                                                                  <  1:244146/139‑1 (MQ=255)
                                                        gggATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAgtg                                                                             <  2:531004/139‑1 (MQ=255)
                                                        gggATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAGCCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAgtg                                                                             <  1:393077/139‑1 (MQ=255)
                                                        gggATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAgtg                                                                             <  2:185979/139‑1 (MQ=255)
                                                             tgctgcAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGaa                                                                        >  2:436089/1‑139 (MQ=255)
                                                              gctgcAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACa                                                                                <  1:222067/130‑1 (MQ=255)
                                                              gctgcAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACa                                                                                >  2:222067/1‑130 (MQ=255)
                                                              gctgcAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTaa                                                                                                                   >  1:558408/1‑95 (MQ=255)
                                                              gctgcAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTaa                                                                                                                   <  2:558408/95‑1 (MQ=255)
                                                                         cTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACagcagc                                                                                                       <  2:310721/96‑1 (MQ=255)
                                                                         cTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACagcagc                                                                                                       >  1:310721/1‑96 (MQ=255)
                                                                         cTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTaca                                                            <  1:398255/139‑1 (MQ=255)
                                                                         cTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTac                                                             >  2:492700/1‑138 (MQ=255)
                                                                          tGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTaca                                                            >  2:357583/1‑138 (MQ=255)
                                                                          tGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTaca                                                            >  1:293404/1‑138 (MQ=255)
                                                                          tGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTaca                                                            >  2:578160/1‑138 (MQ=255)
                                                                            gTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTacaaa                                                          <  2:45036/138‑1 (MQ=255)
                                                                                 cGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCAtgtg                                                                                          <  1:177673/101‑1 (MQ=255)
                                                                                 cGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCAtgtg                                                                                          >  2:177673/1‑101 (MQ=255)
                                                                                 cGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTaca                                                            >  2:395528/1‑131 (MQ=255)
                                                                                 cGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTaca                                                            <  1:395528/131‑1 (MQ=255)
                                                                                 cGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAatta                                                    >  2:275268/1‑139 (MQ=255)
                                                                                    gAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCAtgtg                                                                                          <  1:24811/98‑1 (MQ=255)
                                                                                    gAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCAtgtg                                                                                          >  2:24811/1‑98 (MQ=255)
                                                                                     aaTTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTaa                                                <  1:579070/139‑1 (MQ=255)
                                                                                      aTTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAAc                                               <  1:32366/139‑1 (MQ=255)
                                                                                          ttCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGccc                                           >  1:279476/1‑139 (MQ=255)
                                                                                               gTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTa                                      >  2:238533/1‑139 (MQ=255)
                                                                                                  aaCGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTAc                                     >  2:122851/1‑137 (MQ=255)
                                                                                                  aaCGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTAc                                     >  1:239344/1‑137 (MQ=255)
                                                                                                    cgcgGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTACCAgg                                 <  2:486609/139‑1 (MQ=255)
                                                                                                    cgcgGTATTAGGTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTACCAgg                                 <  1:275268/139‑1 (MQ=255)
                                                                                                               gTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCAtgtg                                                                                          <  2:299722/71‑1 (MQ=255)
                                                                                                               gTGGCGCGCTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCAtgtg                                                                                          >  1:299722/1‑71 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACa                                                                                <  1:86789/77‑1 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACa                                                                                >  2:86789/1‑77 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTAc                                     <  1:100664/120‑1 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTAc                                     >  2:100664/1‑120 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTACCAGGCGAAGAGGGGCtgtg                  >  1:319915/1‑139 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTACCAGGCGAAGAGGGGCtgtg                  >  2:102331/1‑139 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTACCAGGCGAAGAGGGGCtgtg                  >  1:546082/1‑139 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGAc                                                                                                 >  1:423771/1‑60 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGAc                                                                                                 <  2:423771/60‑1 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCAtgtg                                                                                          <  2:581970/67‑1 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCAtgtg                                                                                          >  1:581970/1‑67 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGcc                                            >  2:414581/1‑113 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGcc                                            <  1:414581/113‑1 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTAc                                     >  2:277284/1‑120 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTAc                                     <  1:277284/120‑1 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTACCAgg                                 <  1:59937/124‑1 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTACCAgg                                 >  2:59937/1‑124 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTACCAGGCGAAGAGGGGCtgtg                  >  2:156317/1‑139 (MQ=255)
                                                                                                                   cgcgcTGAGCTATCTGACCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTACCAGGCGAAGAGGGGCtgtg                  >  1:453753/1‑139 (MQ=255)
                                                                                                                     cgcTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTACCAGGCGAAGAGGGGCTGTgg                 >  2:117630/1‑138 (MQ=255)
                                                                                                                      gcTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTACCAGGCGAAGAGGGGCTGTGGCa               <  1:63305/139‑1 (MQ=255)
                                                                                                                                    tCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTACCAGGCGAAGAGGGGCTGTGGCAGATTCATCTGCAg  >  2:265520/1‑138 (MQ=255)
                                                                                                                                    |                                                                                                                                         
TTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGAGCTATCTGATCCTTAACCCGATTTTGAATCGTAAAACGACAGCAGCAATGACGCATGTGGAGGCTAACAGTGTCGAATAACGCTTTACAAACAATTATTAACGCCCGGTTACCAGGCGAAGAGGGGCTGTGGCAGATTCATCTGCAG  >  NZ_CP009273/351690‑351959

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