Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,895,566 (C)6→7 coding (277/1557 nt) yoaE ← CNNM family cation transport protein YoaE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,895,5601.C98.5% 237.0 / ‑3.1 66coding (283/1557 nt)yoaECNNM family cation transport protein YoaE
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (33/32);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

AAGACGGCGTCAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCC  >  NZ_CP009273/1895427‑1895684
                                                                                                                                      |                                                                                                                            
aagaCGGCGTCAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATccccc                                                                                                                         >  2:887051/1‑139 (MQ=255)
   aCGGCGTCAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATcccccccg                                                                                                                      <  1:864274/139‑3 (MQ=255)
    cGGCGTCAGGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATcccccccga                                                                                                                     <  1:1066883/139‑4 (MQ=255)
    cGGCGTCAAGGATGTCAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATcccccccga                                                                                                                     <  1:1101944/139‑4 (MQ=255)
    cGGCGTCAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATcccccccga                                                                                                                     <  2:685786/139‑4 (MQ=255)
    cGGCGTCAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATcccccccga                                                                                                                     <  1:772874/139‑4 (MQ=255)
    cGGCGTCAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATcccccccg                                                                                                                      >  1:160556/1‑136 (MQ=255)
       cGTCAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATcccccccgaac                                                                                                                   >  1:226124/1‑133 (MQ=255)
                  aCAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAgg                                                                                                        >  1:282107/1‑139 (MQ=255)
                      tCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTAATGCAGTTCGGTTGTTGATTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTcgc                                                                                                    >  1:783182/1‑139 (MQ=255)
                        tgtgtgACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAGCAACATAATCAGGTCGCGt                                                                                                  >  1:22869/1‑139 (MQ=255)
                        tgtgtgACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGt                                                                                                  >  1:428077/1‑139 (MQ=255)
                        tgtgtgACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGt                                                                                                  >  1:796392/1‑139 (MQ=255)
                         gtgtgACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTc                                                                                                 >  1:115680/1‑139 (MQ=255)
                          tgtgACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTc                                                                                                 <  1:1123118/138‑1 (MQ=255)
                          tgtgACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTc                                                                                                 <  1:1136836/138‑1 (MQ=255)
                                    accCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCggagaag                                                                                         >  1:446702/1‑136 (MQ=255)
                                    accCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCggagaag                                                                                         >  1:387923/1‑136 (MQ=255)
                                      ccAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTc                                                                                                      <  1:699069/121‑1 (MQ=255)
                                      ccAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTc                                                                                                      >  2:699069/1‑121 (MQ=255)
                                      ccAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCggagaaggaga                                                                                     >  2:110633/1‑138 (MQ=255)
                                        aGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAAt                                                                                  >  1:653704/1‑139 (MQ=255)
                                               ggCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGAc                                                                           >  2:28088/1‑139 (MQ=255)
                                                gCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGa                                                                            <  2:401748/137‑1 (MQ=255)
                                                gCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACg                                                                          >  2:140056/1‑139 (MQ=255)
                                                gCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACg                                                                          >  1:915952/1‑139 (MQ=255)
                                                            ccGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCgg                                                               <  1:388305/138‑1 (MQ=255)
                                                               tGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTg                                                           <  1:539568/139‑1 (MQ=255)
                                                                ggCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTgg                                                          >  1:336403/1‑139 (MQ=255)
                                                                ggCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTgg                                                          >  1:648461/1‑139 (MQ=255)
                                                                 gCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAacaaca                                                                                                                <  2:1202475/84‑1 (MQ=255)
                                                                 gCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAacaaca                                                                                                                >  1:1202475/1‑84 (MQ=255)
                                                                 gCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCggagaaggag                                                                                      <  2:113341/110‑1 (MQ=255)
                                                                 gCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCggagaaggag                                                                                      >  1:113341/1‑110 (MQ=255)
                                                                   cGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGt                                                         >  2:1262482/1‑137 (MQ=255)
                                                                          aTGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGa                                                 <  1:1262482/138‑1 (MQ=255)
                                                                          aTGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGAc                                                <  1:185695/139‑1 (MQ=255)
                                                                          aTGAGCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGAc                                                <  2:1115119/139‑1 (MQ=255)
                                                                             aTCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGAcca                                              <  2:943895/138‑1 (MQ=255)
                                                                                gcgGTTTTCTAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGAccatcca                                          >  2:410385/1‑139 (MQ=255)
                                                                                gcgGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGAccatcca                                          >  1:1262045/1‑139 (MQ=255)
                                                                                    ttttCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGa                                       >  2:149915/1‑138 (MQ=255)
                                                                                     tttCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGaaa                                     >  2:256155/1‑139 (MQ=255)
                                                                                      ttCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGaaa                                     <  1:9957/138‑1 (MQ=255)
                                                                                              gTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGa                              >  2:684488/1‑137 (MQ=255)
                                                                                               ttCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGAcag                           <  2:160556/139‑1 (MQ=255)
                                                                                                   tGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACagcagc                       <  1:149915/139‑1 (MQ=255)
                                                                                                     cAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGccc                     <  1:28088/139‑1 (MQ=255)
                                                                                                       gTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCa                   <  2:540692/139‑1 (MQ=255)
                                                                                                       gTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCa                   <  2:648461/139‑1 (MQ=255)
                                                                                                       gTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCa                   <  1:110633/139‑1 (MQ=255)
                                                                                                       gTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCa                   <  2:226124/139‑1 (MQ=255)
                                                                                                        ttCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGAc                                                >  2:1042919/1‑109 (MQ=255)
                                                                                                        ttCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGAcagc                          <  2:604637/131‑1 (MQ=255)
                                                                                                        ttCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGAcagc                          >  1:604637/1‑131 (MQ=255)
                                                                                                        ttCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCCAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGAc                                                <  1:1042919/109‑1 (MQ=255)
                                                                                                          cGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCggagaaggag                                                                                      <  2:508479/69‑1 (MQ=255)
                                                                                                          cGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCggagaaggag                                                                                      >  1:508479/1‑69 (MQ=255)
                                                                                                          cGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTc                                                        >  2:276517/1‑99 (MQ=255)
                                                                                                          cGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTc                                                        <  1:276517/99‑1 (MQ=255)
                                                                                                           ggTTGTTGCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACg               <  2:653704/139‑1 (MQ=255)
                                                                                                             ttgttgCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCa             >  1:962146/1‑138 (MQ=255)
                                                                                                                  gCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCggaga                                                                                           >  1:1062653/1‑56 (MQ=255)
                                                                                                                  gCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCggaga                                                                                           <  2:1062653/56‑1 (MQ=255)
                                                                                                                  gCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGAc                                                >  1:102540/1‑99 (MQ=255)
                                                                                                                  gCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGAc                                                <  2:102540/99‑1 (MQ=255)
                                                                                                                  gCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAAt         >  2:302544/1‑138 (MQ=255)
                                                                                                                  gCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAAt         >  2:355460/1‑138 (MQ=255)
                                                                                                                  gCTTTGAACAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAAt         >  2:706036/1‑138 (MQ=255)
                                                                                                                        aaCAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAgcgc   <  1:818229/138‑1 (MQ=255)
                                                                                                                        aaCAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAgcgc   <  1:1251708/138‑1 (MQ=255)
                                                                                                                        aaCAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAgcgc   <  1:302544/138‑1 (MQ=255)
                                                                                                                        aaCAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAgcgc   <  2:336403/138‑1 (MQ=255)
                                                                                                                        aaCAGCAAGAATATCCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGcc  <  2:1262045/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                            
AAGACGGCGTCAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCC  >  NZ_CP009273/1895427‑1895684

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