Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,081,683 C→T E511E (GAG→GAA frvR ← putative frv operon regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,081,6830CT92.6% 267.9 / 8.0 81E511E (GAG→GAAfrvRputative frv operon regulatory protein
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base T (36/39);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

TGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGCTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAATATGGTGTTGTGCGACTAATTGCGCA  >  NZ_CP009273/4081551‑4081809
                                                                                                                                    |                                                                                                                              
tGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGt                                                                                                                          >  1:957221/1‑139 (MQ=255)
 gCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTg                                                                                                                           <  2:988663/137‑1 (MQ=255)
 gCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGtt                                                                                                                         >  1:1241665/1‑139 (MQ=255)
  ccAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGtt                                                                                                                         <  1:138391/138‑1 (MQ=255)
  ccAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTc                                                                                                                        >  2:650300/1‑139 (MQ=255)
  ccAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGCTCCTGTTc                                                                                                                        >  1:875052/1‑139 (MQ=255)
   cAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCg                                                                                                                       <  1:1093477/139‑1 (MQ=255)
    aGATAGCTAAATAGTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGc                                                                                                                      <  2:196236/139‑1 (MQ=255)
       tAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTcc                                                                                                                   <  1:1035840/139‑1 (MQ=255)
       tAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTcc                                                                                                                   <  2:707391/139‑1 (MQ=255)
         gCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAg                                                                                                                 <  1:373827/139‑1 (MQ=255)
            aaaTATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCaa                                                                                                              >  2:1186192/1‑139 (MQ=255)
             aaTATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCa                                                                                                               <  1:756198/137‑1 (MQ=255)
                  ttttAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTgg                                                                                                         >  1:485246/1‑138 (MQ=255)
                    ttAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGAt                                                                                                      >  1:988816/1‑139 (MQ=255)
                    ttAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGAt                                                                                                      >  2:1142922/1‑139 (MQ=255)
                           aTGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATggcgagg                                                                                               >  1:749756/1‑139 (MQ=255)
                                 cgcATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGAtt                                                                                          >  1:696856/1‑138 (MQ=255)
                                   cATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGtt                                                                                                                         >  2:658639/1‑105 (MQ=255)
                                   cATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGtt                                                                                                                         <  1:658639/105‑1 (MQ=255)
                                     tCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACaa                                                                                     >  2:119430/1‑139 (MQ=255)
                                       ggcggcggcGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACTTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATc                                                                                   <  2:550448/139‑1 (MQ=255)
                                       ggcggcggcGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAAt                                                                                    <  1:537052/138‑1 (MQ=255)
                                             ggcgCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGttt                                                                             <  2:146172/139‑1 (MQ=255)
                                                  aGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCAcc                                                                        <  1:1103874/139‑1 (MQ=255)
                                                       aTCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCAccttcgc                                                                   >  2:597455/1‑139 (MQ=255)
                                                              cATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTcgcg                                                            <  1:957725/139‑1 (MQ=255)
                                                               aTGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCgg                                                           <  1:863754/139‑1 (MQ=255)
                                                                  gTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAg                                                                                                                 >  1:87830/1‑82 (MQ=255)
                                                                  gTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAg                                                                                                                 <  2:87830/82‑1 (MQ=255)
                                                                  gTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGc                                                         >  2:749622/1‑138 (MQ=255)
                                                                  gTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCg                                                        <  2:1241665/139‑1 (MQ=255)
                                                                  gTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGCTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGc                                                         >  1:828088/1‑138 (MQ=255)
                                                                   ttAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCga                                                       >  1:440183/1‑139 (MQ=255)
                                                                      aTGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCgatga                                                    >  2:236845/1‑139 (MQ=255)
                                                                      aTGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCgatga                                                    >  2:1205954/1‑139 (MQ=255)
                                                                      aTGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCgatga                                                    >  2:1023171/1‑139 (MQ=255)
                                                                       tGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGCTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCgatgat                                                   <  1:885832/139‑1 (MQ=255)
                                                                          cTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCg                                                <  1:250946/139‑1 (MQ=255)
                                                                            tCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGtt                                                                                                                         <  1:791247/64‑1 (MQ=255)
                                                                            tCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGtt                                                                                                                         >  2:791247/1‑64 (MQ=255)
                                                                             cGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCgatgat                                                   <  1:692147/133‑1 (MQ=255)
                                                                             cGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCgatgat                                                   >  2:692147/1‑133 (MQ=255)
                                                                             cGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGCTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTtg                                             >  2:293199/1‑139 (MQ=255)
                                                                                     aCCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCAt                                     <  2:696856/139‑1 (MQ=255)
                                                                                      ccTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATc                                    <  1:1205954/139‑1 (MQ=255)
                                                                                       cTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCg                                                <  1:292686/126‑1 (MQ=255)
                                                                                       cTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCg                                                >  2:292686/1‑126 (MQ=255)
                                                                                          aaCTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcggcg                                >  2:1159442/1‑139 (MQ=255)
                                                                                          aaCTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcggcg                                >  1:706438/1‑139 (MQ=255)
                                                                                          aaCTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcgg                                  <  2:749756/137‑1 (MQ=255)
                                                                                          aaCTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcgg                                  <  1:1023171/137‑1 (MQ=255)
                                                                                          aaCTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcgg                                  <  1:1142922/137‑1 (MQ=255)
                                                                                          aaCTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcgg                                  <  1:650300/137‑1 (MQ=255)
                                                                                          aaCTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGCTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcgg                                  <  1:293199/137‑1 (MQ=255)
                                                                                          aaCTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGCTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcgg                                  <  1:812602/137‑1 (MQ=255)
                                                                                           aCTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcggcgg                               <  2:3340/139‑1 (MQ=255)
                                                                                           aCTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcggcgg                               <  1:236845/139‑1 (MQ=255)
                                                                                               gTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTaa                            >  1:143717/1‑138 (MQ=255)
                                                                                               gTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAat                           <  1:459845/139‑1 (MQ=255)
                                                                                               gTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAat                           <  1:573499/139‑1 (MQ=255)
                                                                                                ttGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAata                          <  2:990615/139‑1 (MQ=255)
                                                                                                 tGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcgg                                  <  1:88160/130‑1 (MQ=255)
                                                                                                 tGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcgg                                  >  2:88160/1‑130 (MQ=255)
                                                                                                 tGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAattt                         >  2:1197680/1‑137 (MQ=255)
                                                                                                    gCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcgg                                  <  1:977484/127‑1 (MQ=255)
                                                                                                    gCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcgg                                  >  2:977484/1‑127 (MQ=255)
                                                                                                    gCGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAATATGgt                      >  2:480105/1‑139 (MQ=255)
                                                                                                     cGAGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAATATGgtg                     <  2:745198/139‑1 (MQ=255)
                                                                                                       aGGGTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTTATATGGtgtt                   >  2:855971/1‑139 (MQ=255)
                                                                                                         ggTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAATATGGtgttgt                 >  1:609510/1‑139 (MQ=255)
                                                                                                         ggTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAATATGGgggtgt                 >  2:269252/1‑139 (MQ=255)
                                                                                                          gTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCg                                                <  2:146511/107‑1 (MQ=255)
                                                                                                          gTAATAAAAAATCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCg                                                >  1:146511/1‑107 (MQ=255)
                                                                                                                    aTCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcggcgg                               >  2:145332/1‑114 (MQ=255)
                                                                                                                    aTCCACGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcggcgg                               <  1:145332/114‑1 (MQ=255)
                                                                                                                        aCGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAATATGGGGGTGTGCGACTAATTGCGca  >  2:139028/1‑139 (MQ=255)
                                                                                                                         cGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcgg                                  >  1:946093/1‑106 (MQ=255)
                                                                                                                         cGAAAACGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATcgg                                  <  2:946093/106‑1 (MQ=255)
                                                                                                                             aaCGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCg                                                >  2:92330/1‑88 (MQ=255)
                                                                                                                             aaCGTCGTTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCg                                                <  1:92330/88‑1 (MQ=255)
                                                                                                                                    |                                                                                                                              
TGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGCTCCTGTTCGCTCCAGCAATGTGGGATGGCGAGGCGATTAACAATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAATATGGTGTTGTGCGACTAATTGCGCA  >  NZ_CP009273/4081551‑4081809

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