Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,464,440 T→C pseudogene (147/291 nt) BW25113_RS12275 → tail fiber assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,464,4400TC93.0% 289.8 / 9.7 86pseudogene (147/291 nt)BW25113_RS12275tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base C (35/45);  total (39/47)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

ACGCCATACGATAAATGGGATGGTGAGAAATGGGTGACGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATTCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCG  >  NZ_CP009273/2464309‑2464575
                                                                                                                                   |                                                                                                                                       
aCGCCATACGATAAATGGGATGGTGAGAAATGGGTGACGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTg                                                                                                                                   <  1:201644/138‑1 (MQ=255)
      tACGATAAATGGGATGGTGAGAAATGGGTGACGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTa                                                                                                                             >  1:476344/1‑138 (MQ=255)
       aCGATAAATGGGATGGTGAGAAATGGGTGACGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTa                                                                                                                             >  1:154448/1‑137 (MQ=255)
          aTAAATGGGATGGTGAGAAATGGGTGACGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGc                                                                                                                        <  2:109239/139‑1 (MQ=255)
          aTAAATGGGATGGTGAGAAATGGGTGACGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGc                                                                                                                        <  2:1175148/139‑1 (MQ=255)
            aaaTGGGATGGTGAGAAATGGGTGACGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTg                                                                                                                      >  1:925561/1‑139 (MQ=255)
                 ggATGGTGAGAAATGGGTGACGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCg                                                                                                                  <  1:292423/138‑1 (MQ=255)
                   aTGGTGAGAAATGGGTGACGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCgg                                                                                                                 >  2:698737/1‑137 (MQ=255)
                   aTGGTGAGAAATGGGTGACGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCgg                                                                                                                 >  1:53258/1‑137 (MQ=255)
                     ggTGAGAAATGGGTGACGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGc                                                                                                             >  2:165721/1‑139 (MQ=255)
                      gTGAGAAATGGGTGACGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCt                                                                                                            <  2:948185/139‑1 (MQ=255)
                          gaAATGGGTGACGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATg                                                                                                        <  1:447246/139‑1 (MQ=255)
                               gggTGACGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATgcagg                                                                                                    >  1:213817/1‑138 (MQ=255)
                                  tGACGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCaga                                                                                                <  2:53258/139‑1 (MQ=255)
                                    aCGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATTCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAgaga                                                                                              >  1:1192322/1‑139 (MQ=255)
                                     cGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGAc                                                                                             >  1:209648/1‑139 (MQ=255)
                                         tACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTTGGCGGAAGCTGATGCAGGCAGAGACCTcc                                                                                         <  2:1290563/139‑1 (MQ=255)
                                         tACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTcc                                                                                         <  1:680995/139‑1 (MQ=255)
                                         tACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTcc                                                                                         <  1:557830/139‑1 (MQ=255)
                                         tACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTcc                                                                                         <  2:925561/139‑1 (MQ=255)
                                         tACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTcc                                                                                         <  2:341716/139‑1 (MQ=255)
                                         tACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTcc                                                                                         <  2:745858/139‑1 (MQ=255)
                                         tACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTcc                                                                                         <  1:221970/139‑1 (MQ=255)
                                          aCCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTccc                                                                                        >  1:414599/1‑139 (MQ=255)
                                              aGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAgg                                                                                                                         <  2:1126644/102‑1 (MQ=255)
                                              aGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAgg                                                                                                                         >  1:1126644/1‑102 (MQ=255)
                                                gCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTcc                                                                                         <  1:407169/132‑1 (MQ=255)
                                                gCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTcc                                                                                         >  2:407169/1‑132 (MQ=255)
                                                gCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACtt                                                                                   >  1:1211342/1‑138 (MQ=255)
                                                   caGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATTCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAaca                                                                               <  2:770892/139‑1 (MQ=255)
                                                   caGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAaca                                                                               <  2:476344/139‑1 (MQ=255)
                                                   caGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAaca                                                                               <  2:328787/139‑1 (MQ=255)
                                                   caGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAaca                                                                               <  2:279929/139‑1 (MQ=255)
                                                   caGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAaca                                                                               <  1:698737/139‑1 (MQ=255)
                                                    aGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAacac                                                                              >  2:667925/1‑139 (MQ=255)
                                                        tAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACtgtg                                                                          <  1:1101509/139‑1 (MQ=255)
                                                        tAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACtgtg                                                                          <  1:805912/139‑1 (MQ=255)
                                                          gCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGc                                                                         <  1:165721/138‑1 (MQ=255)
                                                           cGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTg                                                                       >  2:440802/1‑139 (MQ=255)
                                                             tCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGa                                                                     <  2:154604/139‑1 (MQ=255)
                                                              cGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTcc                                                                                         >  2:645991/1‑118 (MQ=255)
                                                              cGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTcc                                                                                         <  1:645991/118‑1 (MQ=255)
                                                                cAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTcc                                                                                         >  2:72163/1‑116 (MQ=255)
                                                                cAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTcc                                                                                         <  1:72163/116‑1 (MQ=255)
                                                                cAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACtg                                                                            >  1:8095/1‑129 (MQ=255)
                                                                cAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACtg                                                                            <  2:8095/129‑1 (MQ=255)
                                                                           gcagAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATTCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCgg                                                       >  2:764370/1‑139 (MQ=255)
                                                                           gcagAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCgg                                                       >  1:4925/1‑139 (MQ=255)
                                                                            cagAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGt                                                      <  2:159012/139‑1 (MQ=255)
                                                                                 gCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTcc                                                                                         >  1:1163649/1‑99 (MQ=255)
                                                                                 gCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTcc                                                                                         <  2:1163649/99‑1 (MQ=255)
                                                                                 gCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACgg                                                 >  1:836742/1‑139 (MQ=255)
                                                                                      gcgcCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCaa                                              >  1:1227680/1‑137 (MQ=255)
                                                                                      gcgcCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACa                                            <  1:201513/139‑1 (MQ=255)
                                                                                      gcgcCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACa                                            <  1:667925/139‑1 (MQ=255)
                                                                                      gcgcCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACa                                            <  2:139526/139‑1 (MQ=255)
                                                                                      gcgcCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACa                                            <  2:829127/139‑1 (MQ=255)
                                                                                            gTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATAc                                       >  1:858911/1‑138 (MQ=255)
                                                                                             tCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCACGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATAc                                       >  1:921631/1‑137 (MQ=255)
                                                                                                       atACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTaa                                                                                 <  1:796726/85‑1 (MQ=255)
                                                                                                       atACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTaa                                                                                 >  2:796726/1‑85 (MQ=255)
                                                                                                       atACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCgcgc                            <  2:836742/138‑1 (MQ=255)
                                                                                                       atACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCgcgc                            <  2:858911/138‑1 (MQ=255)
                                                                                                       atACTGCAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGcc                           >  2:660389/1‑139 (MQ=255)
                                                                                                             cAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCagag                                                                                               <  2:695415/65‑1 (MQ=255)
                                                                                                             cAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCagag                                                                                               >  1:695415/1‑65 (MQ=255)
                                                                                                             cAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCaga                                                                                                >  1:376726/1‑64 (MQ=255)
                                                                                                             cAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCaga                                                                                                <  2:376726/64‑1 (MQ=255)
                                                                                                             cAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGt                     >  2:541377/1‑139 (MQ=255)
                                                                                                             cAATGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGt                     >  2:822940/1‑139 (MQ=255)
                                                                                                               aTGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAaca                                                                               <  1:272765/79‑1 (MQ=255)
                                                                                                               aTGGCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAaca                                                                               >  2:272765/1‑79 (MQ=255)
                                                                                                                  gCTTCCATTAGTCTGATTCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACt                                                                                    >  2:990076/1‑71 (MQ=255)
                                                                                                                  gCTTCCATTAGTCTGATTCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACt                                                                                    <  1:990076/71‑1 (MQ=255)
                                                                                                                  gCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTaa                                                                                 <  2:446507/74‑1 (MQ=255)
                                                                                                                  gCTTCCATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTaa                                                                                 >  1:446507/1‑74 (MQ=255)
                                                                                                                       cATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACgg                                                 >  1:848407/1‑101 (MQ=255)
                                                                                                                       cATTAGTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACgg                                                 <  2:848407/101‑1 (MQ=255)
                                                                                                                            gTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTaa                                                                                 <  2:76597/64‑1 (MQ=255)
                                                                                                                            gTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTaa                                                                                 >  1:76597/1‑64 (MQ=255)
                                                                                                                            gTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAaca                                                                               <  2:553294/66‑1 (MQ=255)
                                                                                                                            gTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAaca                                                                               >  1:553294/1‑66 (MQ=255)
                                                                                                                            gTCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAAc       <  2:742410/138‑1 (MQ=255)
                                                                                                                             tCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGc                                                >  2:1101914/1‑96 (MQ=255)
                                                                                                                             tCTGATCCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGc                                                <  1:1101914/96‑1 (MQ=255)
                                                                                                                                 aTTCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCg  >  1:992066/1‑138 (MQ=255)
                                                                                                                                   |                                                                                                                                       
ACGCCATACGATAAATGGGATGGTGAGAAATGGGTGACGGATACCGAGGCACAGCATAGCGTCGCAGTAGATGCAGCAGAAGCACAGCGCCAGTCGCTGATTGATACTGCAATGGCTTCCATTAGTCTGATTCAACTGAAATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCG  >  NZ_CP009273/2464309‑2464575

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