Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,829,894 (C)6→5 pseudogene (140/162 nt) BW25113_RS26245 → virulence RhuM family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,829,8890C.97.2% 305.1 / 1.8 71pseudogene (135/162 nt)BW25113_RS26245virulence RhuM family protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (34/35);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

GGCTTCAGGAATACCTGGTCTAAGGTTTTGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTGGCGAAAAAAACTGACTATAAACTT  >  NZ_CP009273/3829753‑3830020
                                                                                                                                        |                                                                                                                                   
ggCTTCAGGAATACCTGGTCTAAGGTTTTGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGccc                                                                                                                                   <  2:44447/139‑1 (MQ=255)
ggCTTCAGGAATACCTGGTCTAAGGTTTTGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGcc                                                                                                                                    <  2:246229/138‑1 (MQ=255)
 gCTTCAGGAATACCTGGTCTAAGGTTTTGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGcccc                                                                                                                                  >  1:1012901/1‑139 (MQ=255)
 gCTTCAGGAATACCTGGTCTAAGGTTTTGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGcccc                                                                                                                                  >  2:1206659/1‑139 (MQ=255)
 gCTTCAGGAATACCTGGTCTAAGGTTTTGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGcccc                                                                                                                                  >  2:1257577/1‑139 (MQ=255)
      aGGAATACCTGGTCTAGGGTTTTGTGATGGCCGAGGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCAGCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGccccctgct                                                                                                                             <  2:80220/139‑5 (MQ=255)
      aGGAATACCTGGTCTAAGGTTTTGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGccccctgct                                                                                                                             <  2:1286244/139‑5 (MQ=255)
            aCCTGGTCTAAGGTTTTGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGccccctgct                                                                                                                             <  2:854398/133‑5 (MQ=255)
            aCCTGGTCTAAGGTTTTGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGccccctgct                                                                                                                             >  1:854398/1‑129 (MQ=255)
             ccTGGTCTAAGGTTTTGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGActct                                                                                                                      <  2:891140/138‑1 (MQ=255)
                ggTCTAAGGTTTTGTGATGGCCGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACACATCCGCGATATTCACGCCAGTATGCGCCGTGATTATTTGCAGGTGATTTGCCCCCTGCTGCATCTCTGc                                                                                                                  >  2:381673/1‑139 (MQ=255)
                  tCTAAGGTTTTGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGcc                                                                                                                 >  1:1152075/1‑138 (MQ=255)
                      aGGTTTTGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCaa                                                                                                               <  1:1257797/136‑1 (MQ=255)
                         ttttGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTAt                                                                                                         <  1:564578/139‑1 (MQ=255)
                         ttttGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTAt                                                                                                         <  1:466841/139‑1 (MQ=255)
                         ttttGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTAt                                                                                                         <  1:261002/139‑1 (MQ=255)
                         ttttGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTAt                                                                                                         <  2:1250367/139‑1 (MQ=255)
                         ttttGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTAt                                                                                                         <  2:1012901/139‑1 (MQ=255)
                           ttGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCg                                                                                                       >  1:169198/1‑139 (MQ=255)
                           ttGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCg                                                                                                       >  1:185376/1‑139 (MQ=255)
                           ttGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCg                                                                                                       >  1:984608/1‑139 (MQ=255)
                            tgtgATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCCTGCTGACTCTCAGCCAACTGTATCg                                                                                                       >  1:820851/1‑139 (MQ=255)
                            tgtgATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGc                                                                                                      >  2:1156550/1‑139 (MQ=255)
                              tgATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCa                                                                                                    >  2:1255912/1‑139 (MQ=255)
                              tgATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCa                                                                                                    >  2:424121/1‑139 (MQ=255)
                                aTGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCAtt                                                                                                  <  2:984608/139‑1 (MQ=255)
                                 tGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTa                                                                                                 <  1:1184029/139‑1 (MQ=255)
                                        tGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTg                                                                                            >  1:562990/1‑137 (MQ=255)
                                              gCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGAtgtt                                                                                    >  2:385537/1‑139 (MQ=255)
                                                      aTCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACt                                                                            <  1:804376/139‑1 (MQ=255)
                                                      aTCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACt                                                                            <  2:1085917/139‑1 (MQ=255)
                                                      aTCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACt                                                                            <  2:893996/139‑1 (MQ=255)
                                                      aTCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACt                                                                            <  1:1212191/139‑1 (MQ=255)
                                                       tCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTc                                                                           >  2:456769/1‑139 (MQ=255)
                                                         cGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCa                                                                                                    <  1:757945/112‑1 (MQ=255)
                                                         cGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCa                                                                                                    >  2:757945/1‑112 (MQ=255)
                                                         cGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGc                                                                         >  1:479622/1‑139 (MQ=255)
                                                         cGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGc                                                                         >  1:475759/1‑139 (MQ=255)
                                                         cGCCTGGCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGc                                                                         <  2:230358/139‑1 (MQ=255)
                                                               tCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAg                                                                    >  1:499910/1‑138 (MQ=255)
                                                                     cttctGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACAcc                                                              <  1:1257577/138‑1 (MQ=255)
                                                                     cttctGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACAcc                                                              <  2:159539/138‑1 (MQ=255)
                                                                         tGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGcccc                                                                                                                                  >  1:24848/1‑67 (MQ=255)
                                                                         tGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGcccc                                                                                                                                  <  2:24848/67‑1 (MQ=255)
                                                                             cTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTta                                                      <  1:984026/138‑1 (MQ=255)
                                                                             cTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTta                                                      >  1:679703/1‑138 (MQ=255)
                                                                                gAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGActct                                                                                                                      >  1:977397/1‑71 (MQ=255)
                                                                                gAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGActct                                                                                                                      <  2:977397/71‑1 (MQ=255)
                                                                                gAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACt                                                                            <  1:925287/113‑1 (MQ=255)
                                                                                gAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACt                                                                            >  2:925287/1‑113 (MQ=255)
                                                                                gAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATAc                                                   <  2:893301/138‑1 (MQ=255)
                                                                                 aaCGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGa                                                 <  2:169198/139‑1 (MQ=255)
                                                                                  aCGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGAc                                                <  1:14042/139‑1 (MQ=255)
                                                                                     cATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACt                                                                            <  2:110220/108‑1 (MQ=255)
                                                                                     cATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACt                                                                            >  1:110220/1‑108 (MQ=255)
                                                                                         cgcgATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGActattcta                                         <  1:829369/139‑1 (MQ=255)
                                                                                           cgATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGActattctat                                        <  2:820851/139‑1 (MQ=255)
                                                                                                ttCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCtt                                   >  1:703635/1‑138 (MQ=255)
                                                                                                     gCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGAc                                                >  2:16638/1‑120 (MQ=255)
                                                                                                     gCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGAc                                                <  1:16638/120‑1 (MQ=255)
                                                                                                        aGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTg                          >  2:794368/1‑139 (MQ=255)
                                                                                                        aGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTg                          >  1:194946/1‑139 (MQ=255)
                                                                                                        aGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTg                          >  1:427834/1‑139 (MQ=255)
                                                                                                          tATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTGGc                        >  2:1020226/1‑139 (MQ=255)
                                                                                                              cgcCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCg                                                                                                       >  2:731092/1‑56 (MQ=255)
                                                                                                              cgcCGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCg                                                                                                       <  1:731092/56‑1 (MQ=255)
                                                                                                                 cGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGa                                                                   >  2:737555/1‑89 (MQ=255)
                                                                                                                 cGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGa                                                                   <  1:737555/89‑1 (MQ=255)
                                                                                                                 cGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTGGCGaaaaa                  >  1:1253284/1‑138 (MQ=255)
                                                                                                                 cGTGTTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTGGCGaaaaa                  >  2:1210701/1‑138 (MQ=255)
                                                                                                                   tgtTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTGGCGAAAAAAAc               <  1:424121/139‑1 (MQ=255)
                                                                                                                   tgtTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTGGCGAAAAAAAc               <  1:1255912/139‑1 (MQ=255)
                                                                                                                   tgtTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTGGCGAAAAAAAc               <  2:562990/139‑1 (MQ=255)
                                                                                                                   tgtTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTGGCGAAAAAAAc               >  1:1238884/1‑139 (MQ=255)
                                                                                                                    gtTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGa                                                 >  2:299087/1‑104 (MQ=255)
                                                                                                                    gtTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGa                                                 <  1:299087/104‑1 (MQ=255)
                                                                                                                    gtTTATTTGCAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTGGCGAAAAAAACt              >  1:1035542/1‑139 (MQ=255)
                                                                                                                             cAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTGGCGAAAAAAACTGACTATaaa     >  1:1204327/1‑139 (MQ=255)
                                                                                                                             cAGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTGGCGAAAAAAACTGACTATaaa     >  1:112939/1‑139 (MQ=255)
                                                                                                                              aGGTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTGGCGAAAAAAACTGACTATAAAc    <  1:1020226/139‑1 (MQ=255)
                                                                                                                                gTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTGGCGAAAAAAACTGACTATAAACtt  <  2:1238884/139‑1 (MQ=255)
                                                                                                                                gTGATTTGCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTGGCGAAAAAAACTGACTATAAACtt  <  2:1204327/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                   
GGCTTCAGGAATACCTGGTCTAAGGTTTTGTGATGGACGATGAGCGGCTGAAGAATCCGCCTGTCGGTTCTTCTGTCCTCGAACGCATCCGCGATATTCACGCCAGTATGCGCCGTGTTTATTTGCAGGTGATTTGCCCCCCTGCTGACTCTCAGCCAACTGTATCGCCATTAAGGTGAAGATGTTGTTTTACTCGCCCCAGAACACCCTTTTTTATACGACTATTCTATTTCTTTAACCGCTGGCGAAAAAAACTGACTATAAACTT  >  NZ_CP009273/3829753‑3830020

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