Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,140,022 T→C I58V (ATC→GTC)  eptC ← phosphoethanolamine transferase CptA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,140,0220TC92.3% 277.3 / 9.8 91I58V (ATC→GTC) eptCphosphoethanolamine transferase CptA
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base C (37/47);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.47e-02

AGAACGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGC  >  NZ_CP009273/4139892‑4140156
                                                                                                                                  |                                                                                                                                      
aGAACGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACcaca                                                                                                                                <  2:962049/139‑1 (MQ=255)
aGAACGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACcaca                                                                                                                                <  2:10091/139‑1 (MQ=255)
 gAACGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAg                                                                                                                               >  1:852670/1‑139 (MQ=255)
  aaCGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGc                                                                                                                              >  1:784565/1‑139 (MQ=255)
   aCGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCg                                                                                                                             <  1:380892/139‑1 (MQ=255)
     gCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGaa                                                                                                                           <  2:62148/139‑1 (MQ=255)
     gCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGaa                                                                                                                           <  2:617729/139‑1 (MQ=255)
     gCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGaa                                                                                                                           <  2:500884/139‑1 (MQ=255)
     gCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGaa                                                                                                                           <  2:243484/139‑1 (MQ=255)
       tctGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGaa                                                                                                                           <  2:860709/137‑1 (MQ=255)
       tctGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGaa                                                                                                                           >  1:860709/1‑137 (MQ=255)
          gCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGaa                                                                                                                           >  1:1004412/1‑134 (MQ=255)
          gCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGaa                                                                                                                           <  2:1004412/134‑1 (MQ=255)
                 cTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGc                                                                                                               >  2:314683/1‑139 (MQ=255)
                     tGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGt                                                                                                           >  1:1196342/1‑139 (MQ=255)
                      gACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTc                                                                                                          >  1:550427/1‑139 (MQ=255)
                         cGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTcgcg                                                                                                       <  2:142435/139‑1 (MQ=255)
                           tAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGaa                                                                                                     >  2:277216/1‑139 (MQ=255)
                            aGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAAt                                                                                                    >  1:497290/1‑139 (MQ=255)
                                  cGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGaa                                                                                                                           >  2:64282/1‑110 (MQ=255)
                                  cGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGaa                                                                                                                           <  1:64282/110‑1 (MQ=255)
                                  cGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGtt                                                                                              <  1:10164/139‑1 (MQ=255)
                                     agtagCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGt                                                                                           >  2:725176/1‑139 (MQ=255)
                                      gtagCGCAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTg                                                                                          <  2:152551/139‑1 (MQ=255)
                                      gtagCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGt                                                                                           <  1:247290/138‑1 (MQ=255)
                                         gCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGcc                                                                                        <  2:966809/138‑1 (MQ=255)
                                         gCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCa                                                                                       <  1:156225/139‑1 (MQ=255)
                                               cgccgcCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAg                                                                                                            <  2:953279/112‑1 (MQ=255)
                                               cgccgcCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAg                                                                                                            >  1:953279/1‑112 (MQ=255)
                                               cgccgcCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGt                                                                                           <  2:274757/129‑1 (MQ=255)
                                               cgccgcCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGt                                                                                           >  1:274757/1‑129 (MQ=255)
                                               cgccgcCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTg                                                                                          <  2:894903/130‑1 (MQ=255)
                                               cgccgcCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTg                                                                                          >  1:894903/1‑130 (MQ=255)
                                               cgccgcCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATaa                                                                                 >  2:550576/1‑139 (MQ=255)
                                               cgccgcCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATaa                                                                                 >  1:671152/1‑139 (MQ=255)
                                               cgccgcCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATaa                                                                                 >  2:632066/1‑139 (MQ=255)
                                                 ccgccAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAAc                                                                                >  1:1208883/1‑138 (MQ=255)
                                                 ccgccAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAAc                                                                                >  2:626820/1‑138 (MQ=255)
                                                  cgccAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCa                                                                              <  1:277216/139‑1 (MQ=255)
                                                  cgccAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCa                                                                              <  2:431776/139‑1 (MQ=255)
                                                  cgccAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCa                                                                              <  1:632066/139‑1 (MQ=255)
                                                     cAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTg                                                                           <  1:626820/139‑1 (MQ=255)
                                                     cAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTg                                                                           <  2:377158/139‑1 (MQ=255)
                                                     cAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTg                                                                           <  2:852670/139‑1 (MQ=255)
                                                         gagGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGt                                                                       <  2:550427/139‑1 (MQ=255)
                                                                   aTAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGc                                                             >  2:721878/1‑139 (MQ=255)
                                                                    tAGCACCCCTCCTTTTCCTTCTGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCTTTAGTGCCACTGTACCCACTGATGTAAATAATGGcc                                                            <  1:685159/139‑1 (MQ=255)
                                                                    tAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGcc                                                            <  2:784565/139‑1 (MQ=255)
                                                                         ccacGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAg                                                       <  2:1208883/139‑1 (MQ=255)
                                                                          cacGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAg                                                       >  2:749667/1‑138 (MQ=255)
                                                                           acGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGaa                                                                                                                           <  1:1214307/69‑1 (MQ=255)
                                                                           acGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGaa                                                                                                                           >  2:1214307/1‑69 (MQ=255)
                                                                           acGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCa                                                                                       <  2:1161607/105‑1 (MQ=255)
                                                                           acGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCa                                                                                       >  1:1161607/1‑105 (MQ=255)
                                                                           acGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAg                                                       >  1:1210350/1‑137 (MQ=255)
                                                                           acGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCa                                                     <  1:576225/139‑1 (MQ=255)
                                                                                 aTTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTg                                                                                          >  1:697420/1‑96 (MQ=255)
                                                                                 aTTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTg                                                                                          <  2:697420/96‑1 (MQ=255)
                                                                                    aCTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTagag                                             >  1:399407/1‑138 (MQ=255)
                                                                                     cTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCATAGTAGAGaa                                           <  1:595785/139‑1 (MQ=255)
                                                                                        cGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACt                                                                                                                         <  1:550121/58‑1 (MQ=255)
                                                                                        cGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACt                                                                                                                         >  2:550121/1‑58 (MQ=255)
                                                                                        cGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGaaaaa                                        >  2:538884/1‑139 (MQ=255)
                                                                                          gCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAAcc                                      <  2:1210350/139‑1 (MQ=255)
                                                                                          gCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAAcc                                      <  1:538884/139‑1 (MQ=255)
                                                                                           cAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAaataa                                                                 <  2:1017594/111‑1 (MQ=255)
                                                                                           cAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAaataa                                                                 >  1:1017594/1‑111 (MQ=255)
                                                                                             ataatTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTg                                                                                          <  1:1038003/84‑1 (MQ=255)
                                                                                             ataatTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTg                                                                                          >  2:1038003/1‑84 (MQ=255)
                                                                                             ataatTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCaa                                    >  1:1039680/1‑138 (MQ=255)
                                                                                                atTTTAATCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAg                                 >  2:723618/1‑138 (MQ=255)
                                                                                                     aaTCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTagag                            <  2:907229/138‑1 (MQ=255)
                                                                                                     aaTCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTagag                            <  1:659577/138‑1 (MQ=255)
                                                                                                     aaTCCGCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTagag                            <  1:749667/138‑1 (MQ=255)
                                                                                                      aTCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCa                          <  2:1120765/139‑1 (MQ=255)
                                                                                                          gCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCa                                                                              <  2:1044363/83‑1 (MQ=255)
                                                                                                          gCTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCa                                                                              >  1:1044363/1‑83 (MQ=255)
                                                                                                           cTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAg                                                       <  1:965428/105‑1 (MQ=255)
                                                                                                           cTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAg                                                       >  2:965428/1‑105 (MQ=255)
                                                                                                           cTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGcc                     >  2:1098295/1‑139 (MQ=255)
                                                                                                           cTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGcc                     <  2:318683/139‑1 (MQ=255)
                                                                                                           cTTCGGAAAGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGcc                     <  1:494153/139‑1 (MQ=255)
                                                                                                                  aaGAGGAATACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCa              <  2:335312/139‑1 (MQ=255)
                                                                                                                          tACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTc      <  1:1098295/139‑1 (MQ=255)
                                                                                                                          tACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTc      >  2:1043908/1‑139 (MQ=255)
                                                                                                                          tACCGGGACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTc      >  1:1109979/1‑139 (MQ=255)
                                                                                                                               ggACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCa                                                                              >  1:556363/1‑62 (MQ=255)
                                                                                                                               ggACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCa                                                                              <  2:556363/62‑1 (MQ=255)
                                                                                                                               ggACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAaataa                                                                 <  2:426799/75‑1 (MQ=255)
                                                                                                                               ggACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAaataa                                                                 >  1:426799/1‑75 (MQ=255)
                                                                                                                               ggACCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGc  >  1:535565/1‑138 (MQ=255)
                                                                                                                                  tCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGc                                                             >  1:1285434/1‑76 (MQ=255)
                                                                                                                                  tCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGc                                                             <  2:1285434/76‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                      
AGAACGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGC  >  NZ_CP009273/4139892‑4140156

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